LUNAR THERMAL MEASUREMENTS IN CONJUNCTION WITH 
PROJECT APOLLO 

by 

Sydney P. Clark, Jr. 


Final Report 
September 1973 


Grant No. NGR-07-004-039 
Yale University 
New Haven, Connecticut 06520 


(NASA-CR-136035) LUNAB 1HEF.MAL 
heascbehenis in conjunction with 
APCIIC Final Report (Yale Umv.) 
HC $1C. 25 


FP.CJEC1 
161 P 
CSCL 03B 


N74-10769 


Unclas 

G3/30 15736 


The Information presented herein was developed from NASA-funded 
work. Since the report preparation was nut under NASA control, 
all responsibility for the material in this document must neces- 
sarily reside in the author. 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


Tin's grant funded feasibility studies and development work on the 
lunar heat flow experiment (HFE). The first task was performed under a 
previous grant (NSG-400), and It consisted, among other things, of the 
Investigation of a novel method of measuring heat flow which has the 
advantage of not requiring a drilled hole. The method was found not i.o 
be feasible, however, and the requirement that a drill be developed for 
lunar use was established. These early results are incorporated in 
Appendix I for convenience. 

Once the necessity of drilling a hole in the moon became clear, it 
became desirable to develop an in situ method of measuring lunar thermal 
conductivity. The alternative, to measure conductivity on a returned 
core, suffers from the disadvantages that the volume available for inves- 
tigation is much smaller than that sampled by the in situ method and the 
physical properties of the lunar material may be permanently altered by 
the coring operation. The temperature field surrounding a heater with 
the geometry of a cylinder of finite length was therefore investigated. 

The results were presented in an interim report dated January, 1967 
(Appendix II). The calculations reported there formed the basis of the 
design of the downhole part of the HFE, which was developed and fabricated 
by A. D. Little, Inc. 

Ourlng the final six years of the grant's duration the main activity 
was travel and consultation with colleagues associated with planning and 
fabricating the HFE. Places most frequently visited were Lamont-Doherty 
Geological Observatory and A. 0. Little, Inc. 

Probably the best measure of the success of a research program is its 
results. Four HFE's have been flown and two of these are in place on the 



moon and returning data of high quality. They have achieved and even sur- 
passed their design objectives, and they have proven themselves to be 
rugged, reliable instruments. They are described in the Apollo 15 and 
Apollo 17 Preliminary Science Reports, attached as Appendix III. The two 
Instruments that failed to return data were victimized by circumstances 
that were unrelated to the HFE Itself. One was lost as a consequence of 
the abort of the Apollo 13 landing. On the Apollo 16 mission, the HFE 
was successfully emplaced, but it was silenced when an astronaut inad- 
vertently tripped and broke the cable connecting it to ALSEP. 
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ABSTRACT 


A number of problems related to the feasibility of measuring lunar 
heat flow at the lunar surface or in a shallow hole have been investigated 
with the following results. Study of the steady periodic temperatures 
in the lunar material will give unambiguous information about its 
properties only if the surface material alone has an appreciable effect 
on the amplitude and phase of the thermal wave. Layering tends to reduce 
the amplitude of the fluctuation at a given depth. High-amplitude 
fluctuations near a place where the poorly conducting surface layer is 
missing do not penetrate far and pose no difficulty. Large perturbations 
of heat flow may be caused by irregular? :iea in thickness of the surface 
layer, and a number of closely spaced measurements at a given landing 
site will be required to minimize this source of error. The "blanket" 
method of measuring lunar heat flow is not considered feasible because 
of the necessity of very closely matching the local albedo with the 
blanket, and because a blanket with properties such that an easily 
measured gradient free from periodic fluctuations can be set up by the 
lunar flux requires a prohibitively long time to come to thermal 
equilibrium. Conversely a blanket with a suitable time constant will 
yield only a small, seriously disturbed gradient that will be difficult 


to measure. 


1. INTRODUCTION 


\ 


A measurement of lunar heat flow will be interesting for a number of 
scientific reasons. Heat flow gives more direct evidence about the 
internal thermal regime of a planet than any other measurement that can 
be made at the surface. Limits to the total amount of radioactive 
elements in the planet's interior can be set, as well as limits to its 
initial temperature. In the case of the moon, a determination of heat 
flow will help to decide just how "dead" it is, since the source of 
volcanism and mountain building must ultimately be thermal energy, most 
of which is leaked to the surface to appear as heat flow. The small size 
of the moon makes it especially interesting from the thermal point of 
view. Cooling from the surface has affected some 707. of the volume of 
the moon cor-ipared with about 207. of the earth, assuming the two bodies 
are of the same age. As a consequence the relative importance of initial 
heat and radiogenic heat may be very different on the moon as compared 
with the earth, a possibility which makes a comparison of heat flow from 
the two bodies all the more interesting. 

But granting the desirability of a measurement of lunar heat flow, 
a number of obstacles remain in the way,. On the terrest rial land surface 
heat flow is measured in boreholes, mine? or tunnels reaching depths up 
to several thousand feet. Considerable depths are necessary in order to 
avoid disturbances which occur uear the surface. There is no projpect of 
drilling a deep hole in the moon in the foreseeable future, and any 
measurement of heat flow oust be made at the surface or in a shallow hole. 
The temperatures near the lunar surface are in the first place affected 
by the large monthly variation in surface temperature, and secondly by 
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thermal refraction due to variability in the thickness of the lunar 
surface material, which is known to be of very low thermal conductivity 
compared to solid rock. The goal cf the present study is to assess the 
seriousness of these difficulties. 

In the calcula.ions which follow, the assumption that the lunar 
situation can be adequately represented by a linear mooel, i.tf. a model 
in which the thermal properties of the lunar material are treated as 
independent of temperature, is made. This assumption i~ probably very 
wrong for materials near the lunar surface under ambient lunar conditions. 
Temperatures are below the Deoye temperatures of common rock- forming 
minerals, implying a temperature-dependent specific heat. Radiative 
transfer is presumably an important contributor to the thermal conductivity 
of the porous surface material, and it is strongly dependent on temperature. 
Both factors argue for treatment of nonlinear models, but the additional 
complication is hardly warranted in view of the remaining uncertainties 
ir. the details of the properties of the lunar surface material. Thus the 
present study represents a first approach to the problem, aimed more at 
recognizing difficulties than at removing them. 

Four problems are considered in detail in the following sections. 

The first is the cas«_ of one -dimensional steady periodic heat flow ir a 
stratified medium consisting of two layers of differing thermal properties, 
resting on a substratum of infinite thickness which has a third set of 
thermal constants. An exact solution is obtained for the case of 
semisoidally varying surface temperature. 

A second problem again concerns steady periodic temperatures, this 
time in a two-layered medium with the upper layer absent within a circular 



3 . 


region. Numerical results are obtained for this model of a hole in the 
moon’s poorly conducting surface layer. Perturbations of heat flow due 
to variable thickness of the surface layer are investigated under steady- 
state conditions, and finally results are extended to calculations of the 
disturbances associated with the emplacement of a blanket- type thermal 
fluxmeter on the lunar surface. 


/ 



4. 


2. STEADY PERIODIC TEMPERATURES IN A 3-LAYERED MEDIUM 
A. Theory 

Mathematically the problem can be expressed in the following way. 

The region 0 < x < X. contains material with properties K. , , c etc. 

1 * l * 

(see table 1 for notation), the region < x < contains materia., /ith 

properties K^, etc., and the region x > contains material with 

properties K^, etc. At the boundaries X^ and X^ both temperature and 
' T 

thermal flux (= K are continuous, T 0 as x ® and T = A sin ut 
^ x — — — o 

when x = 0, where A is the constant surface amplitude. Within each 

° 2 

\*~ t i gj’ 

region T must satisfy the equation of heat conduction, ■*—. r = — — 

^ t a ot 

This problem is most conveniently solved by the Laplace transform 
method described jy Carslaw and Jaeger (1959). Further details about 
this particular problem are given by Lachenbruch (1959), who obtained the 
solution for the special case x = X^» If we write T for the Laplace 
transform of T, and use subscripts to identify the three regions, we have 
(Lachenbruch, 1959): 


=- Fexp(q 1 x) + Gexp(-q 1 x), 

(1) 

= Hexp(q 2 x) + Jexp(-q 2 x), 

( 2 ) 

*=• Rexp^q^x), 

(3) 


where F, G, H, J, and R are constants, independent of X. The two boundary 
conditions it each interface, x * Xj^ and 2 S = £ 3 * pl us the condition at 
x - 0, provide 5 equations which determine the 5 unknown quantities 
F, G, H, J, and R. 



T able 1. Definitions of Symbols 


T Temperature 
t Time 

x Depth variable 

Depth to base of upper layer 
X ? Depth beneath base of upper layer 

Depth to top of substratum, *=X^ + X^ 
id Angular frequency, * ?.66 x 10 _ 6 sec"* for 1 lui^r day 
Thermal conductivity of the ith layer 
Density of the ith layer 
Heat ceoacity of the ith layer 
Thermal diffusivity of the ith layer, ■ K,/Pi c 
Thermal inertia of the ith layer, * (K^p^c^)i 
p Parameter of the Laplace Transform, T - {*exp(-pt)T dt 

q i ( P/ a i^ 

Q H eat flow 
L Thickness of blanket 

b Subscript denoting properties of blanket 
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(l +« )A 

+(-p 2 e 3 +& 1 f3 3 -f 2 +P 2 p i )exp(q 2 x 2 ) - 

(4) 

g . A ° ,,K * t< ’ lXl ' q3!!3) C(-? 2 e 3 +e 1 P3^ 2 ! -e 1 f> 2 >«-(->2*2> + 

<P+«>>A 

+ (-P 2 P 3 ”P ^3'E 2 ■PiP2 )eXl ' <q 2 X 2 * ■* 

(5) 

h = 2A n^(-q2 X 3‘ C1 2 X 3 )m i 6o-P,P->) 
2 2 

(p +0) )A 

(6) 

2A uiexp (o^X-'q^Xo) 

J - ° 2 3J1 ^W3 ) 

(pW)A 

(7) 

A^xp(-q 1 X 1 -q 2 X 2 +q 3 x 3 ) 

R - 2 2 

(p +<u )A 

(8) 

where 



A « (-p 2 P3-e l P 3 -tf' 2 2 -^2 P l )eXp< ' 2q 2 X 2' 2q l X l )+ 
+(-p 2 P^if > 3" f> 2 2 ' ,p 2 e l )eXp( ‘ 2q l X l )+ 

+0 2 P 3 'e 1 P 3 -? 2 2 -»f2 p i )exp( ‘ 2q 2 x a )+ 

+ <P 2 P'3^ 1 e 3 ^2 2 ' 4f 2 P l ) 

and x 2 " X 3 ’ x i* 



. <HMf 9 *<*« AM*. rp* 
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The transforms of equations (1 )- (3) can h inverted by the contour integration 
method. Lachenbruch (1959) has already shown that the line integrals involved in 
the inversion contribute only to the initial transient state and have nothing to do 
with steady periodic temperatures. Hence for present purposes we need consider 
only the residues at the poles of the transforms, ^cles are located at 
r = -iu.', following Lachenbruch, we assume that the quantity A (equation 9) has no 
zeros in the com*. lex plane. It can be shown that the solution given below does 
indeed satisfy all of the conditions of the problem, which constitutes proof that 
either u has no zeros or that the residues at the resulting poles contribute 
nothing to the steady periodic part of the solution* 

The residues at p = *ic lead to the following expressions for the temperatures. 

Tj * “{exp (V^2a^x)sin(ujt-VV2CI|X)+ 

+ P 2 Cexp (-2Vu/2ci 1 X 1 -2vV2Q£ 2 X 2 +>yu/2a i x)s in (a* -2^^2a^^~lJm/20L^^¥^/2Oi^ x)+ 

+ ex t (-2Va/2a 1 X 1 -2V<V2cx 2 X 2 -^V2a x *)sin (tut+2Vu/2Q 1 X 1 +2 v /u/2(^X 2 3+ 

+ P 3 C exp (-2vV2ctjX 1 +y«»/2a 1 *) s in (ujt-2V<D^o^ 1 -h/u/^a 1 x>+ 

+ exp (-2VV2a i X 1 -.yVia i x)sin(cut+2 v ' a /2a i X 1 V^/2a i x) > 

+ F^[exp (-2 W^^Xj V^O^x) sin (ut -vVi^xH 

+ exp (-2vV^ a 2 X 2 '^ 2a i x ) s in "VtV^Ct^x) 3+ 

+ P 2 P 3 [exp (-4 v 4e/2a i X 1 -2VcV2a 2 X 2 -h/V^Ct 1 *) 8 ln )+ 

+ exp (-4 v V2a 1 X 1 -2/u;/2a 2 X 2 -K/u)/2a 1 x ) s in (tot •^u^a^Vuy^x) j+ 

+ P 2 P 4 [exp ( -2v'uj/2a 1 X 1 -4yuu/2<^X 2 4 s ^/2a 1 x ) s in (t«t - 2 ^ 0 / 20 ^ Vw/2a~x)+ 

+ exp (-2 in ( <ot+2vV2ajX ^ -\/u/2a jX ) Jf 

+ P 3 P 4 [exp (-2V^V2a[x 1 -2Vt/2a 2 X 2 Vu ftoi-f) sin (a* 

+ exp (-2vV23 1 X l -v^a^s in (lofr^/tu/tajXj -2Vu>/2ar 2 X 2 -^Wla^x ) > 


/ 


l 
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+ Pj Cext (- Vu/2a 1 X 1 -4^2c^X 2 vV2Ct l *)»in(wt-»vV2o 1 x) }+ 

+ P 3 [ex i (-Vu/2cr^X 1 + y^a^xJsinCotf-K/^a^x) j+ 

+ P^ [exp (-^Jii./ 2 CC 2%2 "v^>/20(^x)8in(ojt-v l V20i 3 x) ] } (10) 

A 

T 2 “ nT ^2 +P 4 > exp ^ V«/2a 1 X 1 V«/2c^ O^+Xj-x) Jf s in [u* (X^Xj '*) 3+ 

+ p 2 exp [ -W(V2a 1 x x - U-x/la.^ 3*in[t«tWuV’2a 1 x 1 -7u/2a 2 (X 1 -x) ]+ 

+ P 3 expt-V^OjX, ]8in[ujt+v^£^X 1 -v^72^(X 2 +X 3 -x)3+ 

+ P 4 exp[-2vS/2c^X2 ]fl in[u>t Vu^CI^ (X^-x) ] }+ 

+ ^(HP 3 )exp[-^V2^X 1 -^V2a 2 (Xj-x) IfsinCuJt^^a^X^u/ZQ^ (X 1 ~x)+ 

+ P 2 exp [ -2vVZa^X t -2^^75^X 2 ] a in (ufc+vV^O^WuT^ 0^X3 “*> 3+ 

+ P 3 exp[-2v^'2^X 1 ]sin[u*Wi^ ? 2a^X^4v ,, u72^(X 1 -x)}+ 

+ P 4 exp [ -WidjXj ] 8 i n[<ut V u/20l X j VV20j (Xj+Xj-x) ] } (11) 

A J*i 

T 3* ~5~ feX pt-v'«7^X 1 ^/2a 2 X 2 4v^72^ (X 3 -x) ] [ s in (u* 

(X 3 -x)}+ 

+ P 2 exp[-2^u/2^X 1 -2vV2a 2 X 2 ]8irIu«Vw/2a\x l W^20^X 2 WV^(X3'x)> 

+ P 3 exp [ -2Vuv'2a l X t J s in [ ajtV<vtc^X ^ V^C^^-fvW^a^ (X 3 -x ) ]+ 

+ P 4 exp[* 2 «/u)/ 2 C^X 2 ]s io[(ut % / u/2a 1 X l 4^/fa 2 x 2 Vw/ 2 a 3 (X 3 -X) ] ) ( 12 ) 

where 

D » l+2P 2 exp (-2VSj7^X 1 -2v^72o^X 2 )cos (ZjtftefifUu/lc ^X £ )+ 

+ 2P 3 exp (-2V^V r 2a 1 X 1 )cos (27oL/2a,X 1 )+2P 4 exF (-Zjw/ZCtfX^cos 

-■ 2P 2 P 3 exp (-A7u/2a J X 1 -zjv/lctjtyco* (bju/lo^^ y+ 

+ 2P 2 P 4 exp(-2^2^X l -4 v ^7^X 2 )co*(2^72a^X 1 )+ 

+ 2P 3 P 4 exp (-2VV2a 1 x l -2vV2a 2 X 2 )co« (Vu/20^ 2 vV 2 c^X 2 ) 

+ ? 2 2 exp(-WV20^X l -V«/2a^X 2 )+P 3 2 exp (-VV^XjHP^exp (-^a/ZOjXj) (13) 
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and 


V2 

P . f 2 3 3 4e l g 2 4g 2 ~^1^3 

^ 2 p 3 -^ 1 e 2 -^ 2 2 4f 1 p 3 

P -P 2 P3 4 ^lg2~ P 2 2 '^l P 3 

e 2 P 3 -f^iP 2 +p 2 

P „ P 2 f ? 4fi l P 2^2 ~ g l p 3 

4 Vs 4 ’ 3 !^ 2 '*!^ 


(14) 


The temperature is a sinusoidal function of time at all depths. It is useful to 
have the solutions in the form T = A sin (u)t+ ), i.e. in terms of the amplitude 
and phase of the fluctuations. We write A * Vb^+C^ and * -tan'^fB^/Cj), 
and find that the and are given by the following expressions. 

B i “ ex t (Vuy^ctjx) s in (Vu/Sa^x)-^ £ ex i (^ 0 / 20 ^ - Vu/2c^X 2 Vo/ia^x) - 

-exi (-2vV2a^X 1 -2v^X 2 -v^72a[x) ]sin (2^a/2a 1 X 1 +2 A Ai/2o^X 2 -Jw/2a~x')+ 

+P 3 [ exp (-2^/20^-^ V2a lX ) -exp ( -2vW^ 1 X 1 -vV2a 1 x)]sin (lJw/ 201 ^ -yG72a^x)+ 

f 

+2P^exp ( -?vV2a 2 X 2 -vV2a 1 x)cos )s Wa>/2a, x 

-2P 2 P 3 exp (-4 N /tu/2a 1 X 1 -2 A /«y'2a 2 X 2 +VV2a 1 x)cos ) s iny/a/ia^x 

+P 2 P 4^ ex * C-2^V2a l X 1 - V^/2a 2 X 2 4vV2a 1 x) -exi (-2^03/20^ -^/2a^ 2 -vVzOEjX) ]x 
xefn(2v^^X 1 --^72^x)+P 3 P 4 CexF ( - 2^/20:^ - x ) - 
-exp (>2</tu/2a 1 X 1 - 27 ( 1 / 202 -J<u/ 2 (x^i) ]sin (2«/u j/J^X, ~2*/w/i.a 2 X 2 -Jw/la^x') 

-P 2 2 exp (-Vuj/2a 1 X 1 -47m/2a 2 X 2 -KV2a i x)sirKV'.ajx 

-? 3 exp <-V«V'2a 1 X 1 4Vou/2a 1 x) s Wuy^ctj xfP^ 2 exp (- Vu/2 a 2 X 2 -Vu^a^ xJsWu^cijX (15) 
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» ext (Vt4 / 2a^x)coeVu/20!7x +P 2^ exp (" W^/2C^X 2 VV 2 ^* )+ 

+exp (-2^u/2a^X 1 -2V^2a 2 X 2 -vW^C^x) ]cos (2>/2a 1 X 1 +2>/2a 2 X 2 -v^xH 
+P 3 Cexp ( - 2yfy2c^X )+exp (-2yty2a 1 X 1 -^2a 1 x) 3cos (2y£p^X 1 >^i72^x)+ 

+2P 4 exp (-2^"X 2 V«/2a 1 x)cos (2y^X 2 )coa7o/20f 1 x 
+2P 2 P 3 exp (- V<2a 1 X 1 >2yV2a 2 X 2 4y^c> 1 x) cos ) co&J^'la^r. 

+P 2 P 4 Cexp (- 2Vu/2a 1 X l - V^2C^X 2 Vu/Zotj^xHexp -vVZa^)^ 

XC08 (2^1i/2^X^-v^/2a^x)+P 3 P 4 C e *P (- zVoj/Za^^- V V 2a 2 X 2"*’^ 2a l X ^ 

+exp (- 2 yV 2 Gf 1 X i - 2 y^/ 2 a 2 X 2 -Vai/ 2 a 1 x' y ]ccs -Jw/la^x) 

+P 2 2 exp (- 4 vV^X 1 -W 2 C^X 2 Vtt' 72 a^x)co&yt U / 2 a i x 

+P 3 2 exF ( - 4yu^2a 1 X 1 Vw/2a 1 x ) cosVuv^a^x+P^ exp ( - ) c<Wu/ Za^x 

( 16 ) 

B 2 * 0 > 2 +P 4^ ex *' *- ^/2^X 1 Va'/2c^ (X^X-j-x) ]{ sin[«/u)/2a 1 X l +>/2a 2 (X 2 +X 3 -x) ] 

-p^exp (-2*/ua/2d!jX^-2Vui/2c^X 2 )8in[V(ip l '2a^X^-</4/2a 2 (X^~x)3 
-P 3 exp (-2 Vu 7 ^X 1 )s in[ v ^7H^X 1 Vw/2a 2 (X2+X 3 -x) ] 

+P 4 exp(-2Vu/2a 2 X 2 )8in[- v V2a 1 X 1 W2a 2 (X 1 «x) j } 

+ (1+P 3 )exv [ Yl^bjm/Za^ (X^x^fs InCvV 2 ^^ V^2a 2 (X., -x) ]- 
-P 2 ex r <-2^20^- U&afc) efa£</«/2ot 1 X 1 Vu>/2a 2 C* 2 +X 3 -x) ] 

-P 3 exp (-2^a)/2o: 1 X 1 ) s ln[«/5?^X 1 -+V^2a^ CXj^-x ) j 

+P^exp ( - 2^/ 2Ct2X 2 ) 3 in0/<V2O|X ^ **/^/ (X 2 +X 3 ~x) j } ^ 

C 2 - (P 2 +P 4 )ex F [Vw/2a 1 X 1 V^2a 2 v^2 +X 3" x >^ co8 ^ u,/ ' 2a i X l' h ^ to/20£ 2 (X 2 +X 3" x) -‘ 

+P 2 exp (-2^/2a 1 X 1 -2y'V 2 a 2 X 2 ) co«C7a/2at 1 X 1 V'V2c^ <X 1 -X)3 

+P 3 exp (-2y^72a^X 1 )co8 [yu^a^-VuT 2 ^ (X 2 +X 3 “x) ] 

+P 4 exp (- 2 /u)/ 2 a 2 X 2 )ros[«yu 72 a 1 X 1 -H ^'20, (X^x) 3 } 

+ (1+P 3 ) exp [ -s/w/la^ (X 1 -x)]{co8[ v /u/2Qt 1 X 1 -vV2a 2 (Xj-x) 3 

■! P,exp (-2^7^x 1 -V^2a 2 x 2 )co8[>/2a 1 x 1 Vav / 2a 2 (* 2 +x 3‘ x ^ 

+P 3 <xp(-2y/57^ 1 )<»«C^ a l X i^ /ia 2 <lt l‘ x);! 

+P 4 exp (-2VS7^ x 8 JcoeCV^O^Xj V«^20^ (X 2 X 3 -x)3 } • '• • (18) 


( 18 ) 
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B 3 = P x exp [ -Vu/ZaijX t (X 3 -X ) ] { s ink/up/tejX ^^55^ >^203 (X 3 -x)] 

-? 2 ex f (-V«/2a l X 1 )sln[^^X 1 V^2a^X 2 4 v ^i72a^(X 3 -x) ] 

-P 3 exp (-aV^ajXj ) s injVaj^OjXj V2c * 3 (x 3 -x) ] 

+P 4 exp(-2vV2a 2 X 2 )sin[Vj7^X 1 -V^72a^X 2 (* 3 -*) 3 ) (IS) 

C 3 * P 1 exp[-vV2a 1 X 1 ^ Vap^a^ Vu/2a^ (X 3 -x)]{cos[Vtt/2a 1 X 1 Vtt/2a^X 2 (X 3 -x)] 

+P 2 exp (-2 v V2a i X 1 -2vG72^X 2 ) cos [^ 7 ^X 1 4 vG 72 ^X 2 4 v / 572 a^ (Xj-x) 3 
+P 3 exp (-2 v /u/2a 1 X 1 )cos L/up/ta^ -Jw/lafo+Jt 1755^ (X 3 -x) ] 

+P 4 exp (- 2 ^/(iu/ 2 a 2 X 2 )cos [vV^Xj -v^7^(X 3 “x)3 } ( 20 ) 

An alternative way of expressing the solutions in the middle layer and in 
the substratum leads to results which are simpler in appearance. In the first 
case, one may use Lachenbruch ' s (1559) solution for the two-layer problem, with 
amplitude and phase at the surface calculated from equations (17) and (18) at 
x = X^. In the substratum one may use the simple solution for a uniform half 
space (Cars law and Jaeger, 1959, p. 65), with surface amplitude and phase 
calculated from (15) and (20) at x * X 3 . The apparent simplification achieved 
i.n this way proves to be of little value for practical calculation, however. 

A number of terms which are independent of x, such as those on the right side 
of (13), must be evaluated in order to obtain numerical results in the upper 
layer and at the Interfaces. Once this la done it seems simpler to continue 
t jse the three-layer theory rather than evaluating new expressions which 
appear in the two-layer theory, and which differ from those already evaluated. 

The extreme simplicity of the expression for temperature in a homogeneous medium, 
however, makes the alternative procedure more attractive than the use of 
equation ( 12 ) in the substratum. 
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The three>layer theory leads to expressions which are far too cumbersome 
for hand calculation. Numerical results are easily and rapidly obtained by 
a digital computer, however. Use of the exact theory insures that no unwanted 
initial transients affect the results. If finite difference methods are used, 
assurance of freedom from transients is secured only by repeatedly cycling the 
calculation, a procedure which is far more costly in machine time than is 
evaluating the exact theory. 


12 . 


B. Applications 

In order to apply the theory developed above to the lunar surface, 

the parameters of the problem must either be fixed, or their ranges must 

be restricted by estimate or by lunar observations. There are eight 

independent parameters (two thermal constants for each layer plus the 

thickness of the upper two layers), since density and heat capacity always 

occur in the equations as the product £ r and can be considered a single 

parameter. Nevertheless a very large number of permutations of values 

remains, and it is important to fix as many parameters as possible. 

We shall take _c_equal to 0.2 cal/gm °C in all models; this value is 

appropriate to all common silicate materials under lunar surface conditions. 

Fixing £ does not of course reduce the number of parameters unless £ is 

also fixed. Perhaps the best-known lunar parameter is the thermal inertia, 

3, of the surface layer, which is known from infrared temperature measure- 

2 k 

ments during a lunation to be about 0.0023 cal/cm e C sec* (see for 
example, S inton, 1961, p. 411). From this result we take the product 
K £ £ for the lunar surface to oe, nearly enough, 5 x 10 cal /cm “C^sec. 

The very low value of the thermal inertia is the principal evidence that 
the lunar surface is composed of granular material. 

Analysis of radar echoes from the moon leads eventually to a determination 
of the product of density and dielectic constant. Since the latter quantity ^ 

varies little among common silicates, the density may be inferred from 
these results. According to Evans' (1961) summary, material with the 
properties of loose sand would fit the radar data, i.e. a density between 

3 

1 and 2 gn/cm would be expected. On the other hand, the radar reflections 


may originate from a level beneath the optically defined surface. Lower 



surface densities would then be possible ana wculd -j the automatic 
consequences of several postulated models of luna*. surtace structure 
(Hibbs, 1963; Warren, 1963; Hapke, 1964). We shall consider models 

3 

with p ranging from 0.1 to 2.0 gn^cm . Since c and £ are regarded as 
fixed by other considerations, a choice of £ alfco fixes K for the particular 
model of the surface layer. 

We have no direct information about the properties of the subsurface 
layers. We shall assume that the substratum consists of unfractured 
basic rock; appropriate properties are shown in table 2. The intermediate 
layer is presumably made up of rubble, with properties between those of 
the surface layer and the substratum. Three possibilities have been 
considered in order to indicate the effects to be expected from such a 
layer. They do not exhaust the possible range of properties; models 
with the surface layer resting directly on a solid substratum or with an 
infinite thickness of surface material may be considered limiting cases. 

The thermal properties that have been considered in the following numerical 
calculations are collected in table 2. 

It is useful at the outset to recognize two limiting types of 
amplitude - depth relations. In a homogeneous medium the amplitude of 
the temperature oscillation decreases with depth according to the relation 
A ■ Ao exp (- 7® /2a *). The exponential dancing law is obeyed far from 
the lower contact of a thick surface layer of low thermal diffosivity. A 
different extreme is encountered if the density of the material becomes 
very small. The term in the equation of heat conduction containing the 
time derivative then becomes negligible, and the amplitude is found to 
decrease linearly with depth. The numerical results which follow contain 
examples of both types of behavior. 
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Properties of laye r s 


K 

P 

c 

a 

3 

cal/cm sec°C 

gny'cnr 

cal/gm e C 

cta^/sec 

cal/cnr°C 

I. Surface layer. 
1. 2.5xl0* 4 

0.1 

0.2 

1.25xl0* 2 

2.24x10' 

2. 5xl0* 5 

0.5 

0.2 

5.0xl0* 4 

2.24x10' 

3. 2.5xl0 -5 

1.0 

0.2 

1.25xl0* 4 

2.24x10' 

4. 1.25x10 

2.0 

0.2 

3.12xl0* 5 

2.24x10' 

II. Intermediate layer. 

A 1x1 0‘ 3 1.0 

0.2 

5.0xl0~ 3 

1.41x10' 

B IxlO" 3 

2.0 

0.2 

2.5xl0* 3 

2.00x10' 

C 2xl0~ 3 

2.5 

0.2 

4. 0x10* 3 

3.16x10' 

III. Substratum. 

5xl0' 3 

3.0 

0.2 

8.33xl0* 3 

5.48x10' 

IV- Blanket materials 
SI-10 2. 69x10* 7 

m 

0.032 

0.2 

4.20xl0* 5 

4.15x10 

SI-91 4.14x10*^ 

0.120 

0.2 

1.72x10*8 

3.15x10 

-4 

Plastic 1.0x10 

1.3 

0.2 

3.85xl0* 4 

5.10x10 


sec 

3 

3 

3 

3 


2 

2 

2 


2 


5 

5 

3 



14. 


In the calculations Ao was given the value 314°C. This is not 
the amplitude of the temperature fluctuation at the lunar surface, but 
rather is twice the amplitude of the fundamental mode in the Fourier 
analysis of lunar surface temperature given by Sinton (1961). This term 
is more interesting than the higher harmonics because it is about 5 times 
as large and because it penetrates the most deeply. Doubling the amplitude 
gives the total range of temperature directly. 

Some typical results are shown in figs. 1 through 6. The curves of 
amplitude and phase vs. depth have characteristic shapes; the sharp 
drops in the curves as interfaces are approached are particularly note- 
worthy. Study of both amplitude and phase seems to give little more 
information than study of amplitude alone, although any program of 
temperature measurement would automatically yield both quantities. 

The amplitudes decrease exponentially near the tops of layers about 
a meter or more in thickness. The law of decrease is the same as in a 
serai-infinite region, and the thermal diffusivity of the layer can be 
obtained from the damping observed. Where the exponential law is not 
obeyed, the properties of more than one layer are involved and it is 
doubtful whether they can ever be uniquely untangled. In the situations 
where a linear law applies (cf. figs. 1 and 2), the properties of the 
lower layers assume special importance relative to the upper layer in 
which the linear damping occurs. 

In a case in which measurements of temperature cannot be made throughout 
the thickness of a layer, the proximity of an interface could be detected, 
if Indeed one were near. No more than this qualitative result can be 
obtained unless the depth of the interface is also known (c.f. figs. 4, 5, 
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and 6). Temperatures must be measured at the interface in order to determine 
the properties of the underlying layer reliably, and little more than its 
thermal inertia can be deduced unless some penetration of the underlying 
layer is possible. 

It is worthwhile remarking again that the above conclusions are 
correct only if lJestfrization of the conduction equation is valid. This 
will certainly not be true close to the surface, and will only become 
valid at depths where the oscillations in temperature are severely damped. 
This depth is critically dependent on the surface material. In a homogeneous 
region of material 4 of table 2(1), the amplitude reaches 1 degree at a 
depth of 30 cm. In a homogeneous region of solid rock (substratum of 
Table 2) an amplitude of 1 degree occurs at a depth of 450 cm. In both 
cases the surface amplitude was taken to be 314 degrees, as before. The 
presence of layering would reduce those depths. In practice, the linear 
theory will probably be valid if the amplitudes are less than 10 degrees, 
but should be regarded with suspicion in cases of higher amplitudes. 
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. . STEADY PERIODIC TEMPERATURES NEAR A ROl-E IN 
THE SURFACE LAYER 

The poorly conducting lunar surface layer may locally be absent, 
and the substratum of higher conductivity may be exposed to high- 
amplitude fluctuations in temperature at the surface. Damping near such 
outcrops will be comparatively inefficient, and large amplitudes of the 
thermal wave many penetrate the substratum both laterally and vertically. 
We require an estimate of the extent of serious disturbance. 

A simple geometrical model of an outcrop is obtained as follows. 
Imagine first a two-layer structure of the sort described in the last 
section, i.e. a uniform layer with one set of properties separated by a 
plane boundary from a substratum of different properties. We then remove 
a piece of the upper layer having the shape of a right circular cylinder, 
and fill the resulting hole with material of the substratum. The result 
is a cylindrical protuberance on the substratum extending to the original 
plane surface. 

Analytical solutions to heat flow problems in heterogeneous regions 
of this degree of complexity are unknown, and recourse to numerical 
methods must be had. The following calculations were made from the 
simplest form of finite-difference approximation to the equation of heat 
conduction in cylindrical coordinates (see for example Carslaw and 
Jaeger, 1S57, p. 468,470). The program written for the computer took 
account of different conductivities in the two layers, but did not allow 
for different densities and heat capacities. This simplification does 
not affect the qualitative conclusions drawn from the calculations. A 
second sirplification was to assume that the surface temperature was 
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independent of position and varied with time in the manner shown by 
Sir.ton (1961, fig. 3). Actually the amplitude of the variation would 
be smaller in the hole, because of the better connection between the 
surface and the lunar interior there, and the extent of the perturbation 
of amplitudes is therefore slightly overestimated because of neglect 
of this effect. 

Results of the calculations are shown in fig. 7 as contours of equal 
amplitudes. The conductivity of the surface lever is taken to be l/10th 
that of the substratum. It is evident crom : gure that the effect of 

the hole is negligible at a distance from the „ . equal to its diameter, 
and that serious perturbations do not extend further than about half this 
distance. The amplitudes decrease monotonically with depth everywhere, 
as is shown by the fact that no contour can be intersected more than once 
by any vertical line. Thus there is no tendency for high- amplitude 
fluctuations originating in the hole to "run under" the surface layer. 

It may be concluded from these results that the Influence of an outcrop 
on amplitudes does not persist for a distance greater than its diameter. 



4. STEADY-STATE PERTURBATIONS OF FLUX DUE TO 
IRREGULARITIES IN THE THICKNESS OF THE 
SURFACE LAYER 

If the thickness of the poorly conducting surface layer is variable, 
heat tends to be funneled towards thin spots in the layer and away from 
thick spots, a pheiomenon sometir :s termed thermal r fraction. Refraction 
causes the flux observed at the surface to be high where the insulating 
layer is thin and low where it is thick. Some studies of terrestrial heat 
flow have revealed irregularities which may be attributable to thermal 
refraction. Errors arising from this effect may be large in cases where 
the conductivity contrasts are large; a good terrestrial example would 
be near a salt dome in poorly consolidated, f ine-g’-a.i.ned sediments. 

The contrast in conductivity near the lunar surface may exceed a 
factor of 10 (tab’e 2), a contrast that is considerably larger than one 
would expect to encounter on earth. The proportional change in flux 
scales according to the ratio of the conductivity of the substratum to the 
conductivity of the surface layer, and hence large perturbations may be 
expected near the lunar surface. The question was investigated 
quantitatively by studying the steady-state temperature distribution 
around cylindrical protuberances on the interface between an upper poorly 
conducting layer and a better conducting substratum. The problem is 
analogous to the investigation of amplitudes near an outcrop discussed in 
the last section, but with constant surface temperature. The same machine 
program was used, steady-state conditions being achieved by allowing the 
calculation to iterate until the temperatures stopped changing. 

In the application of a steady-state theory to the lunar surface, 
it must be assumed ''hat the periodic transients either have been avoided 
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by measuring heat flow in a sufficiently deep hole or have been removed 
by observing temperatures over at least one cycle and calculating 
undisturbed steady-state means. The results of his section show that 
even if one of these ways of removing transient effects can be followed 
(neither will necessarily b; ;sjy to carry out), perturbations leading 
to erroneous measurements of »unar heat flow may still remain. 

A number of typical results are shown in fig. 8. Cases (d) and (e), 
in which the substratum crops out at the surface, lead to the largest 
perturbations, but such localities are obviously atypical and could easily 
be avoided. The perturbations are greatly reduced if the irregularities 
in the interface are completely buried as in the other cases shown, but 
nevertheless they arc appreciable. Local variations up to about 501 may 
be found in all of the cases examined. The results shown in fig. 8 were 
calculated for a ratio of conductivities of 10; reference to table 2 
shows that this value is, if anything, too low. A conductivity ratio of 
20 wou d lead to perturbations of a factor of 2 or more, depending on 
whether one considers enhancement or reduction of the undisturbed flux. 

In order to be useful, a measurement of lunar heat flow must lead to 
an estimate of mean flux in a region with dimensions measured in kilometers 
which is accurate to better than 201. If the error is much greater than 
this the numbers will have little significance for geophysical or cosmological 
theory. The mean value of 10 fairly ciosely spaced measurements would have 
the required accuracy, assuming that the individual values are disturbed 
by no more than 501 end that the disturbances are normally distributed 
with zero mean value. This latter requirement implies that the probability 
of finding a given positive disturbance oust be the same as finding a 


m 



negative disturbance of the same amount. It is not at all clear that a 
system of randomly distributed, small, buried craters would have this 
property. Furthermore systematic error would invalidate this statistical 
method of achieving accuracy if, for example, all of the measurements 
were made within a large buried crater so that all were affected by a 
negative disturbance. 

An alternative approach is to escape the ncar-surface perturbations 
by drilling deeply enough to make the measurement beneath them. It 
should be possible to do this, because porosity will be eliminated or 
greatly reduced by the weight of overburden, and very large contrasts in 
conductivity will no longer be possible. Considerable depths of penetration 
may be required, however, since disturbed temperatures extend to a distance 
beneath the bottom of the anomalous region roughly equal to its diameter. 
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5. THE BLANKET METHOD OF MEASURING LUNAR HEAT FLOW 

A. Introduction. 

If a sheet of material of known thermal conductivity is placed on 
the lunar surface and allowed to come into thermal equilibrium, the heat 
flow through the surface can be determined from c asuvements of the 
thermal gradient in the sheet A device consisting of a suitable insulating 
material and temperature sensors for the determination of the gradient 
is known as a fluxmeter, or blanket. Such devices have found meteorological 
application in the study of heat exchange between the ground and the 
atmosphere, but they have never been successfully used in the measurement 
of terrestrial heat flow except in thermal regions where the heat flow is 
orders of magnitude higher than normal. 

The extreme simplicity of the blanket method makes it appear attractive 
as a tool for determination of lunar heat flow. Associated difficulties 
seem to outweigh this advantage, however, as is discussed below. 

B. Simple steady-state blanket theory. 

Since one has complete control over the geometry of the blanket, it 
is possible to select a shape that is amenable to simple theoretical 
treatment. A circular disc with aiameter greatly exceeding thickness 
proves to be a convenient choice. An approximate method of treating this 
problem has been suggested privately by A. H. Lachenbruch, and much of 
the following discussion is due to him. 

Consider a half space with zero initial temperature. If, starting 
at _t®0, the temperature of the surface is maintained at a constant value 
/iT within a circle of radius R and zero outside the circle, then beneath 
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tl e center of the circle (Lachenbruch, 1S57) 

T - AT[erfc x/JUat - { X/Vx 2 +R 2 }er f Cv^+R2/y£at ] < 21 ) 

where x is depth. Beneath the center, as the c ;pth approaches zer - *, the 
vertical gradient a,,. roaches 

7T « -At[ 1/R erfc R//4at + l/v'not] (22) 

and the heat flow approaches 

AQ = AT[K/R erfc R/y^at + tl/Jnt'j (23) 

(See table 1 for definitions of symbols.) In the steady state (23) ^educes to 

Aq - K AT/R (24) 


As an illustration of the application of these results, consider 
the case of a blanket placed on the plane lunar surface. The upper 
surface of the blanket is supposed to be at zero, as is the lunar 
surface outside the blanketed area. The assumption that the steady 
periodic transient has somehow been removed is implicit. If the lunar 
flux is everywhere Q, then A_T in (24) becomes equal to QX^/Kb , where 
Xj, is the thickness of the blanket and the subscri r t b denotes blanket 
properties. From (24) we find a perturbation of flux 

AQ/Q - 1-Q.A = Ml (25) 

* RK b 

due to the blanket. This result is approximate first because the 
undisturbed flux £ was used to calculate AT, and secondly because 
AT is assumed constant when in fact it varies with radius in an 
unknown way. The first objection can be overcome by substituting 
0^ for 0 in the expression for AT, calculating the new disturbance, 
and iterating the process until it converges. For example, if 
* 10 and * 50, we find A Q/0 *0.2 and 0t$/Q ■ 0.8. 

Substituting AT^ ■ 0,3 4E f° r AT leads to Q^/Q * 0.84, and a second 
iteration gives Q^/Q * 0.842. The process evidently converges 
rapidly. The second objection mentioned above is inherent in the 
method, since edge effects are neglected. The error is small if 
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£/X b is large enough. 

In order to get a quantitative idea of the meaning of ' large enough 1 ' 
consider a second illustration of the method. A blanket is now supposed 
to be buried so that its upper surface coincides with the initial plane 
surface. The geometry is identical to that shown in figure 8(e). Again 
we assume uniform flux as a first approximation. The thermal gradient 
in the blanket is thin (/K^, and elsewhere it is Q/K; the corresponding 
temperatures at the level of the base of the blanket are QX^/K^ and 

QX^/K respectively. Equation (24) then gv'Jef. 

X.. 

AC/O - l-C^Q * ^[K/Kh-lJ (26) 

in this case. Iteration again may be used to improve the result. This 
problem can also be solved by the finite-difference method used above in 
section 3, and a comparison of the results gives some idea of the range 
of applicability of the approximate method (figure 9). The finite-difference 
calculations agree well with equation (26) for R/X b gr^ iter than about 20, 
but iteration does not improve the agreement. The iteration process becomes 
unstable for equal to 10 or less. It appears that some compensation 

between the errors arising from neglected edge effects and those due to 
other approximations in the derivation of (26) takes place, and the use 
of (26) without iteration appears to give the more reliable results. 

Since the finite-difference calculations are probably not accurate to 
better than 5 per cent, the results given by the simple approach outlined 
here are satisfactory. 

There is a second type of disturbance arising from the presence of 
a blanket on the lunar surface which may be treated exactly by the present 
method. If the albedo of the blanket does not match that of the lunar 
surface, the mean temperature of the top of the blanket will differ from 
the mean surface temperature. The disturbance of flux can be estimated 
directly from (23) and (24). For example, if a blanket 100 cm in radius 
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rests on material of cr nductivity 5 x 10 ^ (material 1(2) of table 2), 

then a difference in temperature of only 0.2°C produces a steady-state 

*6 2 

disturbance in flux of 0.1 x 10 cal/cm sec. Such a disturbance may 
already be intolerably large; it becomes worse if the surface material 
is a better conductor or if the radius of the blanket is reduced to a 
more manageable figure. It will be difficult to measure the mean 
temperature of the lunar surface to better than 1°C, so that a serious 
disturbance due to mismatching albedo may go completely undetected. 

C. Time -dependent problems associated with the blanket method. 

It is convenient to consider separately two causes of time-dependent 
temperatures. One is the steady periodic regime prevailing near the 
lunar surface, and the other is the transient disturbance arising from 
the emplacement of the blanket. The latter has two sources. The 
blanket may not be at the same initial temperature as the lunar surface, 
and after emplacement the establishment of the lunar thermal gradient 
within the blanket changes both its temperature and that of the lunar 
material. The first source of disturbance can be avoided by careful 
planning, but the second cannot. 

Steady periodic temperatures in the blanket were investigated by the 
methods of section 2. The blanket, taken to be S cm thick, was assumed 
to rest on a thick layer having the properties of layer 2 of table 2(1), 
on 50 cm of such material which rested in turn on the substratum of 
table 2 (III), or directly on the substratum. Three kinds of blanket 
materials were considered (table 2 (IV). Two of them, SI-10 and 
SI-91, are ''superinsulators 1 ' developed by Linde for the storage of 
cryogenic fluids. The thermal conductivity of these materials is 
extraordinarily low, as is shown in the table. K third blanket material 
was assumed to have properties corresponding roughly to those of ordinary 
plastics (e.g. bakellte or plexiglass). 
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Amplitudes and c hases of the temperature variations at the bottom 
of the blanket are shown in table 3 for the various combinations of 
blanket materials and assumed lunar configurations. The amplitude- 
depth curve in the blanket has the same shape as the curves for the 
upper layer shown in figure la; that is, the amplitude at the center 
of the blanket exceeds the geometric mean of the surface amplitude 
(314 *C) and the amplitude shown in the table. Clearly only the 
superinsulators are capable of reducing the fluctua or to manageable 
proportions (order of tens of degrees or less) in the lower half of the 
blanket. It is doubtful whether the mean tem t erature can be determined 
in the ' plastic 1 blanket to sufficient accuracy. The situation is made 

worse by the fact that the expected gradient is inversely proportional 

to the conductivity of the blanket. In the superinsulators the expected 

gradient is on the order of l-10 # C/cm, whereas in the "plastic" a 
-2 -3 

gradient of 10 -10 °C/cm seems likely. 

Hence we find that the use of superinsulators is indicated in 
order to eliminate the steady periodic fluctuations nost effectively 
and to raise the mean thermal gradient to an easily measured value. 

But now we must consider the transient associated with blanket emplace- 
ment. We assume that the lateral dimensions of the blanket are great 
compared with its thickness, so that the problem can be treated as one 
of 1-dlmensional heat flow. The blanket, occupying the region -L - x < o, 
is assumed to have Initial temperature To, and thermal properties 
indicated by the subscript b. The lunar material (assumed uniform) hss 
initial temperature mx, where x ^ o equals depth, and unsubscripted 
properties. 

Writing T for the Laplace transform of T, as before, we find 

T^ » To/p + A sinh q^x + B cosh q^x (27) 

and 


T * mx/p + C exp (-qx) 


(28) 



Table 3. Amplitude and phaces at 
base of blanket 5 cm thick. 


Blanket material 

Substratum 
(table 2) 

Amplitude 

°C 

Phase 

SI-10 

1-2, III 

4.5 

-60 

SI-10 

1-2 

4.5 

-60 

SI-10 

III 

0.2 

-60 

SI-91 

1-2, III 

0.1 

-253 

SI-91 

1-2 

0.1 

-252 

SI-91 

III 

0.004 

-252 

Plastic 

1-2, III 

274 

-11 

Plastic 

1-2 

274 

-11 

Plastir 

III 

60 

-39 
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where A, fl, C are constants independent of x. Application of initial 

conditions T. =* T and T = rax, and the conditions of continuity of 
— b - o 

temperature and flux at x = o lead to 


A 

* [KqT Q (cosh q^L-1) - Km cosh q^Lj/pD 

(29) 

B 

* [KqT Q sinh q^L + l^q^T^ - Km sinh q^Lj/pD 

(30) 

C 

= [-K^q^^icosh q^L-l) - Km sinh q^Lj/pD 

(31) 

where 

D 

=-Kq sinh q^L - K^q^ cosh q^L 

(32) 


Conversion of the hyperbolic functions in (29) through (32) to 
exponentials, and expansion of D by the binomial theorem then leads to 
the following expressions for the temperatures 

T b - T ° fl -w7 c ? ' er£c • £2 ^ !2] 


+ M t (M) » erfc - t (M)" erfc 


T 


To 


_£b_ 

^b 


CO 


I 


(M) n [2 erfc (- 


(2nrfl)L 



) 


- + + 
k b 


(33) 


- ierfc(- 


(2nrf2)L 




>: + 


mx 


(34) 


Here M * (p-p^/^+p^) and the other symbols are defined in table 1. 

Equations (33) and (34) are most convenient to use for small values of 
time, but they converge *or all times. Numerical values of the flux in the 
blanket divided by the undisturbed .’.unar flux are shown in figure 10 for 
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blankets of material SI-10 and SI-91 on dust (table 2(1)2) and substratum 
(table 2 (III)). In the most favorable case the flux in the blanket is 
less than 157. of the equilibrium flux after 1 year. This result is 
virtually independent of the To term; it arises mainly from the m term. 

Hence no matter how carefu’ly the initial temperature of the blanket i6 
matched to the mean temperature of the lunar surface, a major disturbance 
is caused by emplacement of the blanket, and it persists for years if 
the blanket is made of superinsulating material. The higher the 
conductivity of the substratum, the longer is the time required to reach 
equilibrium. The ''plastic" blanket, on the other hand, achieves equilibrium 
within a year. 

Thus we see that the two classes of time -dependent temperatures pose 
difficulties that appear to require mutually incompatible sets of blanket 
properties for their solution. In the examples given one must face 
either a large periodic fluctuation throughout the blanket, or a 
prohibitively long time for equilibrium to be established. It does not 
appear that the use of a blanket material with intermediate properties 
would solve the problem. One would then be confronted with both a large 
periodic fluctuation and a long time constant. The thickness of the 
banket affects its thermal behavior in much the same way as its thermal 
diffusivity, so that no escape can be found by changing this parameter. 

A final consideration about the blanket tyie of flux meter 
concerns its contact with the lunar surface. In all of the foregoing 
calculations it has been assumed that there is no contact resistance 
between the blanket and the lunar surface, a situation that is difficult 
to achieve in practice. The effect of uniform contact resistance is to 
reduce the effectiveness with which the periodic fluctuation is damped 
out in the blanket and to increase the time required to equilibrate with 
the lunar surface. Nonuniform contact resistance, which is likely to be 
encountered duo to irregularities on the lunar surface, will in addition 
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A 

I 

f cause thermal refraction within the dimensions of the blanket. This 

will cause the flux in the blanket to differ from point to point, 
necessitating a large number of temperature sensors to give a proper 
mean gradient. Readout is not necessarily complicated by such a 
requirement, since a single readout of many resistance elements in 
series and/or ; arallel to give an appropriate mean value wouli in all 
probability be feasible. 
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6. CONCLUSIONS 

A. Or.e-dimensional steady periodic temperatures. 

A limited amount of information about the thermal properties of a 
layer can be obtained from a study of amplitude or phase of the thermal 
wave as a function of depth, if the effect of other layers is small. The 
latter condition can be recognized by the exponential decrease in amplitude 
with depth. Study of both amplitude and phase gives little or no 
information in addition to that provided by study of amplitude alone. 

When the properties of more than one layer influence the temperatures to 
an important degree, it may be possible to determine t! ^ properties of 
those layers penetrated completely by a hole. Extrapolation beyond the 
deepest observation of temperature is not reliable unless the depth to 
the next interface is accurately known independently. 

B. Propagation of the thermal wave near a hole in the surface layer. 

A hole or thin spot in the surface layer will let high-amplitude 
fluctuations leak into the substratum, where they may propagate laterally 
to some distance. This effect does not appear to be serious, however. 

The amplitudes are essentially unaffected by the presence of the hole a / 

few meters away. 

C. Thermal refraction due tc irregular thickness of the surface layer. 

This steady-state phenomenon is far more serious than the periodic 
disturbance discussed under B. Conditions very probably exist near the 
lunar surface which cause differences in flux of 50% or more because of 
thermal refraction. Such anomalies can be avoided by meesuring heat flow 
at depths below regions causing refraction. Errors due to this effect can 
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largely be removed by taking the means of several closely spaced observations. 

It seems best to try ti take advantage of both techniques, and cO measure. 

temperatures in the deepest holes practicable at several points at a given 

lunar site. 

D. The blanket method of mtasuring lunar heat flow. 

The following difficulties are recognized as standing in the way of a 

measurement of lunar heat flow by a blanket-type fluxmetcr. 

1. The flux is disturbed by thermal refraction due to the presence of the 
blanket. This effect can be kept small by choice of proper geometry 
for the blanket, and the correction is calculable. 

2. The flux is disturbed if the albedr of the blanket does not match that 
of the lunar surface and a difference in mean temperature between the 
blanket ana the surface is thereby created. This disturbance is serious 
if the mismatch in temperature exceeds a few tenths of degrees. 

3. The blanket must be made of poorly conducting material in order to 
damp out the steady periodic temperature fluctuation in a reasonable 
thickness, and also to have a readily raeofedrablc thermal gradient set 
up by the lunar flux. But a blanket satisfying these requirements 
takes years to come into equilibrium with the lunar flux. A blanket 
having a manageable time constant associated with its emplacement 
does not satisfy tne requirements imposed by the steady periodic 
fluctuations and the small value of flux to be measured. 

4. The flux through the blanket may vary from point to point because of 
variable contact resistance with the lunar surface. A large number of 

temperature sensors would be necessary to measure a meaningful average flu; . 
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Difficulties (2) and (3) in particular seem insuperable and make tne 
blanket method unattractive for the measurement of lunar heat flow. 
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Figure Ig. Effect of properties of the surface material on ampli- 
tudes in c layer 20 cm. thick. Numbers beside the curves identity 
the rows in table 2(1). Lower layer is substratum of table 2 (10). 
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b. Effect of properties of the surface material on 
les. Same as figure la except ampiitude scale is iinear. 










Figure 3 a. Effect of properties of the intermediate layer on 
amplitudes. Upper layer 20 cm. thick of material I of table 
2 (I). Curves A, B, C for 200 cm. thickness of materials 
A t B, C of table 2 (II). Curve 0, no intermediate layer, 
surface mcterial rests directly on substratum. 
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Figure 4 a. Effect of properties of the intermediate 
layer on amplitudes. Upper layer 20 cm. thick of 
material 3 of table 2 (I). Curves A, B, C for 200 cm 
of materials A, B, C cf table T (It). Curve 0, no 
intermediate layer, surface material rests directiy on 
substratum. 
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Figure 5. Effect of thickness of the intermediate 
layer on amplitudes. Upper layer 20 cm. thick of 
material 3 of table 2 (I), intermediate layer of 
material A of table 2 (ID. Numbers beside 
curves give thickness of intermediate layer. 
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Figure S. Effect of thickness of the interme 
diate layer on amplitudes. Upper layer 20 
cm. thick of material 3 of table 2 (I), 
intermediate layer of material B of table 2 
(II) . Numbers beside curves give thickness 
of intermediate layer. 
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Figure 8. The effect of c number of conf'guntions of the interface be.weer the surface layer 
(stippled) and the substratum on heat flow. In all cases the protuberance has the shop? of a right 
circular cylinder. The line at the top of each figure shows the ratio of flux at a Qiv^n point to the 
flux at great distance from the irregularity for a conductivity contrast of lu. The configuration of 
Icyers is shown at the bottom of the figures. 









Figure 9. Comparison between the simple method of calculating the disturbance due to the blanket (equation 26) 
and the finite-difference calculation. The points were obtained from the equation with and without iterating the 
calculation and the line shows the finite -difference results. 


Time, years 

Figure 10. Ratio of transient to steady -state fluxes for times up to a 
year after blanket emplacement. Curves ! and 2 for blanket of SI -10 
material, curves 3 and 4 for 31-91 (table 2{il 7 )). Curves l and 3 
for blanket resting on dust {tctic 2 (1)2), end curves 2 and 4 for 
blanket resting on substratum {table 2 (HI)). 
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Introduction 


The measurement of heat flow at a lunar site requires knowledge of both 
the vertical thermal gradient and the local thermal conductivity. The former 
quantity can be measured more or less straight forwardly by a suitably 
instrumented probe emplaced in a drilled hole, but the latter presents 
special complications. In normal determinations of terrestrial heat flow, 
the conductivities of samples cored from the hole are measured in the 
laboratory. It is undesirable, and may even be impossible, to rely solely on 
this technique for lunar heat flow, since the sample may either be destroyed 
or may have its thermal properties seriously altered by the operations of 
collection and return to earth. Hence the determination of thermal con- 
ductivity in situ on the moon is clearly desirable and perhaps essential. 

This report deals with a preliminary study of a method of making this measure- 
ment which utilizes a cylindrical ring source. The results presented here 
form some of the fundamental criteria used in the design of a subsurface 
thermal probe for ALSEP by Arthur D. Little, Inc. 

Theory 

Consider a cylindrical hole of radius R, infinite in length, containing 

a cylindrical probe, also of radius R. Between -Z and Z the probe consists 

of a heater of thermal conductivity k^, density p^, and haat capacity c^. 

For > Z the probe has thermal properties k 2 , p ? , and c^, and there is 

no thermal resistance at z * + Z. The lunar material surrounding the hole 

has thermal properties k^, p^, and c^ and there is contact resistance at 

k n dT 

r » R such that a temperature drop AT occurs, given by AT - 

so-called raotation boundary condition). k Q would be k^ at the outer 
surface of the heater, k^ at the outer surface of the probe, and k^ at the 
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inner surface of the hole. The temperature is initially zero everywhere, 
and heat i9 supplied uniformly over the surface of the heater at rate 
£ for time o < £ < t Q . We must find the temperature as a function of r, 
i, and J. . 

The conditions set forth in the preceding paragraph completely specify 
a boundary value problem in heat conduction, but since they involve both 
radial and axial flow in a heterogeneous medium, they are Intractable 
analytically. The problem was solved by finite differences in the following 
way. Consider intervals in space and time 6r, 6z, 6t» and intergers i, j, 
and k st.ch that z = j6z, £ = k6£, and r * ib£ for i <. 1^ and r = (i - l)6r 

for i > = 1^ + 1. The temperature may be regarded as a function of 

i, J[, and k. 1^ 6r * I^r = R, the radius of the hole. However 
T(I^, j, k) / ^(*2 > k) because of the contact resistance, although 
the two points are only infinitesimally separated in space. On the other 
hand at z = Z <= Jfiz, the temperature is continuous. Since the temperatures 
are symmetric about the axis of the cylinder and also about the plane 

z * o, we need consider only positive values of r and z. 

The equations used in the finite-difference calculation depend on 
the points at which the temperature is to be obtained. Referring to the 
schematic space grid shown in Figure 1, let be the dermal 

diffusivity in region l, the heater,^ be the diffusivity in region 2, etc. 
Also, let * c* n 6t/6r^ and ^ * Cl n 6t/6z^, where n * 1, 2, 3. Then we have 
on the axis 

T(o, j, k + 1) - T(o, J, k) (1 - 4M* - 2M*) + T(l, j, k) • 4 m£ 


+ [T(o, j + 1, k) + T(o, J - 1, k)]M* j / J, n - 1, 2, 


( 1 ) 
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c k. + k„ 

,6t . „6t „ , 1 2 


T(o, J, k + 1) « T(o, J, k) [1 - <4^ + 2—,) ( * - + T(l. J, k) 

6r 6*' p l c l + P 2 C 2 


e * k- + k« c . k ? 


+ T(o, J - 1, k) • 2-^, - 1 


6r 2 p l c l + P 2 C 2 


( 2 ) 


In the interiors of regions 1 and 2 


T(i, j, k + 1) - T(i, J, k) (1 - 2M* - 2M*) + [T(i, j + 1, k) + T(i, j - 1, k)]rf 


n n 


+ C (1 - T(i - 1, j, k) + (1 + i-) T(i + 1, j, k)]M* (3) 


n b 1, 2, o < i < I , J / J. 


and in region 3 


T(i, j, k + 1) = T (i, j, k r - 2^ - 2M^) + [T(i, J + 1, k) -f T(i, J - 1, k)]M^ 
+ C (1 - 2i^)T(i - l, j, k) + (1 + 2i^)T(i + 1, j, k)]^, 
i i I 2 (4) 


21 

- 1 r 

Along the outer skin of the heater and probe, we have, setting f * 'r~~Z~T/it **n 

n l, ~ v 4 * 


2IjH6t 

8 n " (l l - l/4)p n c n 6r * 

T(I l , j, k + 1) - T(I l# j, k) (1 - 21^ - f n - g n ) 

+ [Tdj, J - 1, k)+ T{l v j + 1, k)]M* 

+ T(I 1 - 1, J, k)f n + Tdj + 1, j, k)g n , n - 1, 2, j / J (5) 


4 


and 

T(I| > J, k + 1) 


T(I l , J, k)(l - 


2 (JCj + K^St 2I X - 1 Kj + K, 2 


(p i C l + p 2 c 2 )6r2 Xl " ^ ^ + ^ 


6t 

* 2 
6r 


4I-h6t 2H.6t 

'(1^- 1/4) (p^ + P 2 c 2 )6r ) + T(I 1» J ' U k) 


( p i c i + p 2 c 2 ) 6z 


+ T(I., J + 1, k) 


2r^t 


<P 1 C 1 + P 2 c 2 )6z 


21. - 1 , + - .. 

+ T (I. - 1, J, k) 

1 I t - V4 PjCj + P 2 c 2 6r 2 


41 H5t 

+ T(I 2’ J ’ k) (I 1 - l/4)(p iCl -T p 2 c 2 )6r 


( 6 ) 


At times when the heater Is on, terms accounting for its effect must be added 
to the right sides of (5) in region 1 and (6). We write 

r6 1 


26.2896r 6zJ(I 1 - 1/4) 


( 7 ) 


where the numerical factor includes the conversion from total power input, 

<i, in watts to the units of c.g.s. and calories in which the thermal properties 
were expressed. Then a term q/Pj£^ must be added in (5) and a tern 
+ P 2 £ 2 ) mu8t be added in (6) to account for the heat input. 

Along the wall of the hole in the lunar material we have, setting 


21 , + 1 

1 ..r 


2 IjH6t 


I x + I/* and 8 (Ii + i/4)(p 3 c 3 6r 
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T(I 2 , j, k + 1) . T(! 2 , j, k)(l - 2M^ - f*- g') + [T(I 2 , j + 1, k) + 

+ T(I 2 , j - 1, k)]Mj + T(I 2 + 1, j, k)f + T(l v j, k)g (3) 
Finally, along the junction between the heater and the rest of the probe 


T(i, J, k + 1) * T(i, J, k) [1 


(P 1 C 1 + P 2 C 2 )6r2 (P l C l + P 2 C 2 )6 * 2 


J 


+ T(i, J + 1, k) 


2t^t 


2K,6t 


(P 1 C 1 + P 2 c 2 )6z 


■+ T(i» J-l, k) 


( p i c l + p 2 c 2 ) 6s 


+ [T(i - 1, J, k) (1 - ~>) + T(i + 1, J, k)(l + ~)J 


**1 * K 2 )5t , o < i < (9) 

(P 1 c l + P 2 c 2 )6r 2 

Numerical stability proved to be a serious problem. In the Interiors 
of the three regions, the stability criterion is 

l - 2M* - 2M r > o, n « 1, 2, 3. (10) 

n n 

Depending on the relative thermal properties of probe, heater, and moon, 
a more stringent requirement may occur along the axis i ■ o» since here the 
criterion is 

1 2M 2 - 4M* > o, n « 1, 2 (11) 

n n 

3ut even with (10) and (11) satisfied, instability, which always originated 
at i * 1^ and Ig, was sometimes encountered, particularly for relatively 
large values of H. Imposing the additional constraints that 

1 - 2M* - f -%>’>, n « 1, 2 
n n ’n ’ 
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and 

1 - 2M^ - f' - g ! > o, 

did not remove the difficulty. This instability may result from the fact 
that the space step 6r is effectively halved at i *= 1^ and 1 ^, but the 
matter remains unresolved. The time step, 5^, was simply reduced until the 
calculation became stable. 

A second form of numerical difficulty, which may be termed semistability, 
was also encountered occasionally. Immediately after the heater was turned 
on or off, thus disturbing the system, the calculations oscillated, sometimes 
rather violently. The oscillations were damped, however, and the results 
gradually returned to a smooth trend with further cycles of Iteration. This 
semistability could also be eliminated by reducing 6t, thus approximating 
more closely a smooth input of heat. 

Models 

A number of models of probes and of the lunar material have been 
subjected to numerical analysis. The results are extensive and only the 
more relevant ones have been selected for inclusion here. Thermal properties 
of 3 of the probes are shown in Table 1. The thermal conductivity of 
Probe 1 is too low to be practical from an engineering standpoint, but the 
lunar probe is expected to have properties in the range of Probes 5 and 6. 
Further calculations will be necessary when the final configuration of the 
lunar probe is establiahed and its properties sre measured. 
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Table 1 . Thermal characteristics of probes. 


No. 

C 

1 

5 

6 

Heater 





k, cal/cm sec°C 

0 

3 x ID' 7 

IQ" 4 

io' 3 

p, gm/cm 3 

- 

4 x 10 " 2 

0.5 

0.5 

c, cal/gm°C 

- 

0.2 

0.2 

0.2 

Probe body 





k, cal/cm sec°C 

0 

3 x 10 ‘ 7 

io’ 4 

io’ 3 

3 

p, gm/cm 

- 

4 x 10 ' 2 

0.5 

0.5 

c, cal/gm°C 

- 

0.2 

0,2 

0.2 

Moon models are 

shown in Table 2 . Three 

different 

thermal conductivities 

differing by factors 

of 10 were 

used, and for 

the lower 

conductivities, densities 


and hence diffusivities, differing by factor of 4 were considered. These 
models cover the range of values considered likely fcr material close to the 
lunar surface. The ability of a probe to discriminate between them is thea 
a measure of its suitability. 


Table 2 . 

Thermal models of moon. 



No. 

k, cal/cm sec'C 

3 

p, gca/cm 

c, cal/gm°C 

1 

io - 5 

C .5 

0.2 

2 

io - 5 

2.0 

C.2 

6 

io " 3 

1.6 

0.2 

7 

1 

O 

H 

0.5 

0.2 

8 

i 

o 

H 

2.C 

0.2 


\ 

ilz 



Another parameter entering the calculations Is the contact resistance, 

-12 3 

measured by the quantity H. For purely radiative contact H * 5.5 x 10 ET , 
where E is the emissivity. With blackbody conditions H « 4.4 x 10 ^ at 
200°K which is close to the mean lunar temperature. This is about the lowest 
value that H can attain, and it is an interesting case to consider because 
the probe may be designed to assure purely radiative coupling. H can then 
be calculated with confidence, whereas it otherwise remains an unknown 
parameter the value of which must somehow be extracted from the temperature- 
time curve. The effect of varying H was examined by making some runs with 
it set at 10 times the radiative value. 

The lunar probes are to be about 1.9 cm in diameter. The quantity 6 r 
was taken to be 0.475 cm, which places the probe skin at i * 2, and 62 was 
taken equal to 6 r. This is a rather coarse grid, but no refinement of it 
was made in these preliminary studies. The simulation of a 14-hour lunar 
experiment required over 3C minutes on a 7094 in unfavorable cases, and it 
is not worthwhile to choose smaller space steps (which requires reduction 
of the time step as well to maintain stability) until more than hypothetical 
values of the probe parameters are available. 

The length of the heater was taken equal to its diameter, 1.9 cm. In 
rough design calculations it may be desirable to approximate the probe 
configuration using the exact solution for radial flow from a spherical 
heat source, and the "square" shape chosen for the heater gives the closest 
possible approximation to a sphere. Thus the results of this work may be 
compared directly with those obtained from the spherical approximation. It 
should be noted that in the latter approximation nc account of different 
thermal properties between the body of the probe and the lunar material 


can be taken. 
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Numerical results 

It is helpful at the outset to consider the solution for an infinite 
cylindrical source of heat in an infinite medium. In this case the 
temperatures depend on the thermal conductivity and thermal diffusivity 
of the medium, and on the contact resistance. One could hope that the 
dependence on diffusivity could be removed by heating until the temperatures 
became steady, but with this geometry there is no steady state. The 
temperature of the source continues to rise indefinitely. With a heater of 
finite length a steady state is reached; this was an init '.al reason to 
prefer the geometry adopted here to the "line source" geometry* because 
the possibility exists of eliminating the diffusivity as a fa-tor upon which 
the temperature depends. Another attractive feature of tb present 
configuration is its relatively low power requirement. A line source 
demands a certain amount of power per unit length to produce a given 
temperature rise. Hence a long source requires high power. In the present 
case, it was found that 2 milliwatts input power gave adequate temperature 
rises at the heater, and this value for the heat input was used in all the 
calculations . 

The first calculations were aimed at Investigating the possibility of 

achieving a steady state. Results are shown in. Figure 2. In this figure 

and those following, the temperatures are these of the outer surface of 

the probe. In actual lunar probes the temperature sensors will be located 

on the axis, but the temperature difference between these 2 points is 

Insignificant for present purposes. It is clear from the figure that for 

the lower values of jj; the steady state is not achieved after 14 hours, 

and several days of testing may be required to attain it if K is less 
-4 -3 

than 10 . If K * 10 a fev hours suffice. The probe is evidently 



capable of discriminating between various values of K, particularly if the 
heater is operated at low power levels for a long time. The discrimination 
is best at low X, and heating should last for the order of a day or more 
fcr optimum results. 

Similar curves for the case of 1/10 as much contact resistance are 
shown in Figure 3. The discrimination is somewhat better than in Figure 2, 
and the curves have a different shape. The sharp initial ri^e in 
temperature is much reduced. In Figure 4 * he results for a probe of 
higher conduc* vity are shown; the discrimination is not as good as in 
Figure 2. Clearly the thermal conductivity of the probe should be kept as 
low as possible. 

These results show that it is likely that the temperature rise recorded 
during the lunar experiment will depend or. the 3 quantities R, a, and H. 

Some process of curve fitting must be used to determine their values. This 
may be unsatisfactory s_nce many combinations of parameters may give virtually 
identical results. It is therefore important to try to extract more information 
from the experiment, and an obvious way to do this is to record the temperatures 
at more than one point along the probe. The temperature rise at a point on 
the surface of the probe 8 cm from the center of the heater is shown in 
Figures 5 and 6. Figure 5 is for e probe of unrealistically low K, but it 
shows che large differences in rise time that result from the different moon 
models. Intuitively one would expect the curves to be highly sens' 1 tive to 
a and this is born out by the difference between curves 1 and 2 of Figure 6. 

The rise times are about the same for the casei shown there, in which the 
conductivity of the probe is realistic. 8".t if the moon is a better conducto; 
than the probe discrimination still exists at short tines, although it is 
not well-shown on a plot to the scale of Figure 6. Since this is Just the 
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range of conductivity at which the temperatures at the center of the heater 
lose discrimination, complementary Information can be obtained from the 
second sensor. 

So far we have confined the discussion to times when the heater was 
turned on. But a number of short-term numerical experiments have been done 
in which the heater was turned on for only half the time. The durations of 
tue tests were about % hour. The results were that the appearances of the 
cooling curves were virtually identical to the heating curves, but of course 
inverted and displaced in time. Thus there is no new information to be 
obtained from the cooling curves. On the other hand, following the cooling 
curve in effect constitutes repeating the heating experiment, but without 
the aecesslty of expending heater power. It is always desirable to repeat 
experiments if only to get better statistical control. 

Operations on the moon 

All lunar experiments must wait until drilling disturbances have died 
out near the hole. The thermal gradient will be determined nc^t and then the 
heater will be turned on at lex* power (~ 2 milliwatts). The duration of the 
heating cycle will be determined by tbe conductivity encountered. The heater 
will then be turned off and the cooling curve folltrr^d until ambient conditions 
have essentially reestablished themselves. Then, especially if s high lunar 
thermal conductivity is indicated by this experiment, a second heating period 
will be initiated. e heater pcru.?r will be higher (20 milliwatts or more) 
so that the second sensor, displaced along the probe from the heater, will 
record a readily measured temperature rise. By a process of curve fitting, 
which is not coaq>lelely thought out as yet, the quantities K, a, and H will 
be determined. The first 2 of these automatically yield a value of pc, which 
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can be compared vich Che value measured on returned material to give a 
rough check on the internal confistency of the results. An alternative scheme 
would be to assure that H Is known Independently e.g. by making certain of 
radiative coupling alone; then only K and a need be obtained from the 
temperature curves and the accuracy of the measurements will be increased. 

Conclusions 

1. It appears feasible to measure lunar thermal conductivity using a 

cylindrical source o* 

2. It is desirable to have 2 heating cycles, the first at a power level of 
a few milliwatts and the second at 10 or store times that power. 

3. The duration of each heating will range from a few hours to a few days, 
depending on the lunar conductivity. The use of 2 sensors and 2 power 
levels could materially reduce the amount of heating time required. 

4. There is something to be said for assuring radiative coupling to the moon 
so that the contact resistance can be calculated with confidence. Other* 
wise it represents a third unknown parameter to be determined from the 

r 

temperature curves. Some discrimination of lunar conductivity is lost 
by this procedure, but nevertheless more accurate results will probably 
be obtained. 

The beat way of reducing the lunar data remains to be determined. 
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The grid used in the finite-difference approximation. 




Temperature rise at the center of the heater. Probe 5, radiative coupling. Numbers beside the curves refer to moon models (Table 2). 
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11. Heat-Flow Experiment 

Marcus G. l.angxcth.Jr.** SjtJncy P. Chirk, Jr. * John L. Chine, Jr.* 
Stephen J. Kcilun," and Alfred tl. Wechsler c 


The purpose of the heat-flow experiment is to 
determine the rate of hejt flow from the lunar 
interior by temperature and thermal-property mea- 
surements in the lunar subsurface. Heat loss is 
directly related to the internal lemperatuic and the 
rate of internal heat production; therefore, measure- 
ments of these quantities enable limits to be set on 
long-lived radioisotopic abundances (the chief source 
of interior heating), the internal temperature, and the 
thermal evolution of die Moon. 

Preliminary’ analysis of the data from one heat- 
flow probe indicates that the heat flow from depth 
below die Hadley Rilie site is 3.3 X 10 -6 W/cm* 
(±15 percent). This value is approximately cne-half 
the average heat flow of the tardi. Further analysis 
of data over several lunations is required to demon- 
*trafe that this value is representative of the he it ilow 
at the Hadley Rilie site. Subsurtacc tempe-ature at a 
depth of 1 m below the surface is approximately 
352.4 K at one probe site and 250.7 K at the other 
site. These temperatures are approximately 35 K 
above the mean surface tcmneratuie and indicate that 
the thermal conductivity in the surficia! layer of the 
Moon is highly temperature dependent. Eetwcen 1 
and 1.5 m. the rate of temperature increase as a 
function of depth is 1.75 K/'m (±2 percent) at the 
piobc I site. In situ measurements indicate that the 
thermal conductivity of the regolith increases with 
depth. Thermal-conductivity values between 1.4 X 
10 ^ and 2.5 X 10' ’* W/cm-K were determined: these 
values art a factor of 7 to 10 greater than the values 
of the surface conductivity. Lunar-surface brightness 
temperatures during the first lunar night have been 
0 

•LiiW ttt-Mirny ticot.ssivjl Ob'Cr*atoiy. 

**Yak- Funerary 
c Arthur U. l ink-, li..\ 
tpruH-ipal investigator. 


deduced from temperatures of thermocouples above 
and on the lunar surface. The cooldown history aftei 
sunset suggests that a substantia! increase in conduc- 
tivity occurs at a depth on the order of several 
centimeters. Temperature measurements were also 
recorded during the total eclipse on August 6, 1971. 


EXPERIMENT CONCEPT AND DESiGN 

Heat Budget in the Surface Layer of the Moon 

The temperature and the heat flux at the surface 
of the Moon are determined mainly by the solar 
energy impinging on the surface during onc-haiT of 
the 29.5-day lunation cycle. During the ■ mar day. the 
surluCw temperature rises tv — » nmum. . • is. 
which results in heat flow into the subsurface. Aftei 
lunar sunset, the surface temperature drops to mviy 
100 K, and heat flows out of tiie subsurface and is 
lost by radiation into space. These very large tempem- 
turc excursions, in part, arc a result of '.lie exfrt.ir.civ 
low thermal conductivity and volumetric heat ca- 
pacity of the fine rock powders that mantle most of 
the lunar surface (ref. 1 1-1 ) and. in pari, are a result 
of the very tenuous atmosphere of the Moon. The 
low thermal conductivity of the bulk of the regoiith 
(ref. 11-2) strongly inhibits the flow of energy into 
and out of the subsurface. At a depth of approxi- 
mately 50 cm. the large surface variation of 2S0 K is 
attenuated to a nearly undetectable amplitude. 

At low lunar latitudes, the sutfjee Kftnpctaiu.e. 
averaged ever one lunation, is appioxitr.aiely 220 K. 
This mean surface temperature is deieiminc j by the 
balance of solar energy Ho.ving ,n*> ami ene-gy 
radiated out of the surface dming j complete 
lunation. The mean temper, utue m ihe si bsurfa.a fa; 
decimeter depths) may be higher titan the nnati 
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suffice temperature by a few tens of degrees. As will 
be shown in more detail later in this section, the 
mean subsurface temperature at the Hadley Rille site 
is considerably higher than the mean surface tempera- 
ture. This increase in mean temperature is a result of 
the important role of heat transfer by radiation in the 
fine powders on the surface. The efficiency of 
radiative heat transfer between pat tides varies in 
proportion to the cube of the absolute temperature 
(ref. 1 1-3); consequently, during the lunar day, heal 
flows more readily into the Moon than it flows out 
during the night. The subsurface heat balance, over 
one lunation, requires that a substantial steady-state 
heat loss be maintained in the upper several centi- 
meters of lunar material to eliminate the excess heat 
that penetrates during the day (ref. 11-4). This 
outward heat flow causes a steep downward increase 
in mean temperature that extends to a depth of a few 
decimeters, the depth to which diurnal waves pene- 
trate. 

At depths greater than those to which diurnal 
waves penetrate, the thermal regime is dominated by 
heat flow from the junar interior. This Ilow results 
from high interior temperatures and, in the subsur- 
face, is directly proportional to the increase of 
temperature with depth (the vertical temperature 
gradient) and to the thermal conductivity. These 
quantities are related by the equation 

CM) 

where F is the vertical component of the heat flow, 
k is the thermal conductivity, T is the temperature, 
and r is the depth. The average heat flow of the Earth 
has been determined to be 6.2 X 10"^ W/tm 2 hv 
numerous measurements (ref. 1 1-5). Estimates of the 
lunar heat flow, based on microwave emission from 
the Moon, have ranged from i.OX 10“ h W/cin* (ref. 
li-6) to 3.3 X 10'“ \V /cm* (ref. 1 1-7), or one-sixth 
to onc-half the Earth heat Ilow. Thermal-history 
calculations, based on chronjiiilic- and terrestrial- 
isotope abundances for the Moon (ref. 1 1-8), result in 
heat-flow estimates of l X 10“* to 2.5 X 10"* 
W/cm 2 for tills period in the lunar history. Because of 
the exlicmely low conductivity of the regoliih. even 
lltcse very low heat Hows would result in gnd : ents 
ranging from a few tenths of a degree per meter to a 
few degrees per meter. 


Instrument Design 

The essential measurements for determining heat 
flow arc made by two slender temperature-sensing 
probes that aie placed in prcdrillcd holes in the 
subsurface, spaced approximately 10 m apart. Two 
probes enable two independent measurements of the 
heat flow to be made in order to gain some 
knowledge of the L icral varief of heat Ilow at the 

Hadley Rille site. L_ch probe consists of two nearly 
identical 50-cm-long sections (fig. 1 1-1). Each section 
of each heat-flow probe has two accurate (±0-001 K) 
differential thermometers that measure temperature 
differences between points separated by approxi- 
mately 4"’ and 28 cm. With these thermometers, a 
mee urement (with an accuracy of ±0.05 K) of 
acsoiute temperature at four points on each probe 
section also can be made. 

Additional temperature measurements are pro- 
vided by four thermocouple junctions in the cables 
that connect each probe to the electronics unit. The 
thermocouple junctions are located at distances of 
approximately 0, 0.65, 1.15, and 1.65 m from the 
topmost gradient sensors (fig. 11-1). The reference 
junctiofi for these thermocouples is thermally joined 
to a platinum-resistance reference thcnnome;cr, 
which is mounted on the radiator plate of the 
electronics unit. The temperature measurements ob- 
tained from the heat flow expenment are summarized 
in table 1 1-1. 

The differential thermometers consist of four 
platinum resistance elements wired in a bridge con- 
figuration (fig. 11-2). The bridge is excited by 
successive 2.6-mscc, S-V pulses, first of one polarity 
and then of the other. The output voltage, the 
cxcitatiorf voltage, and the bridge current arc mea- 
sured and used to determine the absolute temperature 
and the temperature difference. The ratio of bridge 
output voltage to excitation voltage and the biiJgc 
resistance are calibrated for 42 diffeunt pairs of 
temperature and tempcraturc-dilfcrcnce values. The 
accuracy of these calibrations is better than the 
.v curacies specified in tahic I l-l. The thermocouples 
are calibrated at four temperature points, and the 
reference bridge is calibrated „t live temperature 
points. 

. Conductivity measurements arc music by means «• 
healers llial sut round each of the eight gradient- 
bridge sensors. The experiment is designed to meastuc 
conductivity in two ranges: a lower range of I X 
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FIGURE l J-l -Heat-flow experiment and dual-purpose hear- 
flow probe, (a) Pcat-tlow-experiment equipment, (b) 
Schematic of heat-flow probe. 


JO -5 to 5 X 1(T* W/crn-K and a higher range of 2 X 

10 - 4 to 4 X 10 -J W/cn;-K. To enable measurements 
in the lov/cr range to be made, a heater is energized at 
0.002 \V and the temperature rise of the underlying 
gradient sensor is recorded as a function of time for a 
period of 36 hr. The temperature rise and the rate of 
temperature rise can be interpreted in to mis of the 
conductivity of the aurrounding lunar material. Mea- 
surements in the higher range of conductivities are 
made by energising the same heater at 0.5 W and 
monitoring the temperature rise ai the ring sensor 10 
cm away for a period of approximately S hr. 

Operation of the Experiment 

During normal operation of the experiment (mod* 1 
1 operation), temperatures of all gradient bridg.'. 
thermocouples, and the reference bridge (as well as 
temperature differences of all gradient bridges) are 
sampled every 7.2 min. When a heater is turned on at 
0.002 \\ to enable measurements to be made in the 
lower conductivity r-ngc, the experiment is said to be 
opouting in mode 2. The mode 3 operation is 
C-. • cd for the measurement of conductivity in the 
h: her range. In this node, tempera ureand tempera- 
ture difference at a seL-ctcti ng bridge are read every 
54 sec. These modes j! operutior; and heater turnon 
arc controlled by commands transmitted from Earth. 

The detection circuitry for measuring tridge volt- 
ages at d thermocouple outputs is contained in a 
housing separate front the Apollo lunar surface 
experiments package (ALSF.P) centicl station (fig. 

1 1 - 1(a)) ami Is designed to compensate for amplifier 
offset and gain change, litis compensation is achieved 
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TABU: 1 1 -I . -Summary of the Hcabl'fow-Txftcriinciit Temperature Measurements 


Thermometer 

Number 

Location 


e thljcrcnee K 
'Accuracy 

ten 

~ fringe 

ipeialthc, K 
Acuiracy 

Gradient bridge* 

1 per section 

Sensors separated by 
47 cm 

■ 

i 0.001 
10.01 

190 to 270 

10.05 

King bridge 

1 per section 

Sensors separated by 
28 cin 

42 

10.002 

190 to 270 

10.05 

Thermocouple 

4 per pro be 

Spaced 65 cm apart in the 
first 2.5 m of the cable 
above the probe 


* * 

70 to 400 

10.07 

Tlrerir.ocoiiple reference 
bridge 

1 per experiment 

Mounted on the radiation 
plate of the 
electrrnics box 

. 

• « 

253 to 363 

10.01 


a Gra<lientbridje .msurw.cnts of lempctature difference are made at 2 sensitivities, v ‘h ; ratio of 10 to I . 


by making all bridge measurements with bipolar 
excitation and by measuring the ratio of the output 
voltage fo the excitation voltage (fig. 11-2). 

EMPLACEMENT OF THE EXPERIMENT 
AT THE HADLEY R1LLE SITE 

Drilling of the holes to emplace the heat-flow 
probes was mote difficult than had been expected. 
The resistant nature of the subsurface at the Hadley 
Killc site prevented penetration to the planned depth 
of 3 m. Instead, at the probe 1 site, the borcstem 
penetrated 1:62 m (fig. 11-3); and, at the probe 2 
site, the boiestem ,.enctiated ap^ioximatcly 1.60 m. 
The configuration of the probe in each hole is shown 
in figure 11-3. An obstruction, which was prob- 
ably a break in the stem at a depth of approxi- 
mately 1 m, prevented probe 2 from passing to 
the bottom of the borcstem. Because of the veiy 
large temp:ratuie differences ovr the upper section, 
which extends above the surface, no valid tempera- 
ture measurements were obtained by the ring and 
gradient biidgcs on the probe 2 upper section during 
most of the lunation cycle. 

Tlic shallow emplacement of the probes resulted in 
five cf the cable thermocouples lying on, or just 
above. tl> .■ lunar surface. Thc'C cable thermocouples 
come into radiative balance with the Ittiui surface 
and space, and the measured tcni|ieratures can be 
interpreted in terms of Itmar-suiface brightness tem- 
peratures. A sixth tltcmn ,-ouplc is in the horvslcm 


projecting above the lunar surface at the probe 1 site. 

When this section was written, surfacc-tcmpcrature 
and subsurf ace-temp rature data had been recorded 
for nearly one and a hi .f iunaiiou cycles. During the 
first lunar noon (August 6). a full eclipse of the Sun 
by the Earth occurred. The thermocouples recorded 
surface-temperature data at 54-scc intervals during 
this eclipse. Six in situ conductivity measurements for 
the low range of values also have been conducted. 
Only three of t'ese measurements are reported 
herein. 

SUBSURFACE TEMPERATURES 

The surface-temperature measurements during the 
lunar night and during the August 6 eclipse indicate 
that the surface layer surrounding the probes has an 
extremely low thermal conductivity. The subsurface 
measurements reveal that the conductivity must 
increase substantially with depth, and vjlucs of 
appioximately 1.5 X lO" 1 W/cm-K arc found at a 
depth of I m. With these values of conductivity, it is 
unlikely that ai.v measurable time variation of tem- 
perature as a result of (lie diurnal cycle existed at 
depths below 50 cm before the hores*eni and the heat 
probe were emplaced. However, after emplacement, 
tlic relatively high thermal conductance of the bore- 
stein and the radiative tiaiisfer along the inside of the 
stem allowed surface-temperature varitions to penc- 
tunc to gtcater depths. After several Innations, a new 
periodic steady-state condition will be established 
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KlClIRE H-J.-Panoi ana of (he h?j|-no»-«\pcrir..enl tirfUocnvo; tile (compovtc . phcootrsphs 
AjiS' 7-1 l$fi and ,\Sl5-S7-< ISIS) irul .» cutaway J/_wiaj of Ilw bcat-piotK- and thermocouple 

Ldljiti. 


arc unci the br-cstcm; )nii:;Rv, . o -sever, the boresiein 
and probes will equi'ibrate te’-va-d temperature* that 
exited in ih; subsurface betore emplacement. By 
applying the theory of the cooling of a cylinder, 
which is discussed in the appendix of inis section, it is 
pos'.ihb to dedu-: die ur.msti-ibed rciupeuiute 
profile at depths where diurnal variations are neg- 
ligible and to obtain some estimate of the thermal 
properties cf the icgolith from the fust few hundred 
hours of cquilic ration. 

The temperature histories of all subsurface ther- 
mometers and the evolution of the profiles of 
teir.neiaiure as a function of depth for probes 1 and 2 
are shown in ligtires ll-t a:td 11-5, resr -ctivcly. AU 
sensors iritiiPy cooled very rapidly, and those 
sensors at depths greater ban 0.7 m continued to 
cool monotonically with time and were still cool 1 .g 
aitf 300 lit. Tlte thermometers at depths less than 
0.7 in responded to the high temperatures of the 
borcitcr.t projecting above live surface The tempera- 
t *e ir the top of the probe 1 borcsiern, which is 
projecting above the surface, was 343 K at htnar 
noon. As sltown in figure 11-5, an nbstntciion 
presented heat -flow probe 2 fiom passing r.: the 
bottom oi tlte hole; consequently, th* platinum- 
resistance incrmomeicrs : n the top sci-thm are nut or. 
sc. ’i. Tw clu.mi; do ‘lion of probe 2 is longer 
ilun :h >t u; pntuc 1. prob. >.ly because of the lower 


conductivity material surrounding the So restem ard 
tlte higher initial heat input that resulted from 
extended drilling. 

Extrapolation, to Err" ,: br!jm VsSues 

To extrapolate the sensor ler.ipetjUuis to equilib- 
rium values, tlte first-order approximation of the 
long-term solution of the cooling -cylinder problem 



is used, where T is the absolute temperature of the 
probe in h. i, the tree cquihb. r.tni temperature of 
the probe, and T a is the ptobe initial temperature; 
and where .5 t and 5. are the therm., I It cat capacities 
per unit length of the inner and outer cylinders, 
respectively, ::t V-secjcnt-E. and / is time in seconds. 
A more complete discussion of the dcriejti.m of 
equal -n (11-2) and delinitions ol the vanahhs can 
be lound in the appendix *i this rectioti. 

If an initial estimate %».' the equilibrium tempera- 
ture is in.ide.ii can be sltown tr-al the error in tile 
c iinutc b is giv n by 


t 
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(1 1-3) 


where r, = T, -T : -7^ and T, and T : 

aw two temperatures selected from the lot g-'.crm 
cooling history at times /, and / 3 . respectively . Then 
the true equilibrium temperature is simply * 6. 
The equilibrium temperature determined in this way 
is independent of the initial-temperature estimate. 
The equilibrium temperatures for all sensors not 
i (Tected by tire diurnal variation are shewn >n table 
11-11. The values also arc plotted as funckms of 
depth in figures 1 1-4 and I i-S.Thc accuracy of these 
equilibrium temperatures is ±0.05 K. 

At tire probe I site, the subsurface temperature, 
which increases regularly with depth, is approxi- 
mately 252.0 K at a depth of SO cm. The .ncrease 
along the lower 60 cm of p.obe I is approximately I 
K. For probe 2, the temperature at a depth cf SO cm 
is approximately 250.5 K. Tire two probe 2 sensors 
that are unaffected by the diurnal variations indicate 
an increase in temperature with depth at a rate 
comparable to the rate detected by the probe I 
rensors. This gradient in ;e npereture is a result of the 
outward flow' of heat from the Moon. 

Equilibrium Temperature Differences Along 
the Heat- Flow Probes 

The gradient and ring bridges enable measurement 


of temperature differences between points 47 cm end 
28 cm apart on the probe with an accuracy of tO.O'J! 
K. wliich is far greatci tlun the accuracy so litC 
absolu tc- tern pera t u rc measurements. An analysis simi- 
lar to that used for the equilibration of die individual 
sensors can be used to extrapolate the temperature 
differences to equilibrium values. By using the ttrst 
term of equation (11-13) in the appendix, the 
equilibrium temperature dilfetence between two 
points on the probes is given by the expression 


AT, - AT- - 
1 *'l 


1--^ 

A 


OM) 


where AT i and AT; are temperature differences 
measured at times f, and f 3 , respectively. Equation 
l!-4 is vilid only for very long times (t > 300 hr); 
consequently. only those differential therr*omerers 
th?t have not been affected by the diurnal variations 
in the first (>w hundred hours of observation can be 
used in this analysis. 

. The only differential thermometers that have not 
beei; affected by the dirrnal variations are those on 
the bottom section of probe I. The calculated 
equilibrium-temperature difference across the gradi- 
ent bridge on the bottom section of probe I is 0.779 
K and. across the ring bridge. 0.4S3 K. These result » 
can be interpreted in terms of the temperature 
difference between adjacent points on the borestem 
wall by taking into account the effect of radiative 


TABLE 1 1 -II.— Lunar Conductivities Determined From /iect-Probc and Cooldown Histories 


TCI IB 
TR22A 
TCI2A 
TG22B 
TR12A 
TRI2B 
TG12B 


Equilibrium j Minim Km-cond'tcrivirr cast J Marimum-eimduetiiity care 

temperature. i Initial temperature. I Deduced conductivity. I initial tcn.perature. | Deduced couduemny. 


251.96 

250.53 

252.28 

250.70 

252.40 

252.87 

253.00 


W{em-X 

0.8 X IO-* 

1.1 

1J 

1.2 
1.2 
15 
21 


Wjcm-1 C 
1.5 x 10-* 


2.9 

3J 


*The boratert at the probe 2 liu was drilled down ;ppmiimilclr I m during the first txttmkkulM activity (KVA'. Pi.'oj 
EVA-2, 22 hr bln, the bureMem was drilled an ajdirionjl V min to a ilrpili uf ap|>rovimslel>' 1-5 ni. the initial li'inpcl ■" . 0 was 
estimated hr catcwblm; I : cooldown lion the pir.-mpbt<n»nl wnpvramnr for 22 hr t.iwwn the KVA pcric.K r.n) rJi'iu; the 
addition— heat of dolling during K.VA-2. this procedure mulled in an olnuK :>f initial temperature (Ji3 K) wry clow n- the 
nlM dew wined hy extrapolation of irmpetatnrc data m the IM hour after insertion. 
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pi o Ik*, nreaxurcd variations wiili a 2 K amplitude. Tlic 


TABLE 1 1 -III. Sunmriry oj f jiiiilibrunn- 
Ttniftcri/tuic lhjjcrciu c McuuitcmcHii 
and Suhmrjui v- 7 cmi/h r,ntirv (inklk'nts 


tX/frrrmtmt 

aunmorseftr 

i'lfmitfrnuiM 

temper* nt*e JiJ/rrrm^. K 

l.vumr 

tcmpi'rmtuif 

graJu'Kt. 

k/m 

Ftobc 

Borrslrm 

Gradient bridge 
R«Hp bridge 

0.779 

.487 

OKI 9 
SOI 

1.74 

1.77 


coupling between tlio walls and the probe ami the 
finite axial conductance of tlic piobc. Tlic tempera- 
ture difference over the probe is always slightly less 
than the temperature difference between tlic adjacent 
points on tlic borestem. The ratio of the two. which 
is called the shorting ratio, was determined experi- 
mentally in the laboratory fur cadi section of the 
probes. The temperature differences at points in tlie 
borestem adjacent to the ditrercntial thermometers 
after the shorting ratio has been applied are listed in 
table I I H1 

Tlie relatively high axial conductance of ihc 
borestem results in some axial shunting of the 
steady-stale heat flow; therefore, to determine the 
undisiurheJ gradient (i.e.. the gradient at large radial 
distances from the borestem), some correction must 
be made. 11 . slum. mg effect can be estimated by 
modeling the bor*>tem as a prolate spheroid sur- 
rounded by a medium with a lower conductivity. By 
using a i effective axial conductivity of 2.3 X 10" J 
W/cm-K. for tlie borestem and a lunar conductivity of 
1.7 X KT 4 W/cni-K, tlie model indicates that a 
plus- 1 -percent correction should be applied to tlic 
borestem temperature gradients. This correction has 
been applied to die temperature gradients listed in 
table I Nil. 

Diurnal Temperature Variations 

Variations in temperature synchronous with the 
solar pliasc were observed at depths as great as 70 cm 
during the first one and a half lunations after 
emplacement. The temperature variations measured 
by probes I and 2 ere shown in figures I 1-6 and ’1-7. 
respectively. Tlic pcak-to-peak amplitude of tlic 
variation at tlic top of probe I is approximately 6 K. 
a 43-to-l attenuation of tlic 2f»0 K temperature 
excuision measured in the pari of the borestem that 
projects above the lunar surface. Tlic ring sensor 
TKI I A. which is located 9 nn below the «op of iIk* 


scnsvMs or. ihc lower section of prolie 2 that detected 
variations are somcwlut deeper (49 and 5S en.)and 
recorded correspondingly smaller amplitude'. There 
are two interesting ..a lures of the obieiveJ varia- 
tions. Tlic jdrase sliift of tlic peaks is extremely small, 
in view of the large attenuation factors, end a 
considerable portion of tiic iugli -frequency com- 
ponent of tlie solar radiation penetrates to these 
depths, a< indicated by tire rapid rales of temperature 
change at dawn and sunsc>. These features sucucri 
tint much of the heat transfer to the probe occurs by 
direct radiative exchange with ihc upper part of tlie 
borestem. Tlic thin, aluminized Mylar disk ilia* is 
located on top of each probe as a radiation slueld 
apparently does not prevent significant radiative 
exclwngc between the top of the probe am! the lunar 
surface. 

Hie heat exdiansc during a lunation cycle is very' 
complex, because the borestem conducts heat from 
and toward tlic lunar surface more efficiently than 
the regothh material. Thus, the low nighttime tem- 
peratures penetrate downward along the borestem. 
which crd.ugcs the low-temperature area viewed by 
the top of the probe and, hence, increases the heat 
loss from the probe to the surface. A similar, bu*. 
opposite, effect occurs during the day. This phenom- 
enon nay. hi part, expbin tire asvmmet.y of the 
plots of temperate.'. as a function of tir.o. 

Diurnal temperature variations that propagate 
along the borestem have an important effect on tire 
mesa temperature in the borestem. Because the 
conductivity of tire borestem is not as temperature 
dependent as tlie adjacent lunar material, heat will ire 
lost more readily along the borestem at night. 
Consequently, tlie heat balance over a full lunation 
will require that the borestem, to depths tlta! diurnal 
variations penetrate, have a lower mean temperature 
at a given depth than the regohth. Thus, a net cooling 
of the borestem in the upper meter can be antici- 
pated, which is an effect that is already apparent in 
the subsianiiai decrease in peak temperature during 
tlic second lunation. Comparisons with the tooling 
curves of deeper sensors show that this difference 
cannot be explained by cooling from initial tcmj'e,-i- 
lures alone. This cooling effect results in a gradient i - r 
mean temperature in the upper meter of the hor-sieni 
that i» unrelated to the I cat flow from the iuterbr. 

It is essential to note that the mean temperatures 
and lc , n r K*i_.urc Jiftetcnccs in those sections of the 
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FIGURE 1 1-6. -Temperature as a function or solir phase anfk foe piobc 1 sensors TR11A and 
1GIIA (the sensors that detect diurnal venations) for the lust one and a hall lunations after 
emplacement. 


bores! cm that see diurnal variations cannot be used to 
determine gradients related to the Ire at flow from tl>e 
!unar interior until the effect of tcmpcralurc- 
depen tien t conduction m the berestem and the 
suricumJir ; lunar material is analyzed and the effect 
qtuntitnlivjly Jctcrmincd. Such an analysis is beyond 
the sco-; of this preliminary report, but the analysis 
will be made oi. future subsurface-icmperai arc data 
after the upper pari of the borcstem has equilibrated 
nearer to the mean periodic steady-state temperature. 
This analysis will add two independent measurements 
of heat flow to the result already reported hete. 

Implications of the Lar.p Mean-Temperature 
Gradients in the Upper 50 cm at the Hadley 
Riiie Site 

By using a finite-difference model to generate 
daytime lunar -surface temperatures (which depend 


almost exclusively on the solar flux) and by using the 
reduced thermocouple tempera l urea to obtain luna- 
*inn nighttime surface temperatures, a mean lunar- 
surface temperature of 217 K (±2 K) was obtained. 
Til's result indicates an increase in mean temperature 
(35 K higher than the mean surface temperatuie) at 
depths beyond which the diurnal variation penetrates. 
This phenomenon can be explained in terms of a 
strong temperature dependence cl the thci trial con- 
ductivity, wtiiwh previously has been iiivwiigated by 
Unsky (ref 1 1 -4) and others. Because of lire near 
lack of an atmosphere on the Moon, radiative transfer 
of heat between and through particles of the lunar 
fines can contribute sicnuiiantly to .Ire effective 
thermal ■ .inductivity. Ibis temperature -dependent 
conductivity lus been found to obey a relation of the 
form (ref. 1 1 -3) 
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FIGURE ll-7.-Tem(wntuie as a function of solar phase ancle for probe 2 sensors TR22A and 
TC22A (the sensors that detect diurnal variations) foe the fust one and a Mf lunations after 
emplacement. 


where k f is the contribution from conduction and 
k f T 3 represents the radiative exchange between and 
through particles. Linskv (ref. 114) has used com- 
puter models of the lunar surface to evaluate this 
effect in the absence of a steady-state heat flow. By 
interpolating from these models, the relative contri- 
butions of the conductive and radiative turns can be 
estimated. For a difference of 35 K in mean 
temperature between the surface and depths at which 
no significant time variations of temperature exist, 
the ratio of radiati'V to conductive touts is approxi- 
mately 2 a: a temperature of 350 K. The relatively 
small steady-stale gradient (1.75 K ,‘m ) ptoduccJ by 
the measured steady -state heat flow will have only a 
slight effect on this tatio. Conductivity measurements 
have been performed for a wide taiii.e of tempera- 
tures on returned lunar samples from the Apollo 11 
and 12 miisiuns (rtf, 1 1-9). The results also indicate 
the significant temperature dependence of conduc- 
tivity. For Apollo 1 1 and 1 2 samples, the data cited 


in reference 11-9 indicateo * the ratios are 0.5 and 
1 .5, respectively. The conductivity of the moie highly 
temperature -dependent lunr: lines from the Apoilo 
12 site seem to be more comparable to the uppet 
regolith conductivity at the Hadley Rillc site Further 
refinement of surface-thermocouple data, combined 
with a more accurate determination of conductivity 
as a function of depth and direct mcasuseiuonts of 
the conductivity of returned Apollo 15 samples, wul 
result in the fiist directly mcasuied profile of regolith 
conductivity to a depth of 1 5 m. 

CONDUCTIVITY OF THF. REGOLITH 

Preliminary Deductions from the 
Heat-Flow-Piobc Histories 

Tlw rale of cquilibratior cf the probes depends on 
the thermal diffusivity k of die sum>u:tding iu..jr 
material, the ratio a of heat capacity per unit volume 
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of the lunar material tc the heat capacity per umt 
volume of Use heat probe, ;;ul the cen'act con- 
ductances //, and Hi , where 

a - 2r.b'ecl(Sy ♦i' J ) 

c = heat capacity per unit inass of the surround- 
ing lunar material, W-sec/cin-K 
p = density of tne surrounding lunar material 
g/cm 3 

ffiJ/j = contact conductances of the inner- and 
oulcr-c> liudcr boundaries, respectively, 
W/cm* -K 

A more complete discussion can be found in the 
appendix of this section. By using an estimate of the 
volumetric heat capacity of the lunar material pc, a 
value for a can be determined, because the thermal 
properties of the prut'' and the borcstemare known. 
From an analysis of the cooling history of the probes, 
an estimate of the diffusivity, and, thus, the conduc- 
tivity, can be made. Measurements of the heat 
capacity of samples that represent a wide range of 
lunar rock tvpes result in very uniform values (ref. 
31-10). The density of the rccoLlh material is quite 
variable; preliminary measure men is of samples taken 
by core tubes at the Hadley Riilc site result in vaiues 
lhai range from 1 JS to i.9! g/cm 3 . At the depth of 
tne probes, the densities are probably ik the high 
end of this range ami not so variable. Foi ti . . v, lysis 
described in this section, a density of 1 2 g . win 3 and a 
heal capacity of 0.66 W-sec/g-K have been assumed. 

It is not possible to determine a value of » from 
the ratios of temperatures at various times during the 
cooldown, because, as the long-term solution indi- 
cates, the temperature ratios •h.ycnd solely on the 
ratio of the times. Bullard ( ref. 11-11) has pointed 
out this property of cooling cylinders in his dis- 
cussion of sea-floor heat -flow measurements. To 
estimate a value for x, the initial temperature must be 
known. E-'i mates of the initial temperature can be 
made by extrapolating data recorded soon after the 
probes were inserted to the time the borcstem was 
emplaced. Tin, estimate is considered to be a mini- 
mum value, because cooling during the first several 
minulcs is faster as a result of enhanced radiative 
transfer at high tempera lures. Alternately, the as- 
sumption can be made that the initial tempcralute s 
the temperature of the borcstem before emplace- 
ment, plus some estimated temperature rise as a result 
of the heat produced during drilling. Temperatures 
recorded before, emplacement by probe 2, wliicii wes 


stored temporarily in the drill rack between LVA-I 
and EVA-2, wci used as estinutes of the borcstem 
tcmpcratuic before emplacement. The temperature 
rise mat resulted from drilling is estimated to be 15 
K/min, based on estimated torque levels. Ihc initial- 
temperature estimates based on these assumptions are 
considered to be maximum estinutes. 

The cooling histories of all subsurface sensors that 
arc not affected by diurnal variations were analyzed 
to determine the conductivity of the surrounding 
luna, material for the two limiting estimates of initial 
temperature. By using the equilibrium temperature 
for each sensor, the ratio (7 - T was 

determined for the lust few hundred hours cf 
equilibration. A typical plot of this ratio as a function 
of time is shown In figure 11-8. The procedure for 
determining conductivity is to make an initial esti- 
mate of the parameters h and A , where 

h = k/bHi 

^ = S i7 k/6 j (S, +Sj) 
yS,PMM t 

a = radius of the inner cylinder (heat-flow probe) 
in centimeters 
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FIGUKL 1 1-8.- Ratio of tern pet a lures measured by senso 
TCil 2 \ during lire inih.J ctpjitihrauon of heal probe 
witli the lunar subsurface compared rvtth the theoretics 
cooldown curve i computed from et|i:alt'm (11-12 
(dailtcd lines). 
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^ ' " heat-flow 

’ j Then, by equaling observed temperature ratios for 

, jevrral times ill the cooling history with the ratios 

; . computed from equation (II 12) m the appendix, a 

^ vaJuc for the dimensionless parameter t can be found 

. that concspouds to each time. Once a value for r is 

. known, k can be determined from the relation t - 

ktfpcb *. Tlic conductivity value determined for a 
t lar^- value of r will be the most accurate, because, at 

j long times. f{A,h, a.r) is neaily independent of h and 

, A, where /E-i, ft. a , t) designates the right-hand side of 

I equation (I I -1 2) in the appendix. By using this initial 

- long time estimate and comparisons of observed and 

: computed ratios at early times, a best value of ft and 

// 2 can be determined with two or three repetitions 
of the procedure. 

Tlieorcticjl curves fitted to the data obtained front 
sensor TC.I2A by this procedure and ratios for both 
limiting ir.ttial-teni|>eraiurc estimates are shown in 
figure 11-8. The parameters /i and A itave been 
chosen to fit data for times greater than 6 to 8 hr. 
The theoretical curve at earlier times lies well above 
the observed data; however, it is not possible to find a 
value for // 3 to fit data obtained at times earlier than 
6 hr without degrading the fit at later times. To 
obtain the most accurate value of Jfc, the curve must 
be fit to the data fo r iarge values of t. 

In table 11-11, die maximum and minimum con- 
ductivity values determined by the procedure are 
shown arranged in order of increasing depth The 
conductivities that were determined for minimum 
and maximum initial-temperature estimates differ, on 
the average, by 50 percent. The more accurate 
conductivity measuieir.enis, which were made by 
using the heaters that surround the gradient sensors, 
resulted in values that lie within the ranges listed in 
tabic ll-lll. The deduced conductivity values arc 
considerably higher than the value obtained from 
measurements on returned lunar fines. The value Tor 
the returned lunar fires is approximately 222 X 10" s 
W/ent-K at 250 *v (ref. 11-12). 11 ic higher con- 
ductivity value ' t wctc obtained may be represent- 
ative of fragn regoiith material in a more dense 
and comprcsseu >ijte than the surface fines. 

In Situ Measurements 

Six in situ conductivity -ncasuicmcnts in mode 2. 
which is the luw-omditctivity mode, were conducted 
at the end of the first lunar night and during the first 
half of the second lunar day. The two hcaicix on the 
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upper section of probe 2 v.ere not tinned on because 

ilie gradient bridges wctc off scale. The inode 2 

measurements indicated the subsurface conductivity 

to be in the lower range of measurement and. ui 

addition, showed that a substantial contact resistance 

exists between the borcstcm and the lunar material. A 

decision was made, therefore, not to run the mode 3 

(lugh-conductivity mode) measurements at tliis time 

because of the possibility that the gradient sensors 

might reach temperaluics potentially dangerous ;o 

the sensor calibration. Mode 3 measurements „re 

planned at some future time after the effects c>f 

heater turnon are examined hv usinc the conduc- 

* 

tivities determined front the inode 2 results. 

Three of the conductivity mcasuicments have been 
analyzed. Two of these measurements were obtained 
by the use of heaters on the lower section of probe 1 , 
the section across which the best temperature gradi- 
ent was uetermineo. The third measurement was 
obtained by the use of the upper heater on the lower 
section of ptobe 2. 

The interpretation of the response of the tempera- 
ture-gradient sensor to heater tumor., in terms of the 
lunar conductivity, is accomplished by using a de- 
tailed finite-difference model (ef. 11*13). A .imple 
analytical model of the gradient jemor lore-ten i (r > 

20 hr) pciforntance deduced from the experimental 
data and die finite-diffcience models will be briefly 
discussed in the following paragraphs. 

The temperature increase as a function of time a! a 
given heater-sensor location upon heater turnon 
depends on the quantity of heat generated and the 
rate at which the generated heat can diffuse outward 
from the heater source. This rate will depend on the 
thermal pioperlics of the material that sut rounds the 
source. The heat will propagate axially along the 
probe and radially from the probe to the drill casing. * 

across the contact-resistance layer outside the casinc. 
and into the lunar medium, itoth radiative tra. sfer 
and conductive transfer are involved in the dissipation 
of heat. Shortly after hcatei turnon, the rate of 
temperature increase at tl.’ gradient sensor will 
depend primarily on the thermal properties of the 
probe arid the borcstcm ui the immediate vicinity of 
Ihc heater and on the resistive gaps between ’.he 
probe and the borcstcm and between the borcstcm 
and the lunar material. As the near-sensor piobc parti 
and the hoiedcm temperature mci cases, a tempera- 
ture drop is established across the resistive gaps. When 
tliis temperature diifcrciicc builds to a relatively large 
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value, heat will flow out from the borestem, across 
the contact-resistance v.ip. .tnd into the medium; am) 
the rate of temperature increase at the sensor will 
level off. At long times (limes giealei than 1000 min 
in this experiment), the temperature increase Ai'(r), 
measured ai the scnsoi , closely fits a relation o f the 
form 

blit) = Cjlnff) ♦ C 2 (11-6) 

where Cj and C 7 are constants that depend on the 
contact conductances H, and f/ 7 and the properties 
of the lunar material. The finite-difference thermal 
model 6f the probe in the lunar material shows the 
same long-term characteristics. This relationship has 
the same form as the lone-term solution to the 
problem rf a uniformly heated infinite cylinder (ref. 
11-14). As in the case of the long-term solution for a 
cylinder, it has been determined from the finite- 
difference models that, at lorg times after heater 
tumon, the constant C, is almost solely a function of 
the conductivity of the surrounding material 3 nd the 
heat input. Tims, the slope !&rfr ; ')- A(f : )]/t|i(/-/fi), 
for times greater than I COO min, is a sensitive 
measure of the conductivity cf the s rounding 
material for a constant heat input. Plots of tempera- 
ture increase as a ‘function of time for the three 
conductivity measurements are shown in figure 11-9 
and compared with best-fitted tneorctica! models. 

Hie magnitude of temperature increase at any 
tine greater than 03 hr after heater turnon is very 
seisitive to the magnitude of the contat. con- 
ductance //,. The value of the contact conductance, 
however, lus no detectable effect on the slope of the 
curve of as a func.ion of ln(r) at long times; 
therefore, the determination of a value for k can be 
made independently of // 2 by matching slopes at 
times greater dun 1000 .u'in. By using this value of k, 
a value for H 7 was dctcnnincd by varying // 2 in the 
finitf-diffe.cncc models until the experimental curves 
of as a function of t ucio bracketed witlun j small 
toleiancc. Examples of models that bracket the 
experimental curves of tie three in situ experiments 
are shown in figure 11-9. 

A rather accurate dctumination of the conduc- 
tivity of the lunar subsurface material that surrounds 
cuJi heater location can he made. Eor example, at 
heater location 1123 (the location of sensor TG22A), 
the long-time slope data could be bracketed by 
moJels of k a I J X 10 * and A = 1.4 X 1 0 “* W/cm- 


K. A linear interpolation between t ,.c models 
resulted in a value k - 1 37(10.02) X i0 ' 4 W,\ii. ”. 
However, the assumption cannot be made that tin. 
models represent the physical situation this accu- 
rately; a value k - l.4(±C.I)X I0 -4 W/cm-K would 
be more realistic. Further examination of die effects 
of the errors introduced by the assumptions about 
the model parameters (probe probities, heat -transfer 
linkages, etc.) must be made so the actual precision of 
the k values can be determined. From previous 
limited parametric studies, a range of ±10 peiccnt 
sire ild represent a maximum bound in the enor of 
the k value determinations. 

The best determinations of conductivity values are 
listed in table 11 -IV. As shown in the table, the 
conductivity determinations from the healer experi- 
ments fall within the range of die k predictions of *!.e 
initial probe cooldown analyses and indicate a sig- 
nificant increase of conductivity w ith depth. 

The contact-conductance valu-.s determined from 
lire in situ measurements vary. The contact con- 
ductance probably corresponds to a thin zone around 
the borestem that is filled with mnar fines. If the 
assumption is made that these fines have a conduc- 
tivity oi 2 X I0 -S W/cm-K, which is similar to the 
conductivity of the surface fines, the widths of the 
disturbed zones would be 2.7, 2.0, and 13 mm for 
the locations of sensors TG22A TGl2A.and TG12B. 
respectively. The larger value of // 2 at the location of 
sensor TG12B might result from greater compaction 
of ihe fines, rather than a tliinnei zone. The thicker 
disrupted zone around probe 2 may have resulted 
from the longer period of drilling. 

STEADY-STATE HEAT FLOW PROM THc 
LUNAR INTERIOR BELOW THE 
HADLEY RILLE SITE 

The conductivity of the regolith is shown by the 
measurements returned from the experiment to be 
significantly variable with depth over the lower 
section of probe I . To compute the heat flow from 
temperature di'feicnces over a finite depth interval, 
he thermal resistance must be known. The thermal 
resistance can be calculated from die equation 



where x ( and x 2 are the end points of the Interval. 
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FIGURE 1 1-9. - Tempcratur . increase as a function of time after heater iumon for heaters located at 
savors TG22A, TG12A, and TG12B. Two computed models that closely bracket the data ate 
shown for each of the three heater locations. Temperature increase as a function of time after ICOO 
min is shown in the inset on an expanded scale. The solid lines in the inset arc best fitting 
computer models. 


TABI.t 1 1 IV. -Conductivity Determinations From 
in Situ ExiKrimcuts 


Healer 

sensor 

Depth. 

cm 

Thcrrul conductivity, 
k. W/cm-K 

Cont.t i conducunce 
of It,' borcsterr., 
M,. H'/fm 1 -K 

TG22A 

49 

1.40 HO. 14) X 10~* 

0.7 X 10^ 

TC12A 

91 

1.70 C 0.17V 

1.0 

TGI211 

13S 

2.501- 0.25) 

1.5 


Thus, for the flux to be determined, k -value variation 
in llie interval between the gradient sensors, which 
are located at depths of '21 and 138 cm, must be 
known. Accurate measurements of k were made only 
at the end points; however, a consiuim can be 


applied on the variation with depth frorr the ratio of 
the temperature differences measured by the ring 
bridge and the gradient bridge. If the beat How is 
uniform with dei 'It, the constraint re |uirc;l by the 
ratio of temperature differences is 



Three possible conductivity profiles are shown hi 




• r* 


figure 11-10. Profile B is based on the trend of 
conductivities from the cooldown cuivcs and obeys 
the constraint of equation (I I -8). Profile A also 
obeys the constraint of equation (11-8) hut includes a 
uniform conductivity or 1.7 W/cm-K to a depth of 
136 cm and. then, a thin layer with a conductivity of 
25 X 10 W/cin-K in which the bottom sensor is 
embedded. Profile A would result in a lower limit for 
the heat flow; profile C indicates a uniform increase 
in conductivity over the probe section. Profile C docs 
not obey the constraint of equation (11-8), but 
defines an upper limit for the neat-flow value. The 
trend of conductivity up to a depth of 50 cm that is 
indicated by the probe 2 measurement makes eases 
with higher conductivity than shown in profile C 
unreasonable. Based on these three profiles shown in 
figure 11-10, the temperature difference over the 
lower section of probe 1 results in the heat-flow 
values listed in table ll-V. The uncertainty of the 
conductivity measurements (*I0 percent) should be 
considered as error bounds on each of the heat-flow 
values listed in table I I-V. 

Analysis of data obtained during a full year will 
enable the previous determinations to be refined 
considerably. In addition, a comparison of the value 
obtained from the bottom section of probe I with 
the values obtained floin the upper section of probe I 
and the lower section of probe 2 can be nude once an 
analysis of the effects of the diurnal variations 'tas 
been completed. 


® Values determined by the 
in situ experiment 
* Values deduced Irom 
cooling curves with 
minimum estimates of 
initial temperature 
O Values deduced from 
cooling curves with 
maximum estimates of 
initial temperatu'e 
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Conductivity, \Vfcm-*K 

MCURI-; 1 1-10. -Conductivity as a function of doth, with 
throe ix.v.iHc conductivity profiles (A. B. amt Q (table 
ll-Vj. 
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ased on the trend of TABLE I l-V.- Hcjt-Ftow Data Obtained 



Profile 

tlcut flow, 
IV/ an* 

Crmmnii 


A 

2.99 X I0- 4 

Lower limit 

~ 

B 

3.31 

Best value 


C 

3.59 

Upper limit 



SURFACE TEMPERATURES DEDUCED 
FROM THE CABLE THERMOCOUPLES 

Of the eight thermocouples designed to measure 
the temperature profile in the upper 1.5 in of the 
heat-flow borehole, six presently measure tempera- 
tures that may be used to deduce the variation of 
lunar-surface brightness temperature tlnoughout the 
lunation period. Of particular interest is the determi- 
nation of temperature during total eclipses and lunar 
nights, which arc measurement difficult to obtain by 
Faith-based telescopic observation. The thermo- 
couples in the cable of the heat-flow experiment, 
lying on or just above the lunar surface, provide u 
means by which these measurements can be obtained 
at a sampling rate previously unattainable (one 
measurement set each 54 see). 

During the lunar night, the thermocouples come 
into radiative balance with the lunar surface and 
space. To determine the relationship between the 
cable temperature and the lunar-surface brightness 
temperature, the heat balance for a small cylindrical 
cable element of radius a and length til can be 
considered. The heat balance for such an element 
arbitrarily oriented above the lunar surf.ee during the 
lunar night can be represented by 

F C - y 2Wh rlk T* - 2*oJU c , T* - VfC . 0 (11-9) 

where the first term is the energy received by the 
cable element from the Moon per tuit time, the 
second term is the energy lost from the cable element 
per unit time, and the lltird teim is the energy 
rcqmrcd to change the temperature of the cable 
element per unit time; and where 

F c-m m v ’ cw * ac, ot of the cable element to the 
lunar surface 

• o » Stefan Uoll/.m;r i constant (5.67 X I0~* 
W-scc/m* 1 , 

t M -hmai-suif .remittance 
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( c * cable-element cinillancc 
a cin a c^lc-dcmcnt inftarcd absorptancc 
V * cable-clement volume 
pc = volumetric heat capacity of the cable ele- 
ment 

T » * lunar-surface brightness temperature 

T c = 'able-clement temperature 

For « M = 1 and a flat lunar surface, equation (11-9) 
reduces to 

’« * *c>c) (" >») 

a ClH \ 0 ° / 

The term ^ which accounts for the thermal 

time constant of the cylinder, is retained only for 
eclipse calculations, be nuse the constant is on the 
older of minutes. 

The orientations of the thermocouples that are 
outside the borestems arc unknown. The lime at 
which a p’vcn thermocouple reaches the maximum 
temperature and the value of that maximum are 
strong functions of the o' rotation of the cable 
section in which the thermocouple is embedded. This 
effect, which is shown in figure 11-1!, is a result 
primarily of the variation in the incidence angle of 
the solar radiation; therefore, lunar-surface brigh tness 
temperatures deduced during t'te lunar day are 
subject to error that results from the orientation 
uncertainties. 

For the calculation of the lunar-surface brightness 
temperatures, the assumption was made that the view 
factors from the cable to sp.ioc and from the cable to 
the surface arc identical. The irregular horizon 
formed by the Apenninc Mountains increases the 
effective view factor to the lunar surface, and the 
view factor to space is reduced coi respond ingly. 
Some of the thermocouples may be close enough to 
the surface so that local topographic irregularities 
affect the horizon seen by the cable. For a 10-pcreciit. 
increase in the effective view area to the topographic 
surfaces, the calculated surface brightness tempera- 
tures would be reduced by 2.25 percent during the 
lunar night.' 

A more serious error in lunar-night calculations 
results from the uncertainty in the values of the cable 
absoiptancc auJ omittance, lor the tcmpeiatuies 
given in figure I! 12, the cntitiancc-to-absoi piano? 
ratio was assumed to be unity. A 20-p«.recnt increase 


in this ratio would result in an increase of 4,5 percent 
in the calculated value ot the surface brightness 
temperature. A value of 0.97 was chosen Ibr both the 
cable infrared cmiltaiiec and aosoiptuncc, where the 
assumption was nude that the sections of cable in 
which the thermocouples arc embedded are covered 
by a significant amount of lunar-surface material. 
Photogiaphs show that most of the cable areas arc 
coated with Umar material. 

The surfacc-brightncss-tempcrature history for the 
first iunar night, _s deduced from thermocouple 
temperatures*, is shown in figure 11-12. A very rapid 
cooling of the surface is indicated for the first SO hr 
after sunset; subsequently, the rate of cooling slows 
significantly. A cooling curse (ref. 11-15) based on 
astronomical observations at two different latitudes is 
shown for comparison. Two theoretical curves based 
on finite-difference calculations of thermal models of 
the lunar surface and subsurface ate aiso shown. 
Curve A, which is derived from a model with a linear 
conductivity increase starting at a depth of S cm 
(inset, fig. 11-12), duplicates the rate of cooling for 
tunes greater than 80 hr, whereas curve B, which is 
based on a tcmpcrature-dependcnt-conductivity 
model (ref. 11-12), duplicates the earlier part of the 
observed curve. For both models, the heat capacity 
was defined by the following -ouation (from ref. 
11-16) 

dT) = -0.034T 1/2 + 0.008T - 0.0002T 3 ' 2 

( 11 - 11 ) 

The assumed densities are 1.2 g/cm 3 for cune A and 
l.C g/cm 3 for curve B. The flattening of the observed 
curve could, in pari, be a result of a significant 
increase in density in the upper several centimeters. 

The temperature dependence of conductivity that 
is indicated by flic cooling curves is in agree ml with 
the large increase in mean temperature nesenbed 
previously. The substantial increase in r-mciuclivity 
and density with depth that is suggc.dij by »' 
liaitr ”••• of the cooling curves is in agieenirn 
earlie. ..inclusions that were based on astro- . , 
obseivations (e.g.. ref. 11-16). The mechanical prop- 
erties of the soil measured near the Apollo 15 A LSI . I* 
site revealed that the shear '.rcngih increased rapidly 
with depth in the upper 20 cm (see. 7). This increase 
in shear strength probably is related to a ncar-surfacc 
density gradient. 

Suifaec brightness temperatures during the umbral 
stage of the August 6 eclipse were dedu 1 from ;i . 
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thermocouple temperatures. For these calculations, 
the thermal time constant of the cable must be 
considered. The thermal mass of the cable per unit 
length in which the thermocouple is embedded was 
estimated by summing the properties of the 35 
conductors. Uncertainties result from the fact that 
noise in the thermocouple measurements makes an 
accurate determination of the slope difficult. A 
sample of thermocouple data for the later part of the 
umbra was reduced by determining lire rate of change 
of the cooling curve by graphical techniques. These 
very preliminary results indicate that the umbral 
temperatures reached at the lunar surface correspond 
well with the temperatures predicted by the theoreti- 
cal curve based on the relationship between conduc- 
tivity and temperature for lunar fines (ref. 11-12). 

In summary, the surface-temperature data for 
lunar night and the umbral part of ilic August 6 
eclipse support tiie conclusion that the upper few 
Centimeters of the lunar-surface ma;;rlal have 


thermal conductivity-versus-tcmpcrature relationship 
similar to that found for samples of lunar fines 
measured in the laboratory. The lunar nighttime 
observations reveal a substantial conductivity gradient 
with depth that probably results from increasing 
material density with depth. 


DISCUSSION OF HEAT-FLOW- 
EXPERIMENT RESULTS 

Local Topographic Effects 

The heat-flow determinations at die Hadley Rille 
sin* art susceptible to a number of disturbances. 
Corrections for some of these disturbances (such as 
the 'Jvumal perturbations caused b> visible lopo- 
giapliic feat ures) can readily be made with sufficient 
accuracy , but other distmLunccs may result from 
refr return associated with sloping ink. faces between 
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FIGURE 1 l-I2.-Surficr cooling as a function of time for 
the fir*t lunar night, as dcuiiced fjom llicrmocouplv 
meast'rcmc.ui Two theoretical curses arc given as refer- 
ence*. Olive A is the <• x-ct'J curve for a conductivity 
increase with depth, sia; ‘ng j» c depth of 8 cm; cuivc DU 
the theoiencal curve for a t • i./x.i'turc-dcpendc.it con 
ductivity. ’he data from refer -tit . 1 1 !6 are included for 
comparison. Conductivity as a function of depth for the 
curve A data is shown in the inset. 


materials of dili'eiing thermal conductivity. Suclt 
interfaces, if present at all. ate hidden in the lunai 
interior, 3nd onl> a qualitative discussion of the 
effects can be given. 

Probe I, from which data were used for the 
preliminary iletc. iination o'* heat How, is located in 
a small crater that has been aln tilled ai.J nearly 
obliterated by later boinbaidmepi. The topography 
associated with this old crater is so subdued that no 
correction for the lopuguphy is required. However, 
the possibility remains that the event 'hat produced 
fills crater also locally all ur'd both tin. thickness and 


the physical properties of the rcgohtlt. The irul't of 
this hypothesis cannot be verified because no observa- 
tional evidence is avail 'ble; however, after lire data 
from probe 2 have been .nul/zcd, a comparison of 
•he two heat-flow values may provide more informa- 
tion. 

• The two most conspicuous and important topo- 
graphic features near the hcat-flow-CNpcrMieut site 
arc Hadley Rdle and the Apcnninc ! : tont. The 
topographic effect of Uk riilc was calculated by 
fitting the riilc profile with a two-dimensional Lees- 
type valley (refs. 11-17 and 11-18). This procedure 
results in a correction of 4.5 percent, ai d further 
allowance tor the Elbow, where the riilc abruptly 
charges direction, reduces the correction to 3.5 
percent. 

A seed'd effect of the topography around the 
Hadley Rille site is that the surface in titc area is 
shaded during pari of the day; cc jucittly, rhe 
average surface temperature is low., man ilat por- 
tions of the lunar surface, litis local cold-spot cff"ct 
has been estimated quantitatively in two ways. The 
radiation balance for a poLt halfway down the rille 
wall (including factors for incoming solai radiation 
and radiation front tite opposite wall), at appropriaie 
tunes of lb . lunai day, wms calculated and used to 
derive the mean tempo. ature in the riilc. Alterna- 
tively, the temperature at the vertex of the riilc. 
which was as'umcd to have a symmetrical V-shape, 
was calculated from the sola, input alone; for this 
geometry, the vertex secs neither wall. The tempera- 
ture w , as assumed to var/ linearly from the vertex to 
the lop of the riilc. The tadiation-uala .c method 
resulted in a correction of 5 to 10 percent, and the 
vertex method vsah.J in a correction of 10 to 20 
percent. Th former va.ut is considered to be the 
more 'lial i ■ mainly because the profile of the rtlie 
does not particularly resemble a V-shape. A cc sec- 
tion of 10 percent is considered to be reasonably 
close to the tippet limit of th. effect that results front 
the cold riilc wall. 

The thermal effect of the Ap .mine Front was 
estimated in a preliminary way from the two-dimen- 
sional-slope method of Laclicnbrvch (ref 1 1 -19 ) Uy 
assuming two-dimensional symmcny, me curve; £>«. i 
H** acltc-rbruch indicate that *he ei -icuion. ■ .Inch is 
negative i*i this ease, is 10 pvt cm a i c.imurn. 
Am.nllv, tie heat How-cxperii ci.f site is in an 
emha, .nent in the Apcr.itiiiC l imit, which invalidates 
the •wo-duitersronai approximation and reduces the 
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correction. A correction of 5 percent would be more 
realistic. 

In sununaiy. the topographic effect of Hadley 
Kille is a ppr>\j. n.itely canceled by the effect of tl»c 
Apennine Front, which leaves only the cold-spot 
effect of the ntlc as a remaining correction resulting 
from the visible topography. This latter correction is 
the largest correction, in any ease, and seems most 
likely to approximate 10 percent of the measured 
heat How. 

Implications of the Results 

In this section, the view was adopted that the heat 
tlow obscived at the Hadley Rille site is representa- 
tive of the moonwide value, in spite of the possibili- 
ties of local and regional disturbances and large-scale 
variations in the 'unar heat flow. Thus, the measure- 
ment is considered at face value, with full realization 
that future measurements may produce major 
changes in the conclusions. Because of the prelimi- 
nary nature of the result:, simple models were used 
for this report. Only ihe linear equation of heat 
conduction was considered, for which the thermal 
diffusivity was constant. 

The value that is considered to be an upper bmit 
to the heat flow resulting from the initial heat is 
calculated by assuming that, at the end of a convec- 
tive stage in the early history of the Moon, tempera- 
tures throughout the Moon lay along the solidus for 
lunar basalt (ref. 11-20). After 3 X 10 9 yr, the heat 
flow resulting from these very high initial tempera- 
tures is in the range of 0.2 X 10 -6 to 0.4 X 10"* 
W/cm*. If, as suggested bv the greater ages of all 
known lunar rocks, partial melting throughout the 
Moon lock place earlier, still lower values of heat 
flow would be associated with the initial heat. It is 
likely tha: the assumed initial temperatures arc too 
high; however, by revising the assumed temperatures 
downward, the flux flo:« the initial beat is further 
reduced. The initial he it contributes little to lire 


present luiiai heat flow, and the eontiihutiim can be 
neglected for picsent purposes. 

The major portion of the heal flow from the Moon 
probably results from radioactive heat generation in 
the interior. It is possible to construct an infinite 
number of models based on nniumiform distributions 
of radioactivity; however, in this report, the discus- 
sion was con lined to consideration ol the Moon as a 
sphere with uniform and constant internal heat 
generation. The ratio of the surface heal flow q 
(expressed in I0 - * W/cm : ) and the heat production 
Q (expressed in I0~' 3 W/cm 3 ) is shown for several 
times in the following list. 


Time, yr 


q/Q. cm 


3 X 10* 

4 

5 

Infinite 


3.22 X 10 7 
3.53 
3.89 
5.78 


With a lunar heat flow of 3J X I0“* W/cm 1 , the 
value of Q must be in the range of 0.57 X !0~ |J to 
IDX 10*“\V/cnt 3 .This number :s far lower than the heat 
production of lunar basalt, which has a value of 
approximately 3.5 X KT 1 3 W/cm 3 . However, the 
basaltic rocks ate presumably differentiates that are 
far more radioactive than the parent material. On the 
other hand, ordinary choudritc and type- l carbona- 
ceous chondrites generate heat at rates of approxi- 
mately 0.17 X 10"' 3 and 9.22 X 10“' 3 W/ent 3 , 
respectively. The icspcciive average rates of heat 
generation over the last 4.5 X I0 9 yr arc 0.61 X 
10-' 3 and 0.67 X lO"' 3 W/cm 3 . Even these last 
figures are barely sufficient to provide the necessary 
flux. The conclusion is that, if the observed lunar 
heat flow originates from radioactivity, then the 
Moon must be more radioactive than the classes of 
meteorites that have formed the basis of Earth and 
Moon models in the past. 
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APPENDIX 

Equilibration of an Infinitely Long Cylinder 
in a Homogeneous Medium by Conduction 


lo this appenJix. the theory of the cooling of a 
cylinder by conduction is reviewed and applied to the 
heat-flow probe. Cylindricnlly slraped probes are 
commonly used in making geothermal measurements, 
and. as a consequence, the theory of heat flow in 
cylindrical coordinates has been thoroughly investi- 
gated (refs. 11-11, 11-14, and 1 1-21). The investiga- 
tion of tire effect of a finite contact resistance 
between a cylinder and the surrounding medium (ref. 
1 1-14) is particularly applicable to problems of tire 
cooling of the lunar probes. In this appendix, the 
solution (ref. 1 1-14) has been extended to a some- 
what more complex model that includes a solid 
cylinder inside a thin-walled concentric cylinder 
which, in turn, is surrounded by an infinite medium 
with a conductivity k. The cylinders are assumed to 
be perfect conductors: that is. each cylinder is 
isothermal, which is very nearly true in the case of 
tire heat-flow probe. Contact resistances exist at the 
two cylindrical boundaries. In this model, the inner 
cylinder is an idealization of the heat-flow probe, 
which is radiatively coupled to the borestem (the 
concentric cylinder). The br'cstcm. in turn, loses 
heat by conduction to lire surrounding regolith, of 
conductivity A. through a thin zone of lunar material 
that has been disturbed by drilling and. hence, has a 
different conductivity. The initial temperature of the 
suriounding infinite medium is zero; and. initially, 
tire two inner cylinders have a temperature r’ 0 . 

By defining the dimensionless parameter r as ni/b 1 
(where k is lire thermal dil’lusivity ol the surrounding 
infinite medium, t is time in seconds, and b is the 
radius of the outermost cylinder), the ratio of the 
inncr-cylindcr temperature i-(r) to the inner-cylinder 
initial lempeiaturc i’ 0 is given by 


a - 2nb 1 pc/(S l + S 2 ) 
h = k/bH 2 

A = SjfK/b'iSi +S,) 

B = yic/b 2 
y = S,/2ira//, 

a = radius of the inner cylinder (heat-flow probe) 
in centimeters 

b = outer radius of the concentric outer cylinder 
(borestem) in centimeters 

S, £ 2 = thermal heat capacities per unit length of tlx 
inner and outer cylinders, respectively, 
W-scc/cm-K 

H,Jh = contact conductances at the inner- and outer- 
cylinder boundaries, respectively, \V/cm : -K. 
p = density of the surrounding material, g/cr.i 5 
c = heat capacity per unit mass of the surround- 
ing material, W-sec/g-K. 

k= thermal conductivity of the surrounding ma- 
terial, W/cnr-K 

k - thermal difTusivity (A/pc) of the surrounding 
material, cm 1 /sec 

Also,y 0 («).-A|(«). ^oOO- 2 nd y t (k) are the zero- and 
first order Bessel functions of the first and second 
kinds. 

By following the method described in reference 
11-11, tire expression on the right-hand side ol" 
equation (11-12) can be redesignated j(A,Iia t)\ 
therefore, equation (11-12) can be written asr(T}/i’ 0 
= A more complete derivation of evira- 

tion (1 1-12) and the resulting tabulated valuer vi!l be 
published at a later date. Two typical pints of 
JJA.h, q.t) are shown in figure I l-S. The value* ol A 
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and // j dial were used for computing die theoretical 
curves arc given in table lit!. Other parameters area 
= 3.33. //, = 3.25 X UT* NY/cin 1 • K, y = 445 sec, 
and a and b ate 0.95 and l .259 cm, respectively. 

By using asymptotic values of the Bessel functions 
for large values of r (after a method outlined in ref. 
11-21). a solution that is valid for long times can be 
derived as follows. 


•M 

* 



CU-13) 


where Jn (5) = 0.5772 (Euler's constant). The 
long-term solution exhibits some interesting features 
of the function/d A<*,r). At long times (t > 20). the 
function becomes nearly equal to (Zor) -1 , or 
(4fcrf/S, + S;) _l , which is independent of pc. A 
log-log plot cf/dAa,*) is shown in figure II -S; for 
long times, the function describes a straight line with 
a slope of very' nearly -1. This property of the 
function enables die temperature histories to be 
extrapolated to equilibrium values by the procedure 
described in the text of this section and a value lor k 
to be determined independently of pc from long-time 
poitions of the equilibration curves. In addition, as 
shown by equation (J M3), the contact conductances 
and // 2 , which are contained in the constants y 
and ft, become less important with time (because they 
arc multiplied by f -1 ). 
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HEAT HUW EXPERIMENT 


Introduction and Background 

The objective of the Heat Flow Experiment (HFE) is to make a direct 
measurement of the vertical component of heat flow from the lunar interior 
through the surface.. A second important objective to 

determine the thermal properties of the upper 3 m of the lunar regolith. 

The age of the Moon is placed at 4.6 AE. For a planetary bocy as small 
as the Moon, much of its initial heat energy has been lest to space since 
formation. Even if the Moon were initially at molten temperatures, the present 
flux at the surface would be small. The major contribution to the surface 
heat flow comes from heat generated by the disintegration of long-lived radio 
isotopes of 235 u, 238^ 40 K an( j 232-j-^ Thus p reser) t surface heat flux 
reflects the abundance of these isotopes to a depth e about 300 km, or 43 " 
of the volume of the Moon. It is now certain that extensive differentiation 
occurred during the early history of the Moon which would concentrate these 
isotopes in its outer shells. In this case the present surface heat flow would 
very nearly indicate the total abundance of these isotopes. 

Over 5000 heat flow determinations have been made on Earth and they 

show the average global flux is 6.3 x 10“ 6 W/cm^ (ref. l). (Throughout this 

report we will use W-sec as the unit of heat energy. Other commonly used units 

are the calorie which equals 4.18 W-sec, and the erg which equals 10"^ W-sec. 

The average Earth heat flow is 1.5 x 10" 6 cal/cm 2 -sec and 63 ergs/cm ? -sec. ) 

rate of 

Urey (ref. ) pointed out that the total/ heat flow from the Earth is essentially 
, . , rate of 

equal to the total/ heat production in the Earth if it were constituted of materials 
with chondritic abundants of U, K and Th. However, Gast (ref. 3) showed that 
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Earth rocks were strongly depleted in potassium relative to solar and chondritic 
abundances, which led Wasserburg et al. (ref. 4 ) to propose that to explain 
the Earth's heat flow, a higher abundance of uranium, nearly 3 times that of 
chondrites, is needed. He estimated a uranium abundance of about 30 ppb for 
the Earth. 

Large variations of heat flow are observed over the surface of the Earth. 
These variations are principally related to the present tectonism of the Earth's 
lithosphere (ref. 5). The largest variations are observed at extensional and 
compressional boundaries of vast litho phene plates that are moving relative 
to each other. Seismic observations (ref. 6 ) and the preservation of ancient 
surface features on the Moon demonstrate that no comparable tectonic movements 
have occurred on the Moon for the past 3 AE. The Moon is tectonically dead 
compared to the Earth and therefore any variations in surface heat flow over the 
surface should reflect either deep-seated changes in abundances of radio isotopes or 
convective patterns in the Moon. The static nature of the Moon's outer crust 
makes it an attractive planet for heat flow measurements, since determinations 
at a single location might be quite representative of a very large region of 
the Moon if local surficial effects such as refraction by conductive inhomo- 
geneities and topography are properly accounted for. 

Numerous attempts have been made to determine the surface heat flow 
from the Moon by detecting thermal radiation from the Moon in the microwave 
band. Due to tie partial transparency of lunar surface material, energy with 
wavelengths greater than 1 mm received at earth-based antennas originates 
in the subsurface and contains information on subsurface temperatures. By 
making estimatesof thermal and electrical properties the heat flow can be 
determined from the change of brightness temperature with wavelength. The 
earliest heat-flow determination was that of Baldwin (ref. 7) who estimated 
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an upper limit of 1 x 10 W/cm . Russian investigators estimated the heat flow 
to be very nearly equal to that of the Earth based on a very careful set of 
radio telescope observations in the wave length band 3-50 cm (ref. 8 ). These 
same measurements were later revised to give heat flows in the range 3 to 4 
x 10 W/cm using different thermal properties and a layered model (ref. 9 ). 

On Apollo 15 the first direct measurement of heat flow through the 
lunar surface was made at Hadley Rille (26°06'N and 3°39'E). Now a second 
measurement has been made on Apollo 17 at the Taurus Littrow site (20°09'N and 
30°45'E). At Taurus Littrow two probes to determine heat flow were emplanted about 
11 m apart. Analysis of data taken during the first 45 days after emplace- 
ment indicates that the heat flow from one probe location is 2.8 x 10“ c W/cm 2 
(0.67 ucal/cm 2 -sec) and 2.5 x 10 -6 W/cm 2 (0.60 ycal/cm 2 -sec) at the second 
probe location. For comparison the value measured at Hadley Rille is 3.1 x 10 
W/uv 2 . The Hadley Rille measurement and the Taurus Littrow probe #1 measurement 
have an estimated error of +20%. The probe #2 measurement at Taurus Littrow has 
a slightly larger uncertainty because the heat flow appears to be locally disturbed 

Experiment Concept and Design 

The concept upon which the Heat Flow Experiment is based is to measure 
directly the vertical flow of heat through the regolith. The measurement 
should be made far enough below the surface so that the time varying heat flow 
resulting from the very large diurnal variations of surface temperature is small 
compared to the flow from the interior. At Taurus Littrow this depth is about 
100 cm. Below this depth the increase in temperature with depth principally 
results from the hotter lunar interior. The outward flow of heat, F z , is 
directly proportional to the rate of temperature increase with depth, dT/dz: 



(Eq.l) 


The constant is the thermal conductivity. The negative sign simply in- 
dicates that heat flows in a direction opposite to the increase of tempera- 
ture. 


The experiment is designed to accurately measure the vertical tempera- 
ture gradient in the lunar soil to a depth of 2.3 m. Surface temperature 
measurements are also made which can be used to deduce the thermal properties 
of the upper 10 to 15 cm of the regolith. In situ measurements of thermal 
conductivity of the regolith at depths where the gradient is measured are 
also carried out. Two measurementsof heat flow at locations separated by 
about 10 m are made to detect possible lateral variations. 


Instrument Design: The essential parts of the heat flow instrument are two 

identical temperature sensing probes. {A drawing of a probe is shown in 
figure 1.) Each probe consists of two sections about 50 cm long. In each 
probe section there are two platinum resistance bridges. Each bridge con- 
sists of four 500-ohm filamentary platinum elements interconnected by evan- 
ohm wire. (A diagram of the bridge is also shown in figure 1.) Opposing 
arms of a bridge are wound together in a single sensor housing, and two sensor 
housings, comprising a complete bridge, are mounted at opposite euds of a 
probe section. 

Voltage measurements on a bridge can be interpreted in terms of 
average bridge temperature, T a , and temperature difference, aT, between sensors, 
by means of accurate calibrations. The temperature at each sensor can be 


determined simply by T a + 1/2 AT. The accuracies of the heat flow experiment 
temperature measurements are given in Table 1. 

The cable thermocouples consist of a string of 4 chromel-constantan 
junctions embedded in each probe cable. The lowermost junction is positioned 
inside the gradient sensor housing at the top of the probe. The reference 
junction for each cable is inside the electronics housing and is mounted in an 
isothermal block with a platinum resistance thermometer (the reference 
thermometer). The thermocouple circuit is diagramatically shown in figure 1. 

The accuracy of the thermocouple measurements has special signifi- 
cance for interpreting the subsurface temperature profiles at the Apollo 17 
site. Therefore certain features of the thermocouple measurement, which 
affect the accuracy, should be described. Firstly only the chrome! - con- 
stantan junction TCI inside the topmost gradient is coupled with the ref- 
erence junction during a measurement sequence. The other junctions TC2, 

TC3 and TC4 are coupled with TCI so that, in effect, TCI becomes the reference 
junction and the top gradient sensor becomes the reference thermometer. Since 
the temperatures of TC3 and TC4 are much closer to that of TCI than to that of 
the thermometer in the electronics box the calibration errors are reduced by 
this technique. Secondly, the placement of junction TCI inside the gradient 
sensor permits a direct comparison between the two sensors on the Moon. Thus 
an in situ calibration of the thermocouple circuit against the much more ac- 
curate gradient sensor is carried out. 

Conductivity experiments are made using 1000-ohm heaters that surround 
each of the gradient sensor housings. The experiments can be operated in 
either of two modes by energizing the heaters at 0.002 W or 0.5 W depending on 
the conductivity of the surrounding regolith. Because of the low conductivity 
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of lunar material at Taurus Littrow, only experiments using the lower power 
have been made on the Apollo 17 experiment. After initiation, power is left 
on for 36 or more hours and the temperature rise of the sen-^r beneath the 
heater is precisely measured by the appropriate bridge. /' more - 'tailed 
description of the operation of the heat flow experiment can be found in ref. 11. 

Deployment of the Experiment at Taurus Littrow 

The heat flow experiment is deployed as part of the Apollo Lunar Surface 
Experiments Package (ALSEP) and is located about 200 m west of the LM landing 
site. The electronics housing is located about 9 m north of the ALSEP central 
station and the two probes are emplanted approximately 5.5 m on ei ther side of 
the electronics housing on approximately an E - W line. The 
ALSEP site is on a local topographic high, perhaps an inter-crater ridge 
between two wide but shallow depressions north and southeast of the site. The 
two probes are emplanted near the rim of the shallow northern depression. 

The borehole for probe #1 was drilled east of the electronics housing to 
a depth of 2.54 m, and the probe was inserted to the maximum possible depth, 

2.36 m. The rim of a small, relatively recent crater about 2 m in diameter 
and 0.3 m deep lies about 1 m northwest of this borehole. One of the thermo- 
couples, TC12, is in a portion of the cable exposed above the surface. 

This portion of the cable, which is covered by a black sleeve, is 
oriented nearly N - S arid is tipped up toward the south about 15° relative 
to the surface. 

Probe §2 is buried in a borehole drilled west of the electronics 
housing to a depth of 2.55 m. The rim of a 3 in crater is about 2 m north- 
west of this probe site. Thermocouple TC22 is exposed above the surface, 
aligned nearly N-3 and is almost level relative to the surface. 
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Summary of the Theory Relevant to the Lunar He al Flow Measurement 

In this section we will review the important aspects of heat transfer 
that are essential to interpretation of our temperature and conductivity 
measurements. For the reader primarily interested in the experiment's results 
this section may be skipped and used for reference only. The discussion will 
be summary in nature and for more detailed description, the reader is referred 
to references 10 end 12. 

The Initial Cool Down of the Probe; Estimates of T he rmal Conductivity of 
the Regolith and Equilibrium Temperatures 

When the bores tern is drilled into the soil and the probe inserted, 
they are at temperatures considerably higher than their subsurface sur- 
roundings. These higher temperatures result from the temperature of these 
components above the surface before emplacement and heat generated during 
drilling. Because of the complexity of the drilling process, it i c y 
difficult to estimate the amount of heat dissipated and its distn *on 

along the borestem. We know from measurements made within minutes after 
the probes were inserted that the probe is about 40 to 50K warmer than 
the surrounding undisturbed lunar soil. Because of the very low thermal 
conductivity of the lunar soil, this heat is dissipated very slowly. Even 
after 45 days some of the thermometers indicate continued cooling. 

In general, the rate of equilibration of each probe sensor to the 
surrounding lunar medium depends on a number of different parameters, including 
thermal and geometrical properties of the probe - borestem system, thermal 
properties of the adjacent lunar medium, the thermal contact between the borestem 
and adjacent zone of disrupted lunar maters! and the total heat energy excess 
of the system upon completion of drilling. For long times (t > 20 hr) after the 
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probes are inserted, however, the probe sensor temperatures above equilibrium 
depend essentially on just two quantities, the total . initial energy excess 
per unit length of the system, ae, and the thermal conductivity of the adjacent 
regolith, : . 


T p (t) . T, ' fiE/^t (Eq. 2) 

where T (t) = a probe sensor temperature at time t since emplacement 
and T ro - the equilibrium sensor temperature. 


The form of equation 2 has been verified by finite difference models 
as well as an analytical solution appropriate to a simplified distribution of 
AE between the probe and bores tern. The analytical solution, which takes the 
form of equation 2 for long times, is described in detail in the appendix of 
section 11 of the Apollo 15 Preliminary science Report (ref. 10). 


Estimates of the regolith conductivity, K m , can be made from the long 

, . , . , ^ . .equation 2 , „ 

time cool-down data using / once a value for AE is assumed. Because of the 

uncertainties of the contribution of the drilling process to the total excess 

heat energy, two different assumptions have been made in regard to the aE term: 

first, that drill-heating effects are negligible and that all excess heat energy 

initially resides in the probe and borestem. The probe and borestem are 

assumed to be at the same initial temperature, 1 Q , following drilling, as measured 

by the probe sensors. For this assumption the initial energy excess can then 

be expressed simply as: 


AE = (S p + S b ) (T 0 - TJ 


(Eq. 3) 
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where S p and S b are the heat capacities per unit length of the probe and bore- 
stem respectively. The resultant problem has an analytical solution which is 
used to determine regolith conductivity estimates from the cool-down data. The 
above assumptions constitute the "minimum initial energy" case. 

As a maximum initial energy estimate, it w«s assumed that a 2.2 mm 
thickness of disrupted lunar material had received enough heat during 
the drilling process to attain an initial temperature equal to that of 
the probe and borestem. The 2.2 mm contact zone thickness was chosen 
to correlate with the difference in drill -bit radius and borestem radius. 

The subsequent cool-down problem was then solved, assuming radial heat 
dissipation, using a finite difference model. The curve matching proce- 
dure was carried out tor times > 20 hr when all parameter effects, with 
the exception of k m ,are negligible; i.e., when equation 2 is valid. 

Because of the form of equation 2, the probe equilibrium tempera- 
tures can be estimated quite accurate. y independent of initial conditions 
and k^. By considering sensor temperatures',.^ and Tp 2 for times and 
t 2 >>50 h» , we may eliminate the unknov.;. factor AE/47tk m from equation 2 
to obtain: 



(Eq. 4 ) 


Equilibrium temperatures differences, measured by the platinum resistance 
bridges, can be calculated in a similar manner. 
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Heat-Activated Conductivity Measurements 

Since all eight conductivity measurements were made with a heater 
power of .002 W, this discussion will be confined to this mode of experi- 
ment. After the heater is energized, the temperature rise at the gradient 
sensor it encloses depends in a complex way on the thermal properties of 
nearby probe components, borestem and lunar material, as well as the 

contact zones between these elements. A detailed finite difference model 
conductivity 

of the /experiment at each heater location is used to interpret the tempera- 
ture history of the gradient sensor in terms of the conductivity of 
material external to the borestem. Numerical model computations, labora- 
tory experiments and lunar experiments indicate that for times > 20 hrs 
the temperature rise, Af(t), is well defined by the simple relation: 

aT( t) = c^ln(t) + C£ ( tq . 5 ) 

where c^ and C 2 are constants. The form of equation 5 is the same as 
that for a heated infinite cylinder in an infinite homogeneous medium 
at long times; i.e. > 20 hr for a cylindrical source with a radius and 
heat capacity per unit length of the HFE probe - borestem system (ref. 13). 

For an infinite cylindrical source 

C 1 = Q/Hr, (Eq.6) 

where Q is the power per unit length in W/cn. Thus c 1 depends s ,ely on 
the heater power per unit length and conductivity. The constant c^ can be 
easily determined since it is the slope of the temperature rise curve when 
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plotted versus In t; therefore cylindrical sources are often used as a 
practical technique for measuring conductivity. 

Conductivity is determined from lunar experiments by comparing ob- 
served slopes on a logarithmic time scale with values of c^ calculated 
with the finite difference models. Parametric studies, in which certain 
thermal properties are varied singly in the numerical model, show that for 
times > 20 hr is very nearly insensitive to: changes in pc of the 
surrounding medium, changes in borestem conductance and changes in the ther- 
mal links between the probe and borestem and the borestem and lunar medium. 
However, c^ is sensitive to changes in conductance in the probe body which 
can alter the flow of heat from the heater axially along the probe. As- 
sumptions of thermal properties in the numerical models that influence 
axial heat transfer along the probe are probably the largest source of error 
in the conductivity determinations. 

The similarity in performance of the lunar conductivity experiment 
and an infinite cylindrical source is principally due to the relatively 
efficient flow of heat axially along the borestem. Even though the probe 
heater is very short, 1.9 cm, it heats a section of borestem that is long 
compared to the borestem diameter. For times > 20 hr the isotherms in the 
surrounding medium are roughly cylindrical in the vicinity of the heater 
as is shown in figure 2. The numerical computations also show that the 
experiment is most sensitive to lunar material within about 5 cm of the 
borestem wall. 

The effective conductance of the contact zone has a very pronounced 
effect on the magnitude of the sensor temperature rise at any given time. 
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Since can be determined independently the conductance of the contact 
zone is the principal remaining unknown parameter, and can be determined 
by matching observed and theoretical temperature curves for times >0.5 hr. 

Variations in Surface Temperature and its Effect at Depth 

Lunar surface temperatures vary nearly 300K from just before lunar 
dawn to lun... noon. This variation induces subsurface variations that 
propagate downward as thermal waves. For a homogeneous medium of diffusiv- 
ity a with a sinusoidal variation, A 0 cos «t at the surface, the temperature 
at a given depth, z, is given by 


T(t, z) = A 0 e“ 6z cos (ut - Sz) 


(Eq. 7 ) 


where P. Q is the amplitude cf the surface variation (degrees), u> the angular 

frequency (rad/sec) (2.5 x 10~ 6 for the diurnal variation and 2 x 10“ 7 for 

the annual variation), and <5 = /^(cnr 1 ). (Eq. 8) 

V 2a 


Equation 7 indicates that the variation decreases in amplitude by a 
factor e" 1 and is delayed in phase one radian for every 5 -1 centimeters of 
depth. 

Propagation of surface temperature variations into the lunar regolith 
is more complex for a number of reasons. Firstly, the surface variation 
is not a simple sinusoid but contains significant higher harmonics. Sec- 
ondly, thermal properties vary significantly with depth; and thirdly, 
radiative transfer, which depends on T 3 , plays an important role in the 
upper few centimeters of the lunar soil. It is necessary to resort to 
numerical calculations which include these complications to determine the 
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expected temperature variations in the subsurface. In figure 3 the peak- 
to-peak attenuation and phase lag of the diurnal variation is shown as a 
function of depth for the conductivity profile at the Apollo 17 heat flow 
site. The upper part of the conductivity profile is derived from surface 
temperature measurements which are described in a later section of 
this report. Notice that for depths greater than a few centimeters the 
amplitude decreases in a simple exponential fashion as evidenced by the 
nearly straight line on a semilog scale. Similarly the phase lag shows a 
nearly linear increase with depth below a few centimeters. Thus the 
simple relation of equation 7 would nearly apply below these depths. 

The temperature at lunar noon varies throughout the year due to the 
varying distance of the Earth - Moon system from the Sun. The noon tem- 
perature increases about 6K from aphelion to perihelion. The mean temper- 
ature, i.e., surface temperature averaged over a lunation, varies about 
3K throughout a year. Although the amplitude of the annual cycle is a 
hundredth of the diurnal variation, the decay constant 6 is /T? times 
smaller; consequently, annual variations penetrate deeper and induce 
significant heat flows to depths of a few meters and must be considered 
ir. interpretation of the experiment's results. The attenuation of ampli- 
tude and increase in phase lag for the annual wave is shown in figure 3 as 
a function of depth. Annual wave effects shown in figure 3 are based on 
the conductivity profile at Hadley Rille. 

From figure 3 it can be seen that temperature fluctuations due to the 
diurnal cycle become virtually undetectable at depths > 100 cm and would 
have had little effect on heat flow below this depth before the probe and 
borestem were emplaced. Once the borestem is emplanted in the regolith. 






its higher conductivity and radiative transfer inside the borestem will en- 
hance the propagation of thermal waves downward. However, thermometers at 
130 cm below the surface do not detect any temperature variation during a 
lunation cycle. 

Corrections for the Shunting Effects of the Borestem and Probe 

The axial conductance of the epoxy borestem is considerably higher 
than a vertical column of lunar soil of equal cross section. This fact, 
combined with the finite length of the borestem, results in some shunting 
of the steady-state heat flow through the borestem to the surface. Certain 
short sections of the borestem such as the bit anc joints are made of 
titanium or steel and sizable disturbances occur near these parts. A 
second related effect results from the fact that the probes are radiatively 
coupled to the borestem walls and have a small axial conductance. This 
causes the probe bridges io read slightly smaller temperature differences 
than points on the borestem immediately adjacent to sensors. 

These effects can be estimated by simplified analytical models and by 
laboratory experiments; both methods were used in our earlier analysis of 
Apollo 15 results. However, for the Apollo 17 analysis and the refined 
Apollo 15 results presented in this report, we have resorted to a numerical 
model of the probe in a medium in which heat is flowing parallel to the 
probe axis. The numerical model besides being more detailed also allows 
examination of certain combined effects that are hard to estimate with 
analytical models. 

The numerical model computations show- that the borestem and probe dis- 
turbances to the steady- state heat flow are small. In the extreme case the 
temperature difference across a probe section is 7 % lewer than the tempera- 


ture difference across the same vertical distance interval far from the 
borestem. The numerical model has been used to apply corrections to all 
probe observations. 
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RESULTS 


Subsurface T emperatures 

The heat flow experiment was turned on while probe #2 was being in- 
serted into the borestem and temperatures were recorded only minutes after 
drilling was complete. These temperatures ranged from 295K to 301 K. The 

early cooling histories of probe #1 indicate similar initial temperatures. 
Subsequent to emplacement the probes cool toward the undisturbed rego- 
lith temperatures. Figure 4 shows the temperature histories of all sensors 
deeper than 65 cm for the first 45 days. After 45 days some sensors are 
continuing to cool: however, the expected future decrease is probably less 
than the error of absolute temperature measurement. 


The equilibrium temperature differences and absolute temperatures of 
each sensor are listed in Table 2. The correction for the steady-state 
disturbance of the heat flow by the borestem and probe system is applied 
to temperature data listed under the headings "corrected temperature diff- 
erence" and "corrected temperature." The appropriate corrections of the 
temperature difference due to the annual thermal wave during January 1973 are 
listed in the far right-hand column, but they have not been applied be- 
cause they are based on the conductivity profile at Hadley Rille. Note 
that the largest correction is about -42. 


Temperature measurements of the thermocouples TCI 4 and TC24 are also 
shown in figure 4. Some randomly sampled representative points are 
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shown and the solid curve is fitted by eye to show the trend. The 
standard deviation of the points around the smoothed curve is about 
0.3°C during the day and about half that value during the night. The 
values shown are calculated by subtracting the temperature at TCI 
from that at TC4 and adding the result to the temperature at the top 
gradient sensor. 

Comparison of the temperatures measured by TC11 and TC21 with 
those at the top gradient sensors show relatively 1 a to? errors in abso- 
lute temperature measurement. The errors are shown in Table 3 below. 

The source of these errors has been traced to copper - kovar (Cu/Ko) 
junctions in each thermocouple circuit chat are mounted on circuit 
boards in the electronics housing. The errors are proportional to 
temperature differences between the Cu/Ko junctions, c o that the errors 
in 7C11 and 7C21 are direct measures or this temperature difference and 
can be used to estimate errors at all junctions. The larger errors dur- 
ing the night result from larger temperature gradie^cs across the Cu/Ko 
junctions at night, and the larger errors at probe #2 junctions are 
caused by a greater distance between the Cu/Ko juntions in the probe #2 
circuit. A preliminary analysis of EMF's produced by Cu/Ko junctions 
was used to calculate the corrections that should be applied to the 
data in figure 4. These corrections virtually erase the apparent varia- 
tion between night and day. The corrections have been applied to the 
data compiled in Table 2. The uncertainty of determining these correc- 
tions is estimated as +0.4°C. Studies of the accuracies of the thermo- 
couple measurements are continuing. 

The values given in Table 2 represent the mean of the last 30 days 
of observation. The amplitude of the diurnal variation at this depth 
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cannot be determined with the present accuracy of the data, +0.5 K. 

In figure 5 the equilibrium temperatures are plotted as a function of 
depth. Temperatures along the body of probe #1 show a steady decrease in 
gradient with depth. The gradient decreases from .016 in the depth range 
130 to 177 cm to .012 K/cm in the range 185 to 233 cm. This decrease is 
principally due to a general increase in conductivity of the regolith over 
the interval of measurement. The thermocouple temperature indicates a 
gradient of .013 K/cm from 66 to 130 cm; however the accuracy of this 
measurement is poor. 

At probe #2 the probe thermometers (131 to 234 cm) indicate a 
rather uniform gradient of 0.0078 ./cm, whereas the gradient between 67 
and 131 cm is 0.021 K/cm, a change by a factor of about 3. This large 
variation of gradient can only be partially accounted for by the variation 
of conductivity of the regolith immediately surrounding the borestem. 

Apollo 15 Subsurface Temperatures 

Subsurface temperatures measured at Hadley Rille below the depth 
disturbed by diurnal variations were reported in the "Apollol5 Preliminary 
Science Report" without correction for the annual wave. In addition, 
corrections for the borestem and probe disturbance were derived from 
much simpler models than those discussed in this report. Temperatures 
measured at , onger times after probe insertion are now available and they 
allow a more accurate determination of equilibrium temperatures. The tem- 
peratures and temperature differences at four sensors on probe #1 at Hadley 



18 


Rillc arc presented in Table 4. More accurate corrections for the bore- 
stem disturbance and a correction for the annual wave effect have been 
applied. These measurements will be the basis for a slightly revised 
heat flow value. 

Conductivity Estimates from Cool- Down Analysis 

Because of the uncertainties in the total heat energy associated 
with drilling, .wo cases assuming different initial conditions have been 
examined. These cases have been described in the preceding section on 
theory. Results derived assuming initial borestem and contact zone 
temperatures equal to the initial probe sensor temperature are listed in 
Table 5 under the heading "With Drill Heating Effects." Conductivity 
estimates derived assuming that only the borestem and probe were initially 
at elevated temperatures are lifted under the heading "Without Drill 
Heating Effects." The two cases are considered to be bracketing assump- 
tions of the actual initial conditions. Cool-down conductivity esti- 
mates were made for each of the eight sensors along each probe. Addi- 
tionally, cool-down analyses were performed assuming drill heating 
effects for the thermocouples located 65 cm above each probe. 

The large noise on the thermocouple data limits the accuracy of con- 
ductivities deduced from their cooling history. Deductions of the conduc- 
tivity at depths from 3 to 15 cm below the surface, which will be discussed 

-4 

later in this report, give values of about 1.2 x 10 W/cm-k. Based on 
these results at shallow depths we estimate the conductivity lies in the 
range 1.0 to 1.6 x 10 -4 W/cm-k at 66 cm. . This range is indicated in 
figure 5. 


It can be seen by comparison with the more accurate heater-activated 
conductivity determinations discussed in the next section that drill 
heating effects cannot be neglected if reliable conductivity information 
is to be extracted from the cool-down data. When substantial drill 
heating effects are included in the cool-down analyses, conductivity 
determinations as well as variations with depth agree well with the 
heater-activated conductivity experiment results. The cool-down con- 
ductivity estimates are particularly valuable in interpolating between 
the more accurate heater-activated conductivity determinations. 

Heater-Activated Conductivity Experiments 

Conductivity experiment at each of the eight heater locations have 
been carried out. Figure 6 shows the sensor temperature rise history to- 
gether with theoretical curves for one such experiment. The conductivities 
and contact conductances H£ are given in Table 6. These results are 
shown in figure 5 together with the cool-down estimates with drill heating 
effects. It is clear that the conductivity does not vary in any simple 
way with depth. 

There is a rough correlation between drill penetration rate during bore- 
stem drilling and the measured conductivity. The more resistant layers where 
the drill penetrated slowly correspond to depths where higher conductivity is 
observed. The more resistant layers likely correspond to more compacted rego- 
lith materials or possibly to a higher concentration of centimeter size rock 
fragments. Either of these phenomenon can increase the bulk conductivity. 
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The relatively high cond; tivity measured at 130 cm on probe "1 lies within 
a zone from 80 to 130 r;n where penetration was slow. Directly below this 
layer, drilling rates were relatively high and the conductivity values are 
correspondingly lower. These correlations are used to interpolate values 
between discrete measurements. In figure 5 the solid line which passes 
through the heater-activated conductivity values represents a layered model of 
conductivity in the regolith based in part on penetration rates and in part on 
cool-down estimates. At probe #2, some of the drilling operation was not 
visually monitored so that correlations with conductivity cannot be made during 
the unmonitored period which includes about half the depth range where the 
probes are emplaced. 

One interesting feature of these conductivity results is a rather large 
difference between the conductivity profiles at probe it 1 and probe # 2. It is 
possible that layers as defined by conductivity have some dip relative to the 
surface. For example, the high conductivity layer at 100 cm at probe #1 could 

correspond to the high conductivity layer between 170 and 230 cm at 
probe H2. 

The contact conductance, arises from low conductivity material 
lying in a disturbed zone just outside the borestem. We estimate this 
zone to be 2.2 mm thick. The conductivity, K c , of material in the con- 
tact zone is given by 


K c = • tb + 1/2 Ar] In ( J ?'* Ar j 


(Eq.9) 
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See figure 2 for definition of parameters in this equation. As an example 
for H = 1.4 x 10"^ W/cm^-K, K = 3.0 x 10“ 5 W/cm-K. This value of con- 

Cm V* 

ductivity is about a factor of 6 less than the surrounding regolith. 


Apollo 15 Results 

Six conductivity experiments using a heater power of 0.002 W were 
run on Apollo 15 probes. The analyses of three of these measurements 
were not described in the "Apollo 15 Prelimin Science Report" because 

it was very difficult to separate changes due ^oter turn-oi. from 

large diurnal variations in temperature. Subsequently, two of the meas- 
urements have been repeated it times in the lunation when the rate of 
temperature change at the heater location was minimal. In addition, 
diurnal temperature variation from preceding and succeeding lunations 
are available to help interpolate trends during the time that the heater 
is on. lastly, some refinements have been made in che finite difference 
models of the conductivity experiments. The newer models indicate small 
adjustments in our previously published values. The revised Apollo 15 
conductivity values at each heater location are given in Table 6. 
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Heat Flow 

The magnitude of the vertical component of heat flow in the regolith 
can be calculated from the temperature and conductivity profiles in 
figure 5. Over each depth interval z^'to z % 


F = K 

Y z *a ve 1 2 


(Eq. 10) 


V 2 1 
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where AT, , is the corrected temperature difference listed in Tables 2 

z l~ z 2 

and 4, and K ave is the average conductivity in the depth interval 
calculated from the layered models in figure 5. 

Gradients, average conductivities and heat flows calculated from 

Apollo 15 and Apc^To 17 results are presented in Table 7. Notice tha' the 

bottom line of Table 7 for Apollo 17 probes if 1 and #2 presents heat flow 

data over the entire depth range of temperature measurement. At probe #1 

the most representative value of heat flow is thought to be that determined 
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by the probe data, 2. 8 x 10 W/cm . At probe H2 the measurement is pos- 
sibly disturbed as will be discussed below and the most representative 
value is that calculated using data between 67 and 234 cm, 2.5 x 10 W/cm . 
The heat flow calculated over the interval 91-138 cm is felt to be the best 
value from the Apollo 15 measurements. 

At the Apollo 17 probe #1 site the heat flow is quite uniform over the 

entire depth range. The variation falls well within the estimated error of 

measurement. Probe #2 results show a uniform heat flow along the length of 

the probe but heat flow between 67 and 131 cm is 70% greater. The large 

change in gradient is only partly compensated for by an increase in con- 

ductivi ty with depth. The overall heat flow, 2.5 x 10 W/cm , is in fair 
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agreement with the probe #1 value 2.8 x 10 W/cm. 

DISCUSSION OF HEAT FLOW RESULTS 

The Probe #2 f'p a surement 

The change in heat flow at probe #2 by a factor slightly less than 2 
over the depth range 67 to 234 cm is most reasonably explained by refrac- 
tion of heat flow in the vicinity of a large buried boulder. A 
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relatively large number of rocks are strewn over the ALSEP area. Lunar 
basalts have conductivities of about 1.2 to 1.8 x 10" 2 W/cra-K at 250K 
(ref. 15). These values are 60 to 90 times the conductivity of the fine- 
grained regolith materia!, thus large blocks of solid rock in the sub- 
surface can result in significant shunting of heat flow. 

To illustrate the shunting effect, the distortion of heat flow lines 
and isotherms around and through a square of material having 60X the con- 
ductivity of a surrounding infinite medium is shown in figure 7. The 
model is t..o-dimensional and symmetric about the left-hand margin of the 
figure. One significant feature of the model is that very little effect 
is felt at distances > 1/2 the width of the rock. Thus heat flow measure- 
ments would have to be made quite close to a rock, < 1/4 the width, to 
detect a disturbance as large as that at probe #2. Or; the other hand, 
probe #2 must be more than 5 cm from a rock in order for the rock not to 
have a detectable effect on the heater-activated conductivity experiments. 

If probe £2 were located relative to a large subsurface boulder in 
a zone defined by the dashed rectangle in figure 7, a temperature profile 
similar to that observed would result. Other features such as the slight- 
ly higher absolute temperatures measured by probe #2 relative to probe #1 
at all depths would also be explained. It should be noted that to ex- 
plain the rather sudden change in gradient shown in figure 5, a rather 
sharp cornered rock is required. 

If the vertical gradient were calculated using the temperature diff- 
erence between poi (a) and (b) in figure 7, the result would not be 
much different from the vertical gradient far from the rock. This is be- 
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cause the measurement interval is large enough to span the zone of dila- 
tion of isotherms adjacent to the rock and the compensating zone of com- 
pression of isotherms above it. If the same is true for prooe #2, then 
the gradient determined by using temperature prints at 67 and 234 cm 
should be close to the regional value. The fact that the average gradi- 
ent at probe #2 is in close agreement with results of probe #1, which is 
thought to be free of such disturbances, supports this possibility. The 
results of probe #1 should be regarded as the more representative value 
of heat flow at Taurus Litcrow. 


The Effects of T opography : 

Detailed assessment of the effects of topography on the Taurus 
Littrow measurement have not been completed in time for this report. How- 
ever, some idea of the magnitude of expected disturbances can be made 
based on very simplified terrain models. On the Moon there are two im- 
portant effects of topography: 1) Surface relief causes a distortion of 

subsurface isotherms to conform to the surface temperature distribution. 
Consequently a lowering of regional heat flow is generally found over 
features that are convex up and an augmentation of heat flow over concave 
up features. 2) Because the mean surface temperature depends on a 
radiative balance of solar flux, radiation to space and reradiation from 
the surrounding lunar surface, lunar topography has a significant effect 
on the surface temperature distribution. In general, depressions such as 

craters will have slightly elevated mean temperature because of the de- 
creased view of the crater floor to space and reradiation from other parts 
of the crater (ref. lg). On the other hand a north-facing slope such as 
that of the South Massif would have a lower mean temperature because of 
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the relatively large angle of incidence of solar rays throughout the lunar 
day. 

Craters of all sizes greater than 1 m in diameter that are within a 
distance of one crater diameter could have a measurable effect on the 
heat f 1 ow measurement. For example, the small craters near probes #1 
and #2 must be considered as well as the 600 m crater called 

Camelot about 600 m to the east of the ALSEP site. Most of these craters 
have an aspect ratio, diameter over depth, of about 6 to 1. The dominant 
effect of such craters is to increase heat flow in areas just outside the rim 
due to the slightly higher mean surface temperature insid". Finite difference 

models show that the excess heat flow falls off very rapidly with distance from 
the crater rim. At one crater radius from the rim of a crater whose floor 

Is maimer by 3*C than the surrounding flat surface, heat flow is increased 
- .1 W/cm 2 . Since the heat flow probes are one or more radii away from 
the rims of all craters of interest, the combined effect of all craters 
would be small. An estimate of .3 W/cm 2 or 10% would be conservative. 

The important result is that a negative correction must be applied to the 
observed values to compensate for nearby craters. 

The valley floor at Taurus Littrow aside from craters is relatively 
smooth and only the mountain ranges north and south of the site would 
have significant effects. The effect of the massifs can be estimated 
using a method developed by lachenbruch (ref. 17). The valley at Taurus 
Littrow is modeled as a flat-floored trough in an otherwise flat surface. 

The walls of the trough have uniform slopes equal to the average slopes of 
the north and south massifs measured from topographic maps, lachenbruch 
has published tables which permit an estimate of the effect on heat flow of 
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such sloping surfaces intersecting horizontal planes. Based on this model the 
topographic effect of the massifs will cause a +20% increase i . heat flow 
at the ALSEP site. This estimate is maximal since it assumes teat the trend 
of the ridges extends infinitely at their maximum elevation. The effect 
of the South Massif in particular is overestimated because the ridge does not 
extend very far to the southeast of the ALSEP site. 

Lastly, the valley at Taurus Littrow would have a slightly greater 
mean temperature (- 1 - 2K) than the surrounding regions because it behaves 
radiatively like a crater. The flow of heat from the warmer valley floor 
to the surface outside the bounding massifs would tend to decrease the 
heat flow in the valley. This effect has not been quantitatively estimated 
at this time. It may be significant and would tend to counterbalance the 
positive effects already described. 

In summary, our preliminary analysis of the effect of topography on 
the Apollo 17 measurement indicates that a negative correction estimated 
at 15 to 25%, should be applied to this observation. However, we icel 
that we should await a more detailed and careful analysis before assuming 
that effects of topography result in a significant difference between the 
heat flow measurements at the Taurus Littrow and Hadley Rille sites. 

Possible correlation with surface radioactivity 

As described above, a correction for the gross effects of topography 
of about -20% might be applicable to observed heat flow at Taurus Littrow, 
which would result in a value about 25% less than that measured at Hadley 
Rille. Because the difference depends on the topographic correction, its 
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significance may be questionable. However, taken at face value there is 
a correlation between the heat flow at the two sites and the surface radio- 
activity, as measured by the Apollo 15 gamma ray spectrometer (ref. 18). 
Relatively high thorium abundances were observed over the southeast corner 
of Mare Imbrium and appreciably smaller abundances over the southeast 
corner nf Serenitatis. The gamma ray spectrometer is sensitive only to 
isotopes in the most superficial layer. This te r+ 'tive correlation with 
the heat flow measurement, which detects the integrated effects of radio- 
genic heating to depths on the order of - 300 km, indicates the surficial 
variations may extend to depth. It further suggests that the variation 
of surface radioactivity may be the best available indicator of the varia- 
tion of heat flow over the lunar surface. We want to reemphasize these 
conclusions must await a more thorough analysis of topographic effects. 

Comparison with Earth-Based Microwave Measurements 

Prior to the in situ measurements of Apollo 15 and 17, estimates of 
the lunar heat flow have been made based on centimeter wavelength obser- 
vations of the Moon’s natural emission. (See, for example, references 
7 - g}. As is well-known, beyond the infrared, the 

lunar regolith becomes increasingly transparent at longer wavelengths. 

Thus, observations at longer wavelengths yield temperature information at 
increasing depths into the lunar regolith. For wavelengths greater than 

about 5 cm, the effective emitting layer is far enough below the lunar 
surface that no variation in apparent brightness temperature over a luna- 
tion can be detected. However, an increase in brightness temperature 
with increasing wavelength has been observed. A plot of these measure- 
ments is given in ref. 9. Between 5 and 20 cm, a nearly linear increase 
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of brightness temperature, T^(\)» wi th wavelength, A, yields an average 
spectral gradient of aT^/?* = O.SK/cm. To interpret the spectral tempera- 
ture gradient in terms of a heat flow, one must estimate the electrical 
properties, as well as the thermal conductivity, k, of the effective 
emitting layer for the 5 - 20 cm waves. In particular, the characteristic 
penetration depth of an electromagnetic wave, t e (A), must be known as a 
function of wavelength. 

Tikhonova and Troitskii (ref. 9 ) used simplified models of thermal and 
electrical property profiles to explain the microwave data over the spec- 
tral range 5 - 50 cm. Their resultant inferred heat flow values of 
2.9 - 4.0 x 10~ 6 W/cni are in remarkable accord with the in situ Apollo 
measurements, considering the approximations and assumptions necessary to 
the remote determination, in particular, the electrical and thermal para- 
meters were estimated from observations only applicable to depths char- 
acterized by diurnal temperature variations. The assumption that near 
surface parameter values apply to meter depths could lead to significant 
errors in the interpretation of the observed spectral gradient. Electri- 
cal property determinations based on remote radar measurements as well as 
measurements on returned lunar samples must be used to interpret the 
measured spectral gradient in terms of a thermal gradient. If the elec- 
trical and thermal properties of the reoolith are considered homogeneous 
for depths greater than about 70 cm, and extend at least 5 - 10 m, the 

thermal temperature gradient can be expressed in terns of the spectral 
gradient by the following equation: 


aT X 3 V 3X 

vrw; 


J 


3 

1 



(Eq. 11) 



where R is the appropriate centimeter wave reflection coefficient for 
the lunar surface - space interface. 

For very low electrical conductivities as are found in the lunar 
regolith, the electromagnetic penetration depth, £ e (x), may be written 
for the centimeter wave spectral region as: 

£ e (x) = tan a) (Eq. 12} 

where e is the dielectric constant and tan A the loss tangent at centi- 
meter wavelengths. 

The average temperature gradient measured in situ at the Apollo 15 
and 17 sites, 0.017 K/m, would produce the observed spectral gradient if 
/e tan A i .003, assuming R = 0.05. The feasibility of such a value for 
/c tan A is supported by direct surface observations of Weaver in the 
0.4 to 3 cm wavelength range (ref. 19). Direct measurements of returned 
Apollo samples over a wide range of frequencies indicate a dielectric 
constant for the regolith material in the range 2.2 - 3.2 that is nearly 
frequency independent. However, loss tangent measurements yield values 
in the range 0.0004 - 0.01, and are frequency dependent (ref. 20, 

21. 22). Additional electrical property measurements and refined analysis 
of the existing data on regolith samples must be made before the thermal 

gradient measured in situ can be supported on a Moon-wide basis by the 

spectral gradient observations. 



The Representati venes s of the Two Heat Flo w Measur em ents 

The regional geological setting; of Hadley Rille and Taurus Littrow 
are quite similar. Both are located on lava-flooded embayments at the 
edge of mascon basins. If the heat flow is influenced by structural or 
compositional anomalies unique, to this type of region, they would affect 
both measurements. To that extent they would not be representative of 
global flux from the moon. On the ether hand, the possible compatibility 
of our results with the microwave emission spectral gradient between 5 
and 20 cm wavelength lends support to the possibility that local anomalies 
at the two sites are not large. 

Despite the reservations set down in the previous paragraphs, the 
existing data concerning heat flow from the lunar interior indicates that 
a significant area of the Moon is characterized by a flux of between 2.5 and 

3.0 pw/cm 2 . Numerous thermal history calculations have shown that the 
contribution of initial heat (e.g. that gained during accretion) to the 
present surface flux is relatively small (ref. 10,23 ) even if the Koon 
were molten throughout initially. Some scientists have suggested that 
at the present time the Moon is thermally atsteady state (e.g. ref . 24 ). 

In either case, it follows that a predominant part of the surface flux, 

2.0 to 3.0 pw/cm 2 , must result from radioactive isotopes in the Moon. 

The geochemical data is convincing that most of these isotopes are con- 
centrated in the outer layers of the Moon. In addition, the abundances 

indicated by the heat flow values would require that the heat sources be 
near the surface to prevent melting in the outer several hundred kilo- 
meters. 

The radiogenic heat production per cm 3 of rock can be expressed in 
terms of the abundance of uranium since the ratios of the other two im- 
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portant long-lived, heat-generating isotopes, K and Th, to uranium is 
well established and quite uniform in the lunar samples. The heat production 
per unit volume at the present .time in W/cm 3 is about 0.71 times the 
uranium abundance in ppm (e.g. ref. 23 and 25). If most of the uranium 
is concentrated within 300 km of the surface, so that it contributes to 
the present flux, then the total lunar uranium abundances required to 
contribute 2.0 to 3.0 yW/cm 2 to the heat flow is about .05 to .075 ppm. 

These abundances are much higher than those of chondrites and significantly 
higher than estimates of the Earth's total abundance of about .03 ppm 
(ref. 4). 

SURFACE TEMPERATURES DEDUCED FROM THERMOCOUPLE MEASUREMENTS 
At each of the two heat flow holes, one of the thermocouples is em- 
bedded in a section of the cable which is approximately 15 cm from the 
top of the borestem and suspended above the lunar surface (see fig. 8). 

These thermocouples are in radiative balance with the lunar surface, the 
solar flux and space and hence provide a measurement of the surface bright- 
ness temperature throughout the lunation. The flux balance equation gov- 
erning the thermocouple temperature is: 

2TTa c dl e c <rT c * * 2TTa c c/J F c . m ©Ccr CT Turk 
+ 2. a c d l £ 

+ 2T racJl SAfCcs Fc-m C&5A5in</> (Eq. 13) 
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where: 
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= thermocouple temperature 
= lunar surface brightness temperature 
= radius of thermocouple cable 
= elemental length of cable 
= infrared emissivity of cable 
= infrared absorptivity of cable 
= absorptivity of cable to solar flux 
= view factor of cable to the lunar surface, including 
the surrounding mountains 

= cosine of the angle between the sun line and the cable 
axis, a function of cable orientation and lunar phase 
angle 


k = selenographic latitude at Taurus littrow 
I* = lunar phase angle measured from local sunrise 
= infrared emissivity of the lunar surface ^t.0 
6 = the mean solar constant = .1352 W/cm 2 
A = the lunar albedo = 0.08 
O’ = the Stefan-Bol tzmann constant 


The first term on the right-hand side of equation 13 represents 
flux into the cable element from the lunar surface; the second term 
represents direct flux from the sun; the third term represents solar 
energy reflected diffusely from the lunar surface and impinging on the 
cable. 
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The radiative properties of the cable, C c > °^cir > rtf cs were 
determined by laboratory measurement prior to the Apollo 17 mission. 
The cable orientations for both probe locations were determined from 
ALSEP photographs. 


Solving equation 13 for the surface brightness temperature yields: 


VTF — 5, A Fc-r* cosA SlM <?> 




During the lunar night, using and £ t s «^ C ir» equation 1$ reduces to 


4 -ifc 


=[£] 


(Eq. 15) 


Equations 14 and 15 assume that the surrounding mountains at Taurus 
Littrow are at the same temperature as the surface throughout the luna- 
tion. The deviation from this assumption, especially during the lunar 

day, may be quite large. However, both of the thermocouples hi.ve view 
factors to the mountains about 1/12 of their view to the surface. Thus, 
even large anomalous temperatures on the slopes of surrounding mountains 
will produce only negligible errors in the surface temperature deteimi na- 
tion. 

Shown in figure 9 is a full lunation plot of deduced lunar surface 
brightness temperatures at Taurus Littrow. Vertical bars represent 
estimated error bounds. The daytime temperatures were determined solely 
from the temperature data of the exposed probe #2 thermocouple for two 
reasons. First, the orientation of the thermocouple at probe #1 was 
much more difficult to obtain from the photographs; second, the probe ii 1 



thermocouple appeared to have a more substantial view of the radiation 
shield atop the bores tem; the radiation shield is a highly reflective 

square of aluminized mylar which could add a substantial unknown factor 

/** 

to equation 14. 

The nighttime temperatures are not subject to the errors due to 
uncertainties of the cable orientations and the data shown is an aver- 
age of the two thermocouple reductions. Nighttime surface temperatures 
deduced from each of the two above-surface thermocouples differed by no 
more than 2K throughout the lunar night. The dat a gaps near sunset and 
immediately following sunrise correspond to times of rapid temperature 
changes. During uhese periods, equation 14 loses its validity as the 
finite time constant of the thermocouple cable must be taken into 
account. 

From the data of figure 9, a mean surface temperature of 216K 

(+5K) was calculated for the Taurus Littrow site, indicating a mean 
temperature rise of - 39K between the surface and the top probe sensors. 

Only a small part of this mean temperature rise (no more than 5K ) can 

be accounted for by the measured heat flow. The mean temperature gradient 

is due mainly to the contribution of radiative heat transfer within the 

highly porous dust layer about 2 - 3 cm thick at the surface. During 

the warm lunar day heat is transferred more effectively into the lunar 

surface than it can be transferred out during the cold lunar night. 

To conserve net flux over a lunation, a mean temperature gradient is 

established, mainly confined to the porous surface layer. A similar 

phenomenon was observed at the Apollo 15 site, where a mean temperature 

rise of 45K was measured between the surface and the top probe sensor. 
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(The value of 35K reported earlier (ref. 10 has been revised to include the 
effect of the occlusion of the early morning sun at Hadley Rille due to 
local topography as well as refinements in tne thcrmocouple-to-surface 
temperature reduction calculation.) The total rise in mean temperature 
actually takes place almost entirely over the first few centimeters. 

The temperature dependence of the diurnal heat transfer in the 
near surface layer can be examined quantitatively by postulating an 
effective thermal conductivity which is a function of temperature: 

k k c 4 K ‘ T J (Eq. 1 6) 

equation 16 

The functional form of / has been verified experimentally for silicate 
powders in vacuum (ref. 26 ). Cremers and Birkebak have also found the 
conductivity of returned lunar fines tr fit a functional relationship of 
the form ofeqiatfon 16(.»Bf.27). The parameter R 350 which equals l< r *350 3 /K c 
first used by Linsky in examining this phenomenon (ref. 28 ), represents 
a measure of the radiative contribution to the heat transfer in the 
upper few centimeters. By using one-dimensional models of the lunar 
regolith, one finds that must be within the range 1.7 - 2.2 for 
the Apollo 17 site and between 2.5 and 3.0 for the Apollo 15 site to 
produce the observed mean temperature gradients. The range of values 
for the Apollo 17 site correspond closely to the value of 1.48 
obtained by Cremers and Birkebak for returned Apollo 12 samples (ref. 27 ). 
Similar measurements on Apollo 11 samples yielded R 3 £q values approximately 
equal to 0.5. It is important to note, however, that even very slight 
disturbances to the in situ configuration of the porous lunar fines 
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may have large effects on their thermal properties. In any case, the 
fact that a large mean temperature gradient has been observed at both 
lunar heat flow sites, separated by 700 km , strongly indicates that 
a porous layer in which radiative heat transfer plays a dominant role 
is a prevalent feature of the lunar surface, at least in the mare 
regions. 

It is fortunate that our more accurate surface temperature determina- 
tions are made during the lunar night for it is the post-sunset surface 
cool-down data which is most strongly constrained by the thermal 
properties within the top 6 - 10 cm of the surface. Figure 10 shows 
on expanded scale the reduced lunar surface nighttime temperatures at 
the Apollo 15 and Apollo 17 sites. The solid curves represent best 
fitting computer models of the thermal properties of the upper 15 cm 

of regolith at each site. The density nd mean conductivity profiles 
used to produce the theoretical ^ .. curves are shown in the inset 

figure. (Mean conductivity is the effective cono'uctivi ty at a given 
depth at the mean temperature of that depth.) The density profile used 
for the Apollo 15 site is based on inferences drawn from drill core 
penetration rates and surface disturbance due to astronaut activity reported 
by J. Mitchell (personal communication) of the Apollo 15 soil mcch-mirs 
team. The density profile for the Apollo 17 site was determined from prelimi- 
nary examination of returned core tube samples (D. Carrier, personal communication). 
In both the Apollo 15 and 17 models the heat capacity as a function of 
temperature was taken from Robie ct al. (ref. 29) data. on returned Apollo 11 
samples. In both models, a low conductivity layer of about 2 cm 
thickness is required to fit the steep drop in surface temperature 
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immediately followiny sunset. The Apollo 15 model then requires a steep, 
hut not discontinuous, rise in conductivity with depth down to 5 cm to 
produce the increased flattening of the cool-down curve through the lunar 
night. The Apollo 17 model, however, requires a very sharp rise in 
conductivity at a depth of about 2 cm to produce the abrupt flattening 
of the cool-down curve at about 190° phase angle. The subsr .ent in- 
crease in conductivity with depth is slight enough so that the Apollo 17 
model may be considered essentially a 2-layer model. The large jump in 
conductivity at 2 cm is also supported by the preliminary density pro- 
file which indicates a fairly high density quite close to the surface. 

The Apollo 15 density prorile, on the other hand, supports the possibility 
that a substantial conductivity gradient exists over the upper 30 cm 
of the regolith. 

It is important to emphasize that the most critical surface tempera- 
ture data required for the purpose of determining thermal regolith pro- 
files is that during the 10 - 30 hours immediately following sunset. 
Surface temperature data during this period have been the most difficult 
to obtain from remote infrared brightness scans. The level and steep- 
ness of the cool-down data immediately following sunset is controlled 
almost entirely by the thermal properties of the upper 2 cm. If the 
very early nighttime data is not sufficiently accurate to constrain the 
thermal properties of the upper 2 - 3 cm of dust layer within +30%, 
then subsequent attempts to determine deeper conductivity values unambigu- 
ously from the flattened part of the cool-down curve will not be possible. 
For example, the broken line curve of figure 10 fits the nighttime data 
after 192° phase well within the error bands of the data. However, 


discrepancies in the early post-sunset fit produced by different 
conductivities within the upper 2 cm lead to discrepancies up to a factor 
of 2 in conductivity determinations i ' depths below 2 cm. (See curve 
17a in inset of fig. 10.) 


CONCLUSIONS 

With this last Apollo mission two heat flow measurements have been 

successfully made on the lunar surface. Both measurement sites are in 

similar regional settings in the northeast quadrant of the Moon. 

Taurus LHtrow and Hadley Rille sites are located in embayments in the 

i 

mountainous rims of the Imbrium and Serenotatis mascon basins which have 
been flooded by mare-type basalts. 

Surface brightness temperatures were calculated from the temperature o 

thermocouples suspended several / cen ^al)Ove S the lunar surface. The mean 

surface temperature at Hadley Rille throughout a lunation cycle is 207 K. 

The mean temperature increases with depth very rapidly in the upper few 

centiireters and is about 252 K at a depth of 90 cm. The main reason for 

this increase of 45 K is the predominant role of radiative heat transfer 

in the loosely packed upper layer. During the lunar night the surface 

temperature at Hadley Rille falls to 93 K. Fran the cool-down history 

after sunset we have deduced that the upper 2 cm of the regolith is char- 
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acterized by a conductivity of 1.5 x 10 W/cm-°K. Below this depth the 
conductivity increases rapidly and probably, in a continuous manner until 
it reaches values of about 1.5 x 10 _ * W/cm°K at depths where the probes 
are emplaced. 


At Taurus Littrow the mean surface temperature is 216 K and, as in the 

case of Hadley Rille, increases a few tens of degrees in the upper 2 cm so 

that at a depth of 67 cm a mean temperature of 254 K is measured. The 

minimum temperature just before lunar dawn is 103 K, 10°K higher than Hadley 

Rille. This higher temperature is primarily due to the existence of a 

relatively high conductivity layer at a depth 2 cm below the surface. From 

the point of view of thermal properties the regolith at Taurus Littrow can 

be described as two layers: an upper 2 cm loosely-packed layer of very low 

conductivity( 1.5 x 10"^ W/cm-K) in which heat transfer by radiation pre- 

-4 

dominates and a lower layer with much higher conductivity (>1.2 x 10 W/cn-K) 

, 3, 

and higher density (1.8 - 2.0 am/cm }. 

Subsurface temperature and conductivity measurements at depths below 

SG cm, where the large diurnal variations are negiiqibly small, indicate 

a steady state heat flow through the surface at Hadley Rille of 
•6 2 -6 2 

3.1 x 10 W/cm and 2.8 x 10 W/cm at Taurus Littrow with an estimated 
error of +20i. These fluxes are deduced from average temperature gradients 
in the regolith between 1.3 and 1.7 K/m and an average conductivity in the 
range 1.7 x 10 - ^ to 2.0 x 10"^ W/cm-K. Conductivity generally increases 
with depth in the regolith although some layering, with high conductivity 
materials overlying lower conductivity materials, is found at both sites 
(see Table 6). A conductivity value of nearly 3 x 10^ W/cm-K was measured 
at the bottom of probe #1 at the Apollo 17 site. Thermal gradients de- 
crease with depth in some cases in response to the increase in conductivity. 

At Taurus Littrow probe site $ 2 a large decrease in gradient with depth is 
possibly attributable to a large subsurface boulder in close proximity to 
the probe. 
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The heat flews at both sites are to some extent affected by local 
topography. Preliminary estimates indicate a correction of -15 to -25'/ 
may be applicable to the Taurus Littrow values as a result of the adjacent 
massifs. However, a more refined analysis is required. 

The heat flow measured at the two sites is approximately 1/2 of the 

-6 2 

Earth's average heat flow (6.3 x 1C W/cm ). If these two values are 
representative of the Moon as a whole, then a heat flow of 1/2 the 
Earth's requires a heat production per unit mass 

for the Moon's interior over twice that of the Earth. This statement assumes 

both planetary bodies are near steady state so that total surface heat loss 

is nearly equal to the present interior heat production. 

40 235 238 232 

Since the long-lived radio isotopes of K, U, U and Th are 
the principal source of heat in the Earth and Moon, the heat flew results 
imply a two to three fold enrichment of uranium in the Moon relative to the 
Earth. Lunar samples show that the abundance of potassium relative to uranium 
is 1/3 to 1/4 that of the Earth so that in the Moon uranium is the main 
contributor to internal heating. These isotopes must be concentrated in the 
outer 100 to 200 km of the Moon to avoid extensive melting at shallow depth 

t 

at present. 

Re-interpretation of Earth-based measurements of microwave brightness 
temperatures using the new data on regolith thermal and electrical properties 
will be important in determining the representativeness of the in situ lunar 
heat flow measurements. Further more-refined microwave observations of the 
Moon, especially narrower beamed measurements over discrete portions of the 
lunar disc, would be valuable in determining possible variations of heat 
flow over the lunar surface. 
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TABLE APOLLO 1? HIE SUBSURFACE TEMi'I TUTURE DATA 
TEMPERATURE Dl FEE PENCE ME AS U RKMFN TS 
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BRIDGE 

INTERVAL 

(cm) 

EQUIL. TEMP. DIET . 1 
(K) 

CORK. TEMP. DIFF. 2 

U) 

ANNUAL EAVT 
CORRECTION 

DTG11 

130-177 

0.707 

0.755 

• -0.027 

DTR11 

139-168 

0.435 

0.467 

- 0.018 

DTG12 

185-233 

0.533 

0.559 

- 0.001 

DTR12 

194-224 

0.322 

0.326 

<.001 

DTG21 

131-178 

0.370 

0.390 

->.027 

DTR21 

l4 0-169 

0.218 

0.223 

- 0.018 

DTG22 

186-234 

0.336 

0.359 

- 0.001 

DTR22 

195-225 

0.206 

0.212 

<.001 


1 The error associated with extrapolating to equilibrium temperature 
differences is ± . 003 K. 

2 The uncertainty introduced by these corrections is estimated to be 

± 2 $. 


| 
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ABSOLUTE TEMPERATURE MEASUREMENTS 


SENSOR 

TYPE 

DEPTH 

EQUIL. TEMP. 3 

CORR. TEMP, 



cm 

(K) 

(K) 

thermocouple tca 

66 

254.2 

254.2 

PT.ATTJ'TITT.'i 

PPSTST. 

130 

255.06 

255.02 

If 

It 

139 

255.19 

255J7 

tt 

1 ! 

168 

255.62 

255*64 

It 

II 

177 

255.76 

255.70 

It 

II 

185 

255.91 

255.91 

tt 

It 

194 

256.03 

256.04 

It 

11 

22-. 

256.36 

256.37 

11 

II 

233 

256.44 

256.47 

THERMOCOUPLE TCA 

67 

254.7 

254.7 

PLATINUM 

RESIST. 

131 

256.07 

256.05 

tt 

11 

140 

256.09 

256.09 

II 

ti 

169 

256.31 

256.31 

tt 

11 

178 

256.44 

256.44 

II 

11 

186 

256.48 

256.48 

II 

« 

195 

256.52 

256.53 


ii 

225 

256.73 

256.73 

•» 

•1 

234 

256.8? 

256.84 


3 The accuracy of extrapolated absolute temperatures is ±0.05 K. 

4 The correction for the annual wave to be applied to the thermocouple 
is + 0.04 K. 
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TABLE 3: THERMOCOUPLE ERRORS AND CORRECTIONS 


ERROR (TCI - TOP GRADIENT SENSOR), K 


SENSOR 

LUNAR DAY 

LUNAR NIGHT 

TC11 

+1.1 

+1.9 

TC21 

+1.4 

+2.5 


ESTIMATED CORRECTIONS* 

SENSOR 

LUNAR DAY 

LUNAR NIGHT 

TC14 

+0.4 

+0.7 

TC24 

+0.6 

+1.2 


■k 

The uncertainty in determining these corrections 


is estimated as +0.4°K. 



TABLE 4: AFOLLO 15 HFE SUBSURFACE TEMPERATURE DATA 
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TABLE 5: CONDUCTIVITIES FROM COOL-DOWN HISTORIES 

CO' ’'’JCTlVm (W x lO“‘7cm-K) 

SENSOR DEPTH HEATER LOCATION* WITH DRILL WITHOUT DRILL 

HEATING EFFECTS HEATING EFFECTS 

PROBE 1 


66 


1.0 

- 


130 

' Hll 

2.3 

1.3 

1 

139 


1.9 

1.1 


168 


1.9 

1.1 < 


177 

H12 

2.0 

1.0 

* 

185 

H13 

1.9 

1.1 


194 


2.1 

1.1 

* 

224 


2.8 

1.4 


233 

H14 

2.7 

1.6 




PROBE 2 



67 


1.0 

- 

{ 

131 

H21 

2.0 

1.2 

t 

140 


2.0 

1.1 


169 


2.4 

1.3 

J 

178 

H22 

2. 7 

1.7 

\ 

| 

186 

H23 • 

2.9 

1.7 

i 

i 

195 


2.7 

1.5 

V 

■> 

i 

i 

225 


2.3 

1.5 

j 

* 

234 

H24 

2.5 

1.5 




TABLE 6: RESULTS OF liFE CONDUCTIVITY EXPERIMENTS 


APOLLO 17 


HEATER 

DEPTH 

CONDUCTIVITY 1 


CONTACT CONDUCTANCE 

LOCATION 

(cm) 

(W x 10-'7cm-K) 


(W x 10'Vc.m 2 -K) 



PROBE 

1 


Hll 

130 

2.50 


1.4 

H12 

177 

1.72 


1.6 

H13 

185 

1.79 


1.4 

H14 

233 

2.95 


1.2 



PROBE 

2 


H21 

131 

2.06 


1.6 

H22 

178 

2.36 


1.1 

H23 

186 

2.64 


1.5 

H24 

234 

2.24 


2.3 




APOLLO 

15 


HEATER 

DEPTH 

CONDUCTIVITY 


CONTACT CONDUCTANCE 

LOCATION 

(cm) 

• (tf x 10 -t */cm-K.) 


(W x 10 -,, /cm 2 -K) 



PROBE 

1 


Hll 

35 

1.41 


0.8 

H12 

83 

2.11 


0.8 

H13 

91 

1.60 


0.9 

H14 

138 

2.50 


1.0 



PROBE 

2 


H23 

49 

1.46 


0.5 

H24 3 

96 

2.43 


0.6 


1 The estimated error of conductivity measurement is +15%. 


2 Estimated error is +20%. In the theoretics’ model the thickness 
of the contact zone is 2 mm. 

3 It is very probable that a section of broken borestem lies jus'; 
outside this location so that the uncertainty of this measure- 
ment is very large. 


TABLE 7: 


HEAT FLOW DATA 


DEPTH 

INTERVAL 

(cm) 


TEMPERATURE 

GRADIENT 

(K/cm) 


AVERAGE 

CONDUCTIVITY HEAT FLOW 

(W x 10~'*/cm-K) (W x 10 -e /cm 2 ) 


APOLLO 17 
PROBE 1 


66-130 

0.013 

1.6 

2.1 

130-177 

0.0158 

1.79 

2.83 

139-168 

0.0163 

1.72 

2.80 

185-233 

0.0118 

2.39 

2.81 

194-224 

0.0113 

2.48 

2.81 

66-224 

0.014 

1.8 

2.5 


PROBE 2 


67-131 

0.021 

1.5 

3.1 

131-178 

0.0082 

2.26 

1.86 

140-169 

0.0078 

2.30 

1 79 

186-234 

0.0076 

2.50 

1.89 

195-225 

0.C074 

2.53 

1.87 

67-225 

0.013 

2.0 

2.5 


APOLLO 15 
PROBE 1 


91-138 

100-129 


0.0175 1.78 3.11 

0.0166 1.68 2.82 
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f&r****tok* w ** w. * . - - . . . 

FIGURE CAPTIONS 

1. At top is a drawing of a heat flow probe. At left bottom is a 
simplified schematic of a typical platinum resistance bridge circuit 
and on the right the thermocouple circuit. 

2. The geometry of the probe, borestem, contact zone and lunar regolith 
in the vicinity of a conductivity experiment. The dashed lines show 
surfaces of equal temperature rise in degrees K after the heater has 
been on for 36 hours. The model parameters are Km = 2.4 x 10“** W/cm-K 
and H 2 = 1.5 x 10' 4 W/cm 2 -K. 

3. The attenuation of the peak-to-peak amplitude of the diurnal and 
annual temperature with depth in the regolith is shown at the left. 

The 6 ' s are the effective decay constants deeper than 20 cm. At the 
right the phase lag with depth is shown. The values of * are the 
wavelengths of the thermal wave below 20 cm. The model used for the 
diurnal variation is from Apollo 17 data. See the last section of 
this paper. The annual curves are calculated from a Hadley Rille 
thermal properties model (ref. 14). 

4. Temperature histories of all sensors 65 cm or deeper at probes #1 and 
#2. The short pulses appearing on some of the sensor traces result 
from heater initiation for conductivity experiments. The numbers on 
each curve refer to the depths below the surface. The lowermost 
curves on each plot are TC4 thermocouple measurements. Some representa- 
tive data points from the thermocouples are shown to indicate the 
scatter of these measurements. Temperatures shewn are calculated by sub- 
tracting TCI values from those of TC4 and adding them to the top gradient 
sensor temperature. Corrections given in Table 3 have not been added. 
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5. Equilibrium temperatures, conductivities and heat flows measured by 
the Apollo 17 probes. The open circles on the conductivity plot are 
calculated from cool-down curves assuming maximum drilling energy 
and the solid circles are heater-activated experiment results. The 
solid line represents a layered model used for calculating heat flow. 
In the heat flow figure the solid lines give heat flow over the 
three largest intervals. The dashed line is the average heat flow 
between 66 and 234 cm. At the far left the geometry of the probe in 
the subsurface is shown. 

6. Temperature rise during a conductivity experiment (solid dots) is 
compared with a theoretical curve derived from a model with k m = 

2.7 x 10~ 4 W/cm-K and H 2 = 1.5 x 10~ 4 W/cm 2 -K. In the inset the 
temperature rise for times greater than lPM minutes is shown on an 
expanded scale plotted against the logarithm of time. The observed 
data is compared with two bracketing theoretical curves. The reduced 
conductivity is 2.64 W/cm-K. 

7. The effect of a square of material on vertical heat flow (shaded 
area) which has a conductivity 60 times that of the surrounding 
material, shown by the distortion of isotherms and flow lines. 

These results are based on a finite difference model computation. 

8. Photograph of probe #2 borestem protruding from the lunar surface. 

The Heat Flow Experiment housing is in the background. The thermo- 
couple is in the black portion of cable about 10 cm from the top of 
the stem. 

9. Surface temperatures determined from thermocouple measur 
Vertical bars are estimates of the error limits. 
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10. Temperatures during the lunar night at Apello 17 (solid dots) and 
Apollo 15 (open circles). Vertical bars are estimated errors. 

The continuous curves are theoretical curves derived from the 
thermal property models shov.-n in the inset. 
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